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A MORE ACCURATE AND MORE EXTENDED COSMIC-RAY 


IONIZATION-DEPTH CURVE, AND THE PRESENT 
EVIDENCE FOR ATOM-BUILDING 


By RosBert A, MILLIKAN AND G. HARVEY CAMERON 
CALIFORNIA INSTITUTE 


(Received December 9, 1930) 


ABSTRACT 

The cosmic-ray tonization-depth curve has been extended at both its upper and 
lower ends and made more accurate throughout. The absorption coefficients obtained 
directly from the slope of the curve run from »=0.35 per m. of water at the top 
(Pike’s Peak) to 1. =0.028 at the bottom (80 m. or 262 ft. of water below the top of the 
atmosphere, thus bringing to light both softer and harder components than the au- 
thors had before found. Strong quantitative evidence is presented, on the basis of the 
Klein-Nishina formula, that the strongest and most absorbable cosmic-ray band 
arises from the act of formation of helium out of hydrogen. Striking qualitative evi- 
dence is found that the three more penetrating bands are due to the formation out of 
hydrogen of the only other abundant elements oxygen (C, N, O) silicon (Mg, Al, Si, S) 
and iron (Iron group). Two independent proofs are given that the cosmic-rays enter 
the earth’s atmosphere as photons, namely, (1) they are quite uninfluenced by the 
earth's magnetic field, and (2) the ionization produced by them in a closed vessel does 
not increase continually in going to the top of the atmosphere but passes through a 
maximum, It is shown to follow that the cosmic rays, in coming from their place of 
origin to the earth have not passed through an amount of matter that is appreciable in 
comparison with the thickness of the earth’s atmosphere and that they must therefore 
originate in interstellar space rather than in the atmospheres of the stars. Some par- 
ticipation of the nucleus in the absorption of cosmic-rays is brought to light. 


1. OBJECTIVES 


HE new series of measurements presented herewith on the relation be- 
tween cosmic-ray ionization and depth in equivalent meters of water be- 
neath the surface of the atmosphere was undertaken for two very specific 
reasons. 
First—Our preceding experiments, published in full in 1928,' had brought 
to light what seemed to us very striking evidence that the cosmic rays have 
their origin in the acts of formation “in the depths of space” of the atoms of 
the celestially common elements helium, oxygen (C, N, O), and silicon (Na, 
Mg, Al, Si, S) out of hydrogen. This evidence consisted : 


1 Millikan and Cameron, Phys, Rev. 32, 533 (1928). 
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236 R. A. MILLIKAN AND G. H. CAMERON 


(a) In our experimental proof that within the limits of our observational 
uncertainty these rays show a uniformity of distribution, i.e., an independ- 
ence of both latitude and of sidereal time. Both of these conclusions had been 
established by our trip to the Bolivian High Andes? in 1926, the first point 
not having been previously tested at all by other observers, the second having 
been so tested but with opposite results, though subsequent more careful ex- 
perimenting by Hoffmann and Lindholm,’ Steinke,* Hess,® and by one of us® 
has confirmed our conclusion; 

(b) In our proof shown unambiguously by the curve itself, of the banded 
character of the rays; 

(c) In the general rough agreement between our observed sequence of 
band-absorption coethcients and intensities, and the sequence of energies re- 
leased, in accordance with Einstein’s equation and Aston’s curve, when the 
celestially most abundant elements helium, oxygen, and silicon are formed 
out of hydrogen; and 

(d) In the rough agreement between the observed absorption coetticients 
and those computed from Dirac’s formula connecting ray-energy with ray- 
absorption, 

If our interpretation of our cosmic-ray results is correct, rays of still 
higher penetrating power should exist corresponding to the formation of still 
heavier elements out of hydrogen. Iron, at least, is abundant enough so that 
our theory suggested that we might find a cosmic-ray band, or cosmic-ray 
bands, of higher penetrating power than any we had thus far definitely ob- 
served. Such rays could be brought to light only by working at still lower 
depths in snow-fed lakes with electroscopes still more sensitive than any we 
had thus far used, and this new series of experiments was started in the spring 
of 1928 in part to test this point. 

Second.— Our 1927-8 ionization-depth curve’ showed characteristics at 
its upper end, i.e., at the highest altitudes at which we had ionization-depth 
readings, which seemed to be a bit out of line with the theory. Thus, our 
analysis, by means of the Gold tables, of our curve into its “monochromatic” 
absorption coefficients yielded as the strongest and most absorbable com- 
ponent « =0.35 per meter of water, while the Dirac formula gave for the coeffi- 
cient of rays due to the formation of helium out of hydrogen np =0.30. While 
this was of the right order of magnitude, indeed quite close in view of our ac- 
curacy, the divergence was in the wrong direction, for our computated value 
0.30 had applied to a monochromatic beam, but the actual beam, even if it 
entered the atmosphere as monochromatic would, where observed, have sec- 
ondary, tertiary, etc., components, due to Compton encounters with elec- 
trons, and all such encounters tend, until the beam has got completely into 
equilibrium with its secondaries, to push down the observed uy, i.e., the uw ob- 


2 Millikan and Cameron, Phys. Rev. 31, 163 (1928). 

3 Lindholm, Gerlands Beitrage zur Geophysik 22, 141 (1929). 
* E. Steinke, Zeits, f. Physik 42, 570 (1927) and 48, 647 (1928). 
5 Hess and Mathias, Wien. Ber. 137, 327 (1928). 

6 R. A. Millikan, Phys. Rev. 36, 1595 (1930). 

7 Millikan and Cameron, Phys. Rev. 31, 925 (1928). 
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tained from the observed slope. The theory, then, shows no way by which 
the com puted absorption coethcient of the pure radiation due to the formation 
of helium out of hydrogen can be Jess than the observed coefficient. It rather 
requires that the computed uw be at least as large as, and if equilibrium has not 
vet been obtained at the observation point, appreciably larger than the ob- 
served uw. We thought that our observations at high altitudes showed some 
slight indications that our uppermost lake-readings were accidentally high 
and we therefore wished to repeat these observations with more sensitive 
instruments and under better conditions, and to extend them, if possible, 
to still higher altitudes, hoping that a crucial test of our theory might come 
out of such more accurate high-altitude readings. In a word, then, we under- 
took the present series of observations to extend and improve our observa- 
tional data at both the upper and the lower ends of the ionization-depth 
curve. 


2. TECHNIQUE 


Our mode of procedure was precisely the same as that used in obtaining 
the last ionization-depth curve published in 1928, save that in order both to 
bring out weak effects at greater depths under water, and to obtain increased 
precision in high altitude observing, we were obliged to increase still further 
the sensitivity of our electroscope. To do this we built a new spherical in- 
strument of steel, wall-thickness 3 mm, internal capacity 1622 cc, and filled 
it to a pressure of 30 atmospheres. (See Fig. 1.) This procedure, according 
to direct observational data to be presented later, multiplied our electroscope 
sensitivity 13.82 fold over that obtained under a pressure of one atmosphere. 
It represented, too, a sensitivity more than double that used in our published 
1928 observations.» We determined the capacity of this electroscope with 
much precision by the method heretofore described,* finding it to be 0.979 ab- 
solute electrostatic units. Although we have not reached with this instrument 
the extreme limit of possible electroscope sensitivity, since increases in vol- 
ume, slight decreases in fiber-capacity and still higher pressures are possible, 
vet, for the purposes for which it was to be used, this instrument came close 
to the limit of possible efficiency as a high frequency radiation detector. 

In order to test the insulating properties of our quartz supports, before fill- 
ing with air under pressure we pumped the chamber entirely free of air and 
found that it then leaked, to cite one particular test, from 226.5 volts to 226.0 
volts in four hours. The support-leak was thus not as much as one-tenth of 
the slowest rate of discharge ever observed in this work, and a small fraction 
of a percent of the average rate. Further, it disappears entirely in our com- 
putations on rate of ion-formation because it is included in what we call the 
zero of our electroscope. Such minuteness of the leak of the supports, how- 
ever, shows the completeness with which this source of error has been pushed 
into the background in these experiments. 

As in all of our cosmic-ray work we followed here the procedure of reduc- 
ing the fiber-deflection to volts at the instant at which it is read. This is done 


8 Millikan and Cameron, Phys. Rev. 31, 922—-5 (1928). 
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by calibrating the scale in the eyepiece just before or after each reading. This 
procedure eliminates completely all temperature effects on the fibers or on 
any of the electroscope parts. Thus in under-water work our procedure is to 
calibrate the scale for the deflection to be used under quiet reading-con- 
ditions af the surface, then to sink the electroscope to the desired depth and 
leave it for say twelve hours, then to bring it quickly to the surface, read and 
calibrate again. A small correction is applied for the discharge during lower- 
ing and raising by going through the operation of lowering and at once rais- 
ing and seeing whether an observable change is detectable, or by doing so re- 





Fig. 1. 


peatedly for the sake of magnification and then computing the correction for 
a single operation. These corrections are small and quite accurately obtain- 
able. In some of our earlier work we used self-registering electroscopes, but 
for this work we considered them of no advantage. 

The electroscope with which all of the results herein reported were ob- 
tained is shown in Fig. 1. It was pumped up to a pressure of 30.14 atmos- 
pheres, one atmosphere being reckoned from the conditions existing in the 
Norman Bridge Laboratory at the time of filling, namely, 24°C, 74 cm pres- 
sure, and this pressure has held with no trace of leak for now two years since 
the time of filling. 








————~ 
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The “electroscope constant” of this instrument in ions per cc per sec. is 


volts 0.979 107!" volts 
I cc/sec = —— X —- XK ———————————————. # 1.172 ———. 
hours 300 4.770 X 1622 K 36000 hours 

For observations on land this electroscope was provided with an accu- 
rately fitted, spherical lead shield 7.64 cm thick consisting of four layers of 
carefully shaped and fitted hemispherical shells, and two layers of shells 
divided into quadrants so as to fit closely the exterior electroscope wall, with 
its flange and bolts, and to facilitate the entrance of the necessary leading-in 
wires. The “water equivalent” of this lead shield as obtained by multiplying 
thickness by density was 85 cm. On account of the space occupied by the 
steel flange about which the inner layers of lead were fitted, the water equiva- 
lent was actually a trifle less than 85 em, but in any case by careful experi- 
ments with radium and thorium in about the proportions in which they occur 
in surface rocks the percentage of local rays getting through this lead shield 
was found to be 2.4 percent, and save for this small correction-factor, always 
determined and applied to all readings, the rays inside the lead were pure 
cosmic rays and could be compared, as shown below, with readings at the 
same equivalent depth in water beneath the top of the atmosphere. 


3. READINGS IN SNOW-FED LAKES 


The results reported in this section were taken partially during the sum- 
mer of 1928 and partially during the summer of 1929 in the two different Cali- 
fornia lakes 250 miles apart, Arrowhead Lake (altitude 5100 feet) and Gem 
Lake (altitude 9120 feet). Since we showed in our first work in snow-fed 
lakes in 1925 that the measured intensity of the cosmic rays is a single valued 
function of the depth of the superincumbant atmosphere, it of course follows 
that it depends upon the barometric pressure, as has since that time been 
many times noted. Accordingly, in Table I the barometer height is given for 
each observation, or group of observations. At the higher levels, i.e., down to 
about sixteen meters beneath the top of the atmosphere, where readings at 
the same depth beneath the surface of the lake are made at different baro- 
metric pressures, the results are reduced to a common pressure by applying 
a small correction, which is computed from the Gold tables,’ for the value of 
the absorption coefficient shown by the ionization-depth curve at the eleva- 
tion considered. This computation, however, checks nicely with the observed 
slope of the curve and this slope is of course independent of any theory, so 
that these small barometric corrections may properly be said to be purely em- 
pirical, and hence free from any other than observational uncertainties. This 
correction is negligible below sixteen meters but amounts at Gem Lake to 
slightly more than 1° and at Pikes Peak to about 2% per tenth-inch of mer- 
cury. The under-water readings, taken in two lakes and over the period from 
August 1928 to September 1929, are all collected in Table I. It will be seen 
from columns 7 and 9 that these readings extend from a level equivalent to 


® Gold, Proc. Roy. Soc. A82, 43 (1909). 
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8.25 meters of water beneath the top of the atmosphere down to a level of 80 
meters, or 262.5 feet, beneath the top, and that in that range of levels, or 
depths, the observed rate of discharge of the electroscope changes by more 
than thirty fold, namely, from a value represented by an intensity of ioniza- 
tion within the electroscope of 64.1 ions per cc per sec. down to 2 ions per cc 


TABLE I. Jonization-depth readings in snow-fed lakes. 





Depth 


Depth in below Ions per Mean 
Year Date Lake used Water Barometer Topof ce per I cc/sec 


Reading atmosphere second 


1928 Sept. 11 Gem .85 m 21.45 8.25 64.1 64.1 
“ “ 9 “ 1.00 “ 21.48 8.42 59.8) 
" " 9 1.00“ 21.48 8.42 60.4} 60.1 
“ “ 4 1.00“ 21.48 8.42 60.2) 
“ . th 2.00“ 21.46 9.41 43.8) 43.8 
“ 11 “ 2.00“ 24.41 9.40 43.8 wie 
1929 Sept. 6 Arrowhead aa“ 24.81 9.40 43.8 43.8 
“ “ 7 . 1.50“ 24.77 10.06 37.5 37.5 
“ “ 6 2.00“ 24.82 10.58 33.5 33.5 
1928 “= 12 Gem 4.00“ 21.65 11.48 30.7) 30.4 
. >. a“ 4.00“ 21.65 11.48 30.2) 
1929 * 4 Arrowhead 3.00 “ 24.94 11.59 29.6) 29.5 
. 5 . 3.00 * 24.94 11.59 20.4) sitions 
‘ “ 6 “ 4.00“ 24.79 12.56 25.5 25.5 
‘ “ 5 “ 5.00“ 24.87 13.56 23.0) 23.1 
‘ “ 6 5.00“ 24.87 13.56 23.2) _ 
“ “ 7 6.00“ 24.74 14.55 21.1 21.1 
1928 Aug. 22 6.25“ 24.92 14.86 20.2) 
‘ «23 6.25 24.93 14.86 20.8) 20.6 
‘ « 2 6.25 24.81 14.82 20.9 
. . 2B 8.25 24.92 16.85 17.10) 47 33 
“ . © 8.25 24.89 16.84 17.57) . 
. ._ 23 10.64 24.81 19.21 14.72) 44 5 
. ._ 3 10.64“ 24.89 19.24 14.33/ on 
. a. = ‘ 15.90 * 24.81 24.47 10.22) 10.23 
“ + & “ 15.90“ 24.83 24.49 10.24) - 
. ‘ 22 ‘ 21.10 * 24.91 29.70 7.81) > 
. . 23 “ 21.10 * 24.93 29.71 7.96) lie 
“ + 2 “ 26.25“ 24.81 34.82 6.07) eo 
“ « 1 “ 26.25“ 24.81 34.82 ry is 
“ ‘ 19 “ 30.35 24.79 38.91 5.34) 521 
“ ‘ 18 “ 30.35 24.80 38.91 5.08) ies 
“ ‘ 9 “ 37.05 24.86 45.64 4.33) 
“ ‘ 10 “ 37.05 24.86 45.64 4.09} 4.25 
“ + “ 37.05“ 24.86 45.64 4.34) 
. Sept. 9 Gem 43.00 21.48 50.42 3.68) 3.62 
“ “10 “ 43.00 21.39 50.41 3.56) hie 
" Aug. 17 Arrowhead 42.78 24.91 51.38 3.90) 
. we 45 “ 42.78 24.91 51.38 3.88 3.79 
“ >. 122 “ 42.78 ¢ 24.89 51.38 3.00 
. Sept. 8 Gem 50.00 “ 21.42 57.40 3.30 3.30 
“ - 22 “ 60.00 “ 21.36 67.38 2.49 2.49 
“ “10 “ 72.5 “ 21.39 79.90 1.95) 2.00 


. + 8 ‘ 72.6 “ 21.45 80.00 





| 
| 


per sec. The absorption coefficients as computed from the Gold tables at suc- 
cessive points along this curve are given in Table II. 

The curve starts a little higher up than does the 1928 curve, analyzed on 
page 927, Physical Review, Vol. 31, and it is significant that the absorption 
coefficient at the top is now a little higher than it was there, namely 0.27 in- 
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stead of 0.22, -1.e., it is larger, not smaller, then before. Also, in keeping with 
this fact the bend, or knee, at about 10 meters is sharper than ever, as shown 
by the change from 4 =0.27 to u=0.16 between 9.5 m and 10.5 m, when on 


TABLE II. Absorption coefficients at various depths, in meters of water, below top of atmosphere. 


























Depth in m. of water Absorption Depth in m. of water Absorption 

beneath top of atmos. coef. u beneath topofatmos. coef. u 
8.25- 9.5 0.27 20-30 0.045 
9.5 -10.5 0.16 30-40 0.038 
10 .5-11.5 0.11 50—40 0.028 
11.5-12.5 0.095 50-00 0.028 
12.5-15.0 0.067 60-80 0.028 
15 .0-20.0 


0.058 











the old curve we got u=0.20 between 9.5 and 10.5. Of course this means that 
the most absorbable cosmic-ray band springs into view from these figures more 
insistently than before. 


4. LAND-READINGS UP TO GREAT ALTITUDES 


The highest altitude snow-fed lake used for the foregoing readings was 
Gem Lake (altitude 9120 feet) and for the sake of being free from the possi- 
bility of effects due to the radioactive emanations of the atmosphere (though 
over large bodies of water these effects are actually very small) we used no 
reading nearer the surface of Gem Lake than 0.85 m, a level corresponding 
to 8.25 m of water beneath the top of the atmosphere. At this level, as in- 
dicated above, the curve was already beginning to show departures in the 
wrong direction for satisfactory explanation from the standpoint of the Dirac 
formula. However, in accordance with the second of the objectives discussed 
in $1, the most significant data were to be expected at still higher altitudes, 
and in order to obtain accurate data at least a meter of water higher up we 
arranged for a series of land observations as follows: 

With the aid of the radiations emitted by known quantities of radium and 
thorium set up at suitably chosen points, from 2 to 10 meters away, and all 
around our 7.64 cm lead screen we took readings when the screen was in place 
and when it was removed from our electroscope, and thus found that about 
2.4°% of the local radioactive rays from surface rocks and soils get through 
the lead screen and produce ionization within our electroscope. We then took 
a series of land observations in various localities, situated in widely different 
latitudes and at varying elevations from sea level up to 14,100 feet (the top 
of Pike’s Peak), half a dozen or more readings being in general taken at each 
locality over a period of several days, first, when the lead screen was in place, 
then when the screen was removed. By comparing these observations with 
those taken at the same levels beneath the surfaces of snow-fed lakes the 
water equivalent of the lead screen was quite accurately determined, as shown 
below, and in this way the depth-ionization curve was reliably extended up- 
ward the equivalent of about a meter of water above the highest point ob- 
tainable in Gem Lake. This last meter proves to be of great significance for 
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the purposes of the second objective ($1) and in general for the interpretation 
of the cosmic radiations. Table III contains the record of these land-readings, 
which were taken during the summer and fall of three consecutive years, 
1928, 1929, and 1930. The two readings on Pikes Peak taken two years apart 
are rather noteworthy. They differ by 27, but in view of the difference in 
the barometer reading this difference is not only in the right direction but of 
the right amount. 

The figures given in the last and the third from the last columns are in all 
cases the means of from three to nine different readings, the fluctuations in 
which are illustrated, for example, by the nine consecutive readings taken on 
Mt. Manitou at about two hour intervals beginning at 10:30 a.m. These nine 
readings, in ions per cc per sec, run 54.7, 52.8, 55.5, 54.3, 52.5, 53.7, 52.4, 52.2, 
54.3. Mean =53.6. 


5. Tut Deptu-lONIZATION CURVE AND ITs SIGNIFICANCE 


The graphical representation of the results in Tables I and III is given 
somewhat inadequately in Figs. 2 and 3,—inadequately because no small 
scale graph can reflect the consistency and precision of these readings. For 
this reason the readings themselves have been given in Tables I and III so 
that the reader may plot his own large-scale graph if he so desires. However, 
two important results stand out immediately and conspicuously from the 
graphs: 

1. At great depths, i.e., between 40 m and 80 m (see Fig. 3) the readings 
are so consistent as to show that even a blanket of 80 m or 262 ft. of water is 
insufficient to absorb completely the cosmic rays. The curve has here reached 
a value of 2 ions per cc per sec., not a fiftieth of its value at Pike’s Peak, but 
it is still falling. In order to find the zero of the instrument, or the ionization 
when all external rays have become absorbed and at the same time the absorp- 
tion coefficient of these hardest rays, we analyzed the curve by the trial and 
error method between 40 m and 80 m with the aid of the Gold tables, and 
found that with a zero of 1.2 J perce ‘sec.and an absorption coefficient of 0.028 
per meter of water, the whole long stretch of curve between 40 m and 80 m was very 
accurately reproduced. In a word, our curve shows at its lower end just such a 
band of hard rays as we had been looking for, and a single coefficient is adequate 
for the whole range between 40 m and 80 m. 

The significance of the value of this coefficient will be discussed presently, 
but the very existence of such a coefficient means that a hundred meters 
farther down, i.e., at 180 m, an ionization of about 0.03 J per cc/sec. should 
be observable by an instrument capable of detecting such an amount. We 
checked this conclusion qualitatively in 1928 by taking our electroscope and 
its lead shield down intoa shaft 185 m deep, beside Lake Arrowhead, and find- 
ing, after allowing for the local rays, a fraction of an ion still left for the cosmic 
rays, though, on account of the location of the shaft beside the lake we could 
not reliably estimate the equivalent water depth. 

However, Regener'® has reported more dependable deep water observa- 


10 E, Regener, Die Naturwissenschaften 15, 183 (1929). 
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tions which are in substantial agreement with the results here given. He re- 
ports his results in volts per hour, which is nearly the same as our ions per cc 
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per sec., since the multiplying factor of 1.171 reduces in our case J volt/hour to 
I per cc/sec. His lowest really detectable reading is at 186 m, though he takes 
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one observation at 230 m. He obtains both his zero and his absorption coeffi- 
cient by essentially the same trial and error method that we use, and we take 
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it that both he and we may have an uncertainty of as much as 0.2 J per cc/sec. 
in our zeros. Such a change in our zeros, if taken in the right direction for 
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both of us, would bring our coefficients into fairly good agreement. His value 
is 0.018 per meter of water. However, for reasons which we shall presently de- 
tail, the precise value of the coefficient for this most penetrating cosmic radia- 
tion is not particularly significant. The important result that appears both 
from Regener’s deep water work and from our own is merely that there exists 
a component of the cosmic rays which can penetrate as much as 180 meters 
of water and which is therefore approximately twice as hard as the hardest 
component we had before directly observed on our former curve namely, 
w=0.05 per m of water. 

2. The second immediately striking result that appears from our curves 
is that the quite accurate observations we have now taken at high altitudes— 
up to 14,100 feet, show that the absorption coefficient now rises far above 
anything the Dirac formula will inany way permit. The absorption coefficient 
shown at the top of our 1928 curve was 0.22 per m of water, that at the top 
of our present curve taken in water is 0.27 per m of water, and that at the 
top of the land-curve is 0.35 per m of water. This new curve then brings out 
more strikingly than ever the banded structure of the cosmic rays. 

We have analyzed this new curve with the aid of the Gold tables precisely 
as we did the former one in our 1928 article. We at that time found that the 
peculiar sharp bend in the curve at between 10 and 12 meters could not be 
reproduced without invoking three definite bands having roughly the relative 
frequencies 1, 4, 8. The new curve revealed the same necessity through the 
shape of its upper and middle portion, while its shape at great depths requires 
the introduction of a fourth band, as indicated above. Fearing that we might 
have become special pleaders for our banded structure, we asked Dr. Bowen, 
who had not thus far helped in this kind of analysis, to start from first prin- 
ciples with our curve and see, without any suggestions from us, what kind of 
structure it demanded. He proceeded without reference to any theory to 
build up with the aid of the Gold tables our observed curve out of four com- 
ponents—no smaller number would do—and in such a way that the syn- 
thetic and the observed curve fitted exceedingly nicely from one end to the 
other. 

The components from which the synthetic curve was thus built up to 
vield the ionization observed in our electroscope were as follows: 


TABLE IV. Assumed absorption coefficients and intensities of synthetic curve. 





Assumed Jo at Total J=Jo/u 





coefficients top of atmosphere 
0.03 33 1100 
.10 80 800 
.20 130 650 


.80 141 ,000 





The foregoing of course assumes that the ionization has its maximum 
value Jo at the surface of the atmosphere (see below, where this is shown to be 
incorrect), and the last column then shows that the total energy of formation 
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of the iron group, the silicon group and the oxygen group is of the same gen- 
eral order of magnitude. This checks reasonably well with that we know of 
the relative abundance of these elements, and it would probably check better 
if suitable corrections could be made for the foregoing incorrect assumption. 
Table V then gives a comparison of the calculated and observed ionization. 


TABLE V. Comparison of synthetic and observed curves. 


Depth Calculated Observed Difference 


in meters 

7.5 89.7 90.8 —1.0 
8.0 70.6 70.6 0 

9. 48.2 48.1 — .l 
10. 37.1 36.9 — .2 
12 26.4 26.5 + .1 
15. 19.1 19.1 0 
20. 12.53 12.55 + .Q2 
25. 8.87 8.75 — .12 
30. 6.58 6.56 — .02 
40. 3.93 3.83 — .10 
50. 2.49 2.62 + .13 
60. 1.63 1.88 + .25 
70. 1.11 1.29 + .18 
80. e 80 + 0.6 


|| 
|| 


Not only can the observed curve not be fitted accurately with less com- 
ponents than four, but also the four must be roughly of the foregoing type, 
though of course the lower bands can be split up into a finer structure if de- 
sired, i.e., 0.03 can easily be replaced by 0.02 and 0.04, for example. Not very 
much liberty, however, can be taken with the upper bands. 

We shall now compare the results of this purely empirical study of our 
curve with the results computed by the Nlein-Nishina formula," the Einstein 
equation mc? = E, and Aston’s curve. This formula has the form 


2rNery 1+ {= + a) 
== gs 


1 
2 log (1 + 2a) | 


m*c4 \ a? 1+ 2a a 
1 1+ 3a \ 
+ — log (1 + 2a) — : 
2a (1 + 2a)?S 


Where a=hv/mc?, N=no. of electrons per cc. The values of a for the atom- 
building processes, according to Aston, are 


IT-TTe = 0.029 + 5.479 X 10-4 = 52.9 
ITO =0.1245 + 5.479 X 10-4 =227 
IT—Si=0.232 +5.479 x 10-4 =423 
II— Fe =0.48+5.479 X10°*=876 


5.479 X 10-4 is the atomic weight of an electron obtained from e/m spectro- 
scropically determined. The numerical value of the constant factor is then 


" Nature 122, 399 (1928). 
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2rNe!t ae a) 
mc4 c4 m 


2X x X 6.004 X 10% & 10 





: ee —— & (4.77 XK 10° !°)? & (5.279 K 10')? 
(2.998)* K 10% XK 18 


= ().16614. 


Then Table VI is the comparison of uw obtained from the Klein-Nishina form- 
ula and those of the synthetically obtained curves. The closeness of the agree- 
ment at the top is most significant. The progressive departures as the atoms 
become heavier look at first sight like a dithculty, but the next section goes 
into possible causes of this behavior. 


TABLE VI. yw tn meters of water. 


Observed 


Computed 
H—He 0.7957 0.80 
H—-O . 2409 .20 
H-—Si .1418 10 


Hi—Fe .0754 .028 


6. CONDITIONS OF EQUILIBRIUM BETWEEN A BEAM OF PHOTONS AND ITS 
SECONDARIES PRODUCED BY THE COMPTON PROCESS 


When a beam of photons strikes matter it is obvious that a certain thick- 
ness of matter must be traversed before the beam gets into equilibrium with 
its secondaries, tertiaries, etc., this condition of equilibrium being attained 
when as many of each kind of secondary is disappearing per second from the 
beam as is forming per second in it. While this process of getting into equi- 
librium is going on, the absorption coefficient of the beam, as measured by the 
rate of change with distance of the ionization produced per cc, is obviously 
smaller than it can be after equilibrium has been reached. Further, the ab- 
sorption coefficient of the pure photon beam when it first strikes matter 
must be the same as that of the beam after it has got into equilibrium with its 
train of secondaries, provided the secondaries are more absorbable than the 
primaries,” for the reason that in this equilibrium condition the percentage of 
these secondaries is not changing at all as the beam moves on, the only ele- 
ment that is so changing being simply the number of primaries; but this is pre- 
cisely the situation in which the beam found itself when it first entered mat- 
ter. 

2 The apparent absorption coefficient when equilibrium is reached is equal to the absorp- 
tion coefficient of the primary or the secondary, depending upon which has the smaller coeffi- 
cient, the general relation being 

I= = ee — eH), 
we — 


Of the last two terms that having the larger coefficient will die out, leaving the effective coeffi- 
cient the one having the smaller value. It is practically certain from Bothe and Kolhiérster’s 
work that the coefficient of the beta-rays at sufficiently high frequencies approaches that of the 
photons, and it is entirely possible that it may even fall helow it for the hardest rays. 
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One of us has recently shown that the cosmic rays enter the earth’s at- 
mosphere as streams of pure photons. This means that the ionization in a 
closed vessel should not be a maximum at the top of the atmosphere, but 
that there must be an optimum position somewhere beneath the top where 
this maximum ionization is reached This is precisely what the 1922 high 
balloon flights of Millikan and Bowen, when taken in conjunction with Hess 
and Kolhérster’s earlier and lower ballon flights proved experimentally to be 
the case. For in the 1922 flights a recording electroscope rose to a hight of 
15.5 kilometers, at which height 0.92 of the earth’s atmosphere had been left 
behind, and the total ionization recorded by the self-registering mechanism 
proved to be only about one-fourth of that calculated from the absorption 
coefficient of 0.57 per meter of water, which had been found by Hess and Kol- 
hérster in rising in manned balloons to from 5 to 9 kilometers. The discrep- 
ancy seemed to be eliminated in 1923 when Kolhérster in experiments in the 
Alps got a coefficient only about 0.25 per meter of water, a value not in con- 
flict with Millikan and Bowen’s high flight. But up to the present day, 
though the lower value has been accepted, it has remained a mystery why the 
earlier European flights yielded such high values. This is now quite clearly 
explained, for the completely reliable curve shown in Fig. 1 is at the top rising 
quite as fast as did the Hess-Kolhérster curves, while the 1922 work shows un- 
ambiguously that it cannot continue to do so up to 15.5 kilometers. In other 
words, the coefficient has passed through a maximum somewhere between 
these two levels, and at 15.5 kms. has fallen back again to low values. This 
is obviously what, from the foregoing considerations, it must do if the cosmic 
rays enter the atmosphere as pure ether waves. 

It will be seen from the foregoing that the Klein-Nishina formula, com- 
bined with Aston’s measurements and Einstein's equation, yields quite accu- 
rately the observed absorption coefficient of the most absorbable band, when 
it is assumed that that band arises from the synthesis of helium out of hydro- 
gen. Also there is good reason to assume that at the level corresponding to 
the top of our curve this radiation has already reached a condition of equilib- 
rium with its secondaries, so that the comparison is probably here legiti- 
mate. The less absorbable the radiation, however, the farther must it tra- 
verse matter thus to get into equilibrium with its secondaries, and it is most 
illuminating to see how the absorption coefficients computed by the Klein- 
Nishina formula (see Table V1) for the formation of oxygen, silicon, and iron 
out of hydrogen are progressively higher than the observed coefficients as ob- 
tained from the curve, thus indicating that these progressively harder rays 
are farther and farther removed from the situation in which they have tra- 
versed enough matter to get completely into equilibrium with their seconda- 
ries. Further, such attainment of equilibrium should become more and more 
difficult the nearer the absorption coefficient of the beta-rays released by 
Compton encounters with photons approaches that of the photons them- 
selves, and Bothe and Kolhdérster’s recent experiments" show that this condi- 


13 R, A. Millikan, Phys. Rev. 36, 1595 (1930). 
4 Bothe and Kolhirster, Zeits. f. Physik 56, 751 (1929), 
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tion is somewhat nearly approached for the harder rays, though it is probably 
not so for the softer components. It may be that for these very penetrating 
radiations the secondary electrons are more penetrating than the primary 
photons and consequently even when equilibrium is reached a lower absorp- 
tion coefficient should be observed than that called for by the Klein-Nishina 
formula.” 


7. EVIDENCE FOR ATOM BUILDING 


In a word, then, the general qualitative evidence that the cosmic rays are 
due to the formation out of hydrogen of the only four abundant groups of ele- 
ments that there are, namely, helium, oxygen (C, N, ©), silicon (Mg, Al, Si) 
and iron (the iron group), which elements, barring hydrogen, constitute more 
than 99 percent of all matter;' is quite extraordinarily good, but this evidence 
only becomes quantitative in the case of helium. That this band, however, 
which contains within itself probably more than 90 percent of all the cosmic- 
ray energy, has so closely the absorption coetticient predicted for it by the 
Klein-Nishina formula, taken in conjunction with Einstein’s equation and 
Aston’s curve, is exceedingly significant. By a process of exclusion we are 
well-nigh forced to adopt the synthesis of helium out of hydrogen as the 
origin of this cosmic-ray band; for the Einstein equation and Aston’s curve 
leave us no other alternative, provided the NKlein-Nishina formula yields a 
result of even the right order of magnitude for the relation between absorp- 
tion coethcient and photon frequency or energy. The only other act that has 
been suggested, namely, the falling together of a positive and negative elec- 
tron, actually releases an energy 35 times that observed as obtained through the 
Klein-Nishina formula. This formula has been proved by Millikan and 
Bowen,'® by Chao,'? by Tarrant!’ and by Meitner'® to be approximately 
correct for the gamma rays from Th C’’ while the cosmic ray band under 
consideration, in accordance with our direct measurement, 7s but five times 
as penetrating as these gamma rays, so that wholly apart from its theoretical 
credentials, the extrapolation from Th C”’ up to the least penetrating cosmic- 
ray band is not a very long one. In other words, the Klein-Nishina formula 
ought to hold reasonably well for this softest cosmic-ray band. 


8. PARTICIPATION OF THE NUCLEUS IN Cosmic-RAY ABSORPTION 


The Wlein-Nishina formula, however, cannot be rigorously correct, for it 
makes the absorption proportional to the number of extra-nuclear electrons. 
We reported at the fall meeting of the National Academy in 1928 our definite 
evidence that the nucleus plays a role in cosmic-ray absorption. This evi- 
dence is found in the two curves of Fig. 2. If the mass absorption law held, the 
water equivalent of our lead screen would be 0.85 cm, this being obtained 


16H. N. Russell, Astro. Phys. Jr. 70, 11 (1929). The enormous abundance of H and He 
is the most striking feature of this article. 

16 Millikan and Bowen, Proc. Nat. Acad. 16, 421 (1930). 

17 Chao, Proc. Nat. Acad. Sci. 16, 426 (1930) and Phys. Rev. 36, 1519 (1930). 

18 Tarrant, Proc. Roy. Soc. 129, 342 (1930), 

1? Meitner, Naturwissenschaften 18, 534 (1930). 
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merely from the thickness of the screen (7.64 cm) and the relative densities 
of water and lead. The actual water equivalent of the lead screen for the 
softer cosmic rays is seen from the distance apart, on the x-axis, of the upper 
parts of the two curves, to be 122 cm. This distance would of course be ex- 
pected to increase for harder rays, and the curves show that this is indeed the 
case, since for rays of the hardness that are found at sea level (10 m) the water 
equivalent of the lead has become 170 cm. Chao" has recently brought to 
light a similar behavior for gamma rays, though of less magnitude, thus sug- 
gesting again the identity in nature of the gamma rays and the cosmic rays. 
The existence of such nuclear influences on absorption means that the Klein- 
Nishina formula must itself be but an approximation. 


9. Tue UNIFORMITY OF DISTRIBUTION OF THE Cosmic Rays 


One of us has recently presented at length new evidence that the small 
fluctuations that have been reported in the intensity of the cosmic rays are 
due merely to changes in the thickness of the absorbing atmospheric blanket 
which surrounds the earth, and that the cosmic rays themselves are streaming 
into the earth with entire constancy and uniformity of distribution over the 
celestial sphere. The correctness of this conclusion with respect to latitude 
could scarcely be more beautifully attested than by the left-hand curve of 
lig. 2. The observations there presented were taken with the same instru- 
ment under identical conditions as to observational technique, but at times 
extending over three different summers and in widely different altitudes and 
latitudes, the latter ranging from 34 to 59 degrees north. Yet they fit with quite 
surprising exactness one single tonisation-depth curve. The observations at 
Lake Louise in Canada, latitude 511%, at Colorado Springs, latitude 391g, and 
at Lake Arrowhead, latitude 34, are especially comparable because taken at 
about the same heights, as are also the near sea-level observations at Pasa- 
dena (latitude 34) and Churchill, Manitoba (latitude 59). 

This entire constancy in distribution of the cosmic rays is their most sig- 
nificant as well as most amazing property, and must mean, when taken in 
connection with their absorption coefficient, first, that the temperature exist- 
ing even in the atmosphere of the sun, whence alone they could get to us 
from this star, and the same is true for other stars, is inimical to the union of 
hydrogen into the heavier elements, for hydrogen is present in enormous 
quantity in the sun’s atmosphere. In the second place, these facts must mean 
that the cosmic rays do not originate in any places in the universe from which 
they are obliged to come to us through any appreciable amount of matter what- 
ever. If they had done so they would on entering the earth be partly beta-rays 
and partly photons, and, on account of the earth’s magnetic field, the beta- 
ray part would of necessity be stronger near the magnetic pole than at lower 
latitudes. But no trace of such an influence can be discovered! That, however, 
all about us “in the depth of space” hydrogen is somehow uniting into helium 
at least seems to us to be convincingly shown by the foregoing data, and that 
it is also uniting into the only other measurably abundant elements, the oxy- 
gen group, the silicon group, and the iron group, is strongly indicated, though 








252 R. A. MILLIKAN AND G. H. CAMERON 


precise quantitative proof becomes here impossible because of the change— 
always a decrease—in the absorption coefficient of a beam of highly penetrat- 
ing photons while it is getting into equilibrium with its Compton secondaries. 
This phenomenon—the existence of which was demonstrated by the 1922 
Millikan and Bowen experiments—renders futile a more precise analysis of 
our curve than we have given above. 


10. SUMMARY 
The foregoing results may be summarized as follows: 


1. Much the most dependable ionization-depth curve which we have thus 
far presented, especially at very high and very low elevations, has been ob- 
tained and analyzed. 

2. This curve shows, as did its predecessor, an unmistakable banded 
structure, and it further presents excellent quantitative evidence that the 
most intense and least penetrating band is due to the particular photon re- 
leased when four hydrogen atoms unite to form an atom of helium. 

3. Excellent evidence has been presented that these cosmic-ray bands 
enter the earth’s atmosphere as ether waves and must penetrate far into it be- 
fore getting into equilibrium with their secondaries; also, that until such equi- 
librium is attained the observed absorption coefficient is smaller than that of 
the initial monochromatic radiation. This means, as shown in the 1922 Mil- 
likan and Bowen experiments that there is a level beneath the top of the at- 
mosphere at which the ionization due to a pure cosmic-ray beam is a maxi- 
mum, this maximum moving rapidly farther down as the frequency of the 
initial photon increases. 

4. The observed cosmic-ray curve has been shown to be consistent with 
the theory that it is made up of four bands due to the formation out of hy- 
drogen of helium, oxygen, silicon, and iron the only atoms of sufficient abun- 
dance to render the radiation released by their formation detectable anyway, 
and it has been shown that the differences between the calculated and ob- 
served absorption coefficients of the photons produced by the formation out 
of hydrogen of oxygen, silicon, and iron are very nicely explained by the non- 
equilibrium theory given in (3), the departures all being in the right direction 
and increasing with frequency in the right way. 

5. The constancy and the uniformity of distribution of the cosmic rays 
over the celestial sphere has been again brought strikingly to light, and the 
significance of this for the place and mode of origin of the cosmic rays has been 
again pointed out. 

6. Some participation of the nucleus in the absorption of cosmic rays has 
been experimentally established. 

We wish to express our appreciation to Professor Bowen for assisting, as 
indicated above, in the analysis of the curve. 
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ABSTRACT 


The radiation of multipoles formed by putting together elementary dipoles is in- 
vestigated according to the classical theory and wave mechanics. It turns out that the 
radiations are of two kinds: 

The first part has the same frequency as the dipoles and is present only if 
the multipole structure is present even without vibration. 

The second part is always present. If 7, the order of the multipole, is even then 
the frequencies are 2v, 4v, up to mv. If 1 is odd then the frequencies are v, 3», up to nv. 

The relative intensities of subsequent multipoles are proportional to the square of 
the ratio: dipole amplitude to wave-length, as has been found in other cases. This cal- 
culation would apply to the infrared radiation of molecules. 


ECENTLY there have appeared a number of papers on the radiation of 

higher poles;! this radiation seems to explain =he presence of forbidden 
lines.2 Although e.g. Rubinowicz’ paper gives rather complete general for- 
mulas for the classical theory, it seems not useless to bring out the physical 
facts as clearly as possible by considering the case of a combination of linear 
oscillators, as they are present in molecules like CO: (two) or CCl, (four). 


I. CLASSICAL THEORY 


We consider quadrupoles, but there is no difficulty in generalizing to 
higher poles. The simplest quadrupole consists of two equal dipoles vibrating 
in the direction of the line joining them but with opposite phases. 

If the direction of vibration is the z-axis, the field of one dipole can best 
be described by a Hertz-vector Z 


Z,=Z,=0 
cos 2av(t — r/c) 
Z,=Z=€4 (1) 
r 
1 az 1 az . 
E= AZ — , H =— rot . (2) 
c at C ot 


If we have now the two dipoles the distance / apart along the Z axis and 
vibrating with the same amplitude but opposite phase, we will get the re- 


11, Placinteanu, Zeits. f. Physik 39, 276 (1926). A. Rubinowicz, Phys. Zeits. 29, 817 
(1928); Zeits. f. Physik 53, 267 (1929); 61, 338 (1930); I. Blaton, Zeits.f. Physik 61, 263 (1930); 
I. Bartlett, Phys. Rev. 34, 125 (1929); A. F. Stevenson, Proc. Roy. Soc. A128, 591 (1930); L. 
Huff and R. W. Houston, Phys. Rev. 36, 842 (1930). 

2 1. Bowen, Proc. Nat. Ac. 14, 30 (1928); Astrophys. J. 67, 1 (1928); R. Frerichs and I §. 
Campbell, Phys. Rev, 36, 1460 (1930). 
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sultant Hertz vector at any point by taking the difference between the con- 
tributions from the two dipoles. The resultant Hertz vector will consist of 
three parts: 

acos 2rv(t — r/c) 1 


Z, = — ——_——_——_—_ — cos. (3) 
r r 


This part will not contribute to the radiation in distant places on account of 
the factor //r 
asin 2rv(t— r/c) l 


2, = 2e———_ ——— —= 08 §, (4) 
r r 


This part does give radiation which will represent a spherical wave of the 
same frequency as the one which would be emitted by a dipole, although its 
intensity varies in a different manner with @ and its amplitude is diminished 
by the factor 27/]/X. The formula is 


; 2rl 4r° sin 2rv(t — r/c) 
+ iH, = Ke = — cos 04 —a 
N iee 





sin ot ; (4a) 
r f 

The total intensity emitted is the fraction (87?/5) (2/X) of that emitted 
by one dipole. 

The third part is the most interesting one; it is always present, even if the 
two dipoles are not apart (/=0) when at rest, that is to say if there were 
e.g. two negative charges vibrating on opposite sides of a positive charge 
with the same position of equilibrium. It arises from the difference of posi- 
tion caused by the vibration itself. That is to say, we have to take 7 now not 
as the distance from the point of observation to the dipole, or quadrupole, but 
to the vibrating charge. If we call the deviation from the equilibrium posi- 
tion in a given moment 2’, we have to put in (1) 


r=ro+3' cosé (5) 
and have for the dipole 


cos 2rv(t — ro/c) + 22/As’ cos @ sin 2xv(t — r/c) 

















L=@ i al 
r 
a cos 2rv(t — ro/c) 2ra* sin 2rv(t — r/c) cos 2xv(i — r/c) 
= + cos 6 —— (6) 
r r 
a cos 2rv(t — ro/c) 2ra* cos @ sin 242v(t — ro/c) 
r r 2 r 


If the two dipoles of opposite phase are present, the first member drops out, 
and we have 


2r cos@ sin 22v(t — ro/c) 
Z3 = —<¢ — a 





(7) 


r r r 
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RADIATION OF MULTIPOLES 2 


That represents a spherical wave, with the same angular distribution as 
Z2, but with the double frequency and an amplitude proportional to the square 
of the amplitude of oscillation. The total amount radiated is (44?/A,?(a?)2/5) 
times the amount radiated by one dipole. 

Generalization to higher poles will give the following result: Assume that 
we have a pole whose order is . 1 is defined in this manner: Let the unit 
vector p; indicate the direction of one of the constituent unit dipoles which 
make up the higher pole. Let r be the unit vector in the direction of observa- 
tion. 

Then a pole of order 7 is one for which 


> (pir)"’ = 0 


for all rand n’<n and > (pir) "() for some r. In fact the left hand side is a 
spherical harmonic of order ». With this definition, a dipole has n=1, a 
quadrupole »=2, a tetrahedron n=3. 

Let us now define a vector which will depend only on the orientation of 
rtoasystem of coordinates fixed in one pole 


L, = > p.(pir) dls 


we will need later besides a scalar, defined by 
o, = (L,r’) 


where r’ is a unit vector perpendicular to r. 

We find then two different kinds of radiation. 
(1) If there is a “statical” multipole structure of order » and strength p, 
present, we get a radiation (Hertz vector) of the same frequency as the fre- 
quency of the constituent dipoles, but with the factor 


1 /2r \"" 
bee See © 
—( r ) P 


(2) If the multipole has its origin only in the motion (that is if the field in 
equilibrium corresponds to a pole of higher order), part (1) is absent. But 
there is always present a radiation with a Hertz vector 


1 a” gn} 


Z = — —L, cos""2av(t — r/c)—— cos 2av(t — r/c). 
nior or" 


If m is odd, we have 


1 a” ar n—1l 
Ee af 90 <1,(= cos” 2rv(t — (r/c)) 
n! r Ao 
(— 1)" 1)/2 a” ar n—1 2s=—n—1 N . 
on aiceneeteente —(= L,2 be ( ) co 2r(n — 2s)v(t — (r/c)) 6 (8) 
2"-!n! r \Xo \ 2-0 \S f 
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If mis even 


(—1)"/* a"/2ar\""! { tenn—8 n n 
as Qn Int or (=) ed 2X (( ") 7 ( — )) 
sin 2r(n — 2s)v(t — (r/o) (8a) 


Accordingly we get radiation which has the frequency (7 — 2s). 

The only similar case where this seems to have been stated explicitly® is 
the case of m equal particles moving equidistant on a circle, which arrange- 
ment of course forms a pole of » fold multiplicity. 


Il. WAVE MrECHANICS 


In this case it is simpler to treat absorption instead of emission. The 
transition to emission is always easily made by using as external field the 
Landé-Dirac ghost field. 

We assume a harmonic oscillator, whose Schrédinger equation is 

dy 8r°M drill dy 
in —2r*Mv*s*y = — . ee (9) 
dz? h h at 





with the solution 


y= Dicwpyen2tik vt ve = E20, (0b, 


with the abbreviation 
Mv\!? 
¢ = 2r{— S. 
h 


H, is a Hermitian polynomial. 


1 /4r Mv\!/\12 
b,, = (2-*{= ~) ), 
ni\ h 


If there is now a perturbation function present of the form 


V 2risvteyh 


= ad 
we find 
dors 8 Qri 


~ Tt ——— oe VW iad. 10 
i h if ViWirgsd6 ( ) 





If we have a light wave polarized parallel to s and propagated parallel to x, 
its electric field will be 


E, = E cos 2avs(t — x/c). (11) 


The perturbation function is then given in the first case (of a statical mul- 
tipole structure) by 
A on-l 
V = —¢.9.--—-E. 
n! axe 


3 J. J. Thomson, Phil. Mag. (VI) 6, 673 (1903); G. A. Schott, Electromagnetic Radiation, 
Cambridge 1912, p. 102. 

4 See f. e. E. U. Condon and P. M. Morse, Quantum Mechanics, New York, 1930, p. 47. 
A. Sommerfeld, Wellenmechanischer Ergiinzungsband, Braunschweig 1929 p. 17. 
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from which we get® 


dénss ri (=) amt E 1 f Viburad 
— = ——¢,/ — nOnE— |} whiWiysdz. 
dl h ’ Xo P n! — 


The latter integral gives of course the same selection rule as for a dipole s = +1 
and only modifies the intensity of the absorption by the same factors as in 
the corresponding classical example. 

More interesting is the case arising from the existence of the oscillatory 
We will put the perturbation function multipole. 


1 o”-} 
V = —3"c,e ——E, 
n! Ox"! 


which leads to 


dCkss 2nri so ecE f2xr\""! 1 n+l 
ee EE eg geo ome o, ne bidiisl ke” 
dt 





h nt! \Xo (=)" 
2x| — 
h “ws 
1C; Onf dv \Or-dDP 1 ) 
=— qt =) — T..” 
(mMvh)!? ni\Me 2*tsl2( kk + s)!)12 
where 
I 2 

Iy.s" = bil — { oY Wigsde - [ire Wy H yisde. (13) 

rOh+s & 


To evaluate this integral, we use first the following formula,® where an 
upper index j indicates j fold differentiation in respect to the argument 


HW, = 2kH,_, 
and therefore 
Hy, = 2k-(2k — 2) +--+ (2k — 27 + 2)My;. (14) 
If we write 


HH, = Ya ci 


i 


we have 


W,.?(0) = flash = R(k — 1)--- (k — jf + 1)Ay_,(0). (15) 
But on account of the formula’ 
i, = (- tyr ett 
dg* 


® Of course, we have been quantizing not the motion of one charge, but of a characteristic 
vibration (normal coordination). As these are independent (no natural interaction), there are 
only transitions between different states of the same vibration possible, as G. Dieke has kindly 
pointed out. 
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we have 
d*-3 ~ ge 
HM, (0) = (— 1)4 i) ( 1-1 -) )! 
Ldc* ? 0 l! r=0 
= O tor k — j odd 
an (RD! | 
= (— 1) : for k — / even 
k—j\, 
( -F 
and finally 
k! | (g— a)! . 
a;*) = ——24(— I)(f-0F for k — i even (16) 
(sz)* k—i 
! 
C; 
= 0 for k — 7 odd. 
We now use Sommerfeld’s method? 
© d*+s 
[= os ——- 2) U7 ,.dé 
stg dé i dl 
a 
= e s* (CLT) dg 
iat d¢' 
x d +s h 
= } es ( Dae “) as 
“—x dé le i=0,1 
=OQOfork+s>k+unors>un 
ee) 
= >. f a®(itaitu-—1)---(itn—k — s\gite-’ ede forss n 
i —20 
if we now put 
itnx-—-k—-s=z! 
we have 
n—s (k) *xX ss 
[= ) yee AL+k+s)-- wf e felds, 
l=0,1 a 
rt . . . . (2/)! oa . . 
rhe integral is 0 if 2is odd and z' -— —if lis even = 2). 
92) . 
Accordingly, we have 
Ty, =Oifs —an odd 
ome (4))! 
=m? SY ———_(2j+k4+s)---2j 
j=0 241(27)! 
(n — s — 23)! 1 
(— 1) 2+ Gs—») i Aa hd) i i 
t—s ; ~rkts—anyP 
est - —j)! 


® Condon and Morse, reference 4, p. 50. 
7 Sommerfeld, reference 4, p. 59. 
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»)! 


4 (4)! 1 
= attt(—1)t-eoegeioms at D> (—1)! = (ths)! 
a 


7=0 (2) og = 
. (17) 
(2g — 27)! 1 
(q—Jj)! (G{tk-—gq)! 

with m—s = 2g. 
That means that for » odd the transitions 1, 3,5 - - -2 and ” even the transi- 
tions 2,4, 6, - - -m are permitted, and accordingly the frequencies v, 3y, 5y - - - 
nv or 2v, 4v - - -nv absorbed and emitted, in analogy to the classical theory. 


In normal cases, the formula are not as terrible as they look. For example, 
for n =2 (quadrupole) we have s = 2 
: (hk v: 2)! , 
Tyg = ()'/72*R!— = (r)'*2*(k+1)(k +2 
(k!)? 
If n =3 (tetraeder) there are possible two transitions 
s=1, frequency v 


j (hk + 1)i2! 4! 1 (k + 3)!) 


Ty 2) = — (x)*/228+2R! —- —— — 
™ Va-pe (202 4 e+ De$ 
. j 1 (k + 2)(k + 3)\ 
= — (x)'/22*+2. 2h2(k + 1) — — 
| 4 k+1 f 
s=3 frequency 3v 
kik + 3 
Tx .3 = (x)! yeh = (w)'/*(k + 1)(k + 2)(Rk + 3) 
In general, we have for s =n, frequency nv 
kiik+s)! 
Tien = (w)'/?224— i - = (w)'/*2*(k + 1)--- (Rk + 2) (18) 
which will make the factor in (12) proportional to |k-(k+1) - - -(k+n)]! or 


k"'?, if kR>>n; this corresponds to the fact that the amplitude of light emitted 
classically is proportional to a" as a?~k ~the energy of one dipole. 

The relative transition probability from k+s to k for a pole of the ” order to 
the same probability of a pole of (7 —1) the order is 


1 on” (—*:) ( ~)- 
N On?-1 ra 1) Me 
hv/ Mc? is the ratio of the energy of vibration to the intrinsic energy and can 


be rewritten 
M/2om: 7'max 1 - ( =) 1 (—") ( a y 
oe al on he | Sh, 
Mec? 2 c 2 Cc No 


Rubinowicz? has found a similar result for atoms. 
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THE MOTION OF ELECTRONS IN ARGON 
By Hl. B. Wantin 
DEPARTMENT OF Puysics, UNIVERSITY OF WISCONSIN 
(Received December 27, 1930) 


ABSTRACT 
The motion of electrons in argon containing 1.53 percent Ns has been investigated. 
For low fields the mobility is a constant, but rises rapidly as the field is increased due 
to the increase of the electron free path with increasing thermal velocity. 
The mean free path of electrons in pure argon and in thermal equilibrium with the 
argon is 0.385 cm at 1 mm pressure. 


N A series of earlier papers! results of investigations on the motion of elec- 

trons in No, He, He and CO have been presented. The method used was 
the alternating potential method of measuring mobilities. Since in the case 
of the electron motion through a gas the mobility is in general not a constant, 
this method does not lead to a direct determination of the mobility. It is pos- 
sible, however, as was done in this earlier work, to make an indirect compari- 
son of the experimental results with the theoretically calculated mobilities. 
In all the cases presented so far, with the exception of N., a good agreement 
with the theory proposed by Kk. T. Compton? was shown to exist, provided it 
is assumed that the energy loss of an electron on impact with a molecule is 
greater than it would be if due to momentum transfer only. In N, the varia- 
tion of the electronic free path with the electron velocity enters into the ex- 
perimental results to such an extent that a check is impossible. Compton's 
theory takes no such variation into account. 

In the present paper results on argon corresponding as much as possible 
to those on the other gases, are presented. 

The argon used in these experiments contained, after being freed of oxy- 
gen, 1.53 percent nitrogen. This was not removed, but instead, as shall be 
seen, corrections were made for its presence. All the measurements were made 
at a pressure of 760 mm. 

When the mobility is a constant, it may be calculated exactly from the 
voltage intercept of the mobility curves by the expression 


K = rnd?/(2)'?V. (1) 
If the mobility is a function of V as it will be for electrons, this equation 
gives an average value determined by the relation 


1/2n | L/2n | 
f —K sin 2rnt dt = { f(z)— sin 2rnt dt. 
0 d 0 d 


Compton’s theory gives 


f(v) = 





a 





(1+ (1+ BV} (2) 


1H. B. Wahlin, Phys. Rev. 23, 169 (1924); 27, 558 (1926); 35, 1568 (1930). 
2K. T. Compton, Phys. Rev. 22, 333, (1923). 
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In Fig. 1 the average values of K calculated from Eq. (1) are plotted 
against the field strength (V/d) in which they were determined. For very 
low values of the field the mobility is a constant, and over this range the 
values calculated from Eq. (1) are nearly exact. In fields above 2 volts/cm 
the mobility rises sharply. The values found for this range serve merely to 
illustrate qualitatively the rapid increase in the mobility. 


260 
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Fig. 1. 


An examination of Eq. (2) shows that according to Compton’s theory one 
should expect a decrease in the mobility with increasing voltage due to the 
increased thermal velocity of the electrons. In order to account for the dis- 
crepancy between theory and experiment it is necessary to assume a very 
rapid increase in the electron free path in argon with increasing thermal veloc- 
ity. This is in agreement with the free path determinations of Townsend and 
Bailey® and others. 

The results shown in Fig. 1 serve in one way as a check on Compton’s 
theory. The equation giving the mobility as a function of the field strength 
(Eq. (2)) demands that for low field the mobility shall be a constant; i.e. 
that for low field the electrons shall be in thermal equilibrium with the gas. 
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The fact that the mobility is nearly a constant as V approaches zero shows 
that apparently the thermal velocity does not vary sufficiently to produce 
any appreciable change in the free path. 

In Fig. 2 two typical mobility curves are presented. The curve taken at a 
frequency of 14,100 cycles/sec is of the same character as those found for 
other gases. The one taken at the higher frequency is characteristic for argon. 
This frequency gives two values for the mobility. The variation in the shape 
of the mobility curve is due to the Ramsauer effect which in argon is much 
more pronounced than in other gases studied. 


3 J.S. Townsend and V. A. Bailey, Phil. Mag. 44, 1033 (1922), 
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The limiting electron mobility in argon containing 1.53 percent Ne» is for 
low fields 59,600 cm/sec/volt/cm. The mean free path calculated from this 
value assuming 

kK = 0. 75eX,; me 


is 0.363 cm at 1 mm pressure. 
To correct for the presence of N» the following expression may be used 
700 i’. ge 


+ 
Am X, he 


where P,, is partial pressure of No, P. the partial pressure of argon, A,, the 
mean free path at 1 mm in mixture, A, the mean free path at 1 mmin N. and 
X,, the mean free path at 1 mm in pure argon. This expression was used by 
Townsend and Bailey with a mixture of H, and A and was justified by them 
for this case. The corrected value for the mean free path when the electrons 
are in thermal equilibrium with the argon atoms is (assuming A, =0.0996 cm) 
0.385 cm at 1 mm pressure. 

In Fig. 3 the author's value of \ (indicated by VY) together with Townsend 
and Bailey's values are plotted as a function of the thermal velocity é. This 
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Fig. 3. 
curve shows the prominence of the Ramsauer effect in argon, and also serves 
to explain the shape of the mobility curves obtained at high frequencies. 
SUMMARY 


The results of the determinations of the mean free paths of electrons in 


gases, together with the values calculated from kinetic theory, are shown in 
Table I. 














TABLE I. 
Obs. KT Obs./KT 
H. 0.041 0.0741 0.553 
He .066 .117 .565 
CO .069 .0391 1.75 
Ne .0996 .0393 2.53 
9.2 


A .385 .0417 














As is seen there is no systematic agreement between the observed and the 
kinetic theory values. 
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DISPERSION AND REFRACTIVE INDEX OF NITROGEN 
MEASURED AS FUNCTIONS OF PRESSURE BY 
DISPLACEMENT INTERFEROMETRY 


By CLARENCE E, BENNETT 
DEPARTMENT OF PHysics, BROWN UNIVERSITY 


(Received December 15, 1930) 


ABSTRACT 

The Lorentz-Lorenz relation has never been satisfactorily checked by experiment 
for pressures above atmospheric, because of a lack of trustworthy data on the refrac- 
tive indices of pure gases as functions of pressure. Such data would also be particularly 
valuable in the study of polar molecules on account of Maxwell's law connecting the 
refractive index with the dielectric constant. Due, however, to the difference between 
optical and electrical frequencies the refractive index must be reduced to zero fre- 
quency before Maxwell's law can be applied. This necessitates a precise knowledge of 
dispersion which is not available for many pure gases. In the present paper a method 
of measuring simultaneously the dispersion and refractive index of a gas is discussed. 
The method is justified on theoretical and experimental grounds and applied in a 
practical way to the measurement of dispersion and refractive index of nitrogen over a 
range of pressures and at two different temperatures with a high degree of precision. 

The results show that over a seven atmosphere range the Lorentz-Lorenz relation 
is followed and that the dispersion is a linear function of pressure. The value of the 
Cauchy constant B is found to be 1.95 X 10-" per cm of mercury pressure. The index of 
refraction for nitrogen is measured for three wave-lengths at 0°Cand 20°C with a limit- 
ing error in (u—1) of one third of one percent. The dispersion information enables the 
values to be reduced to infinite wave-length with the result that a value for the dielec- 
tric constant is obtained which differs by less than two-thirds of one percent from that 
measured directly by Zahn. The values of « under NTP conditions for the wave- 
lengths 4811A, 5893A, and 6362A are 1.0002991, 1.0002969, and 1.0002967 respec- 
tively. 


INTRODUCTION 


HE problem of measuring the index of refraction of gases has been at- 

tacked by many observers in the past, but precise results have been ob- 
tained with only a few substances. The most accurate determinations have 
been made for air, an admittedly complex substance, and in practically all 
cases the observations have been made at atmospheric pressure. The data on 
air obtained by Meggers and Peters! with a Fabry-Perot interferometer over 
the whole visible spectrum range, undoubtedly represent the best work of 
this sort ever attempted. Trustworthy data on the index of refraction of pure 
gases as a function of pressure, however, were scarce at the time the present 
research was started. Gale? in 1902 worked on air, about the constitution of 
which he knew very little, up to a few atmospheres with a Jamin interferom- 
eter, and Magri* in 1905 went up to fourteen atmospheres with air also using 


1 Meggers and Peters, Bulletin of the Bureau of Standards, No. 14, 1918-19, p. 697. 
2H. G. Gale, Phys. Rev. 14, 1 (1902). 
3 L. Magri, Phys. Zeits. 6, 629 (1905). 
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a Jamin interferometer. Cheney,’ using pure gases, in 1927 measured the re- 
fractive indices of nitrogen, ammonia, carbon dioxide, sulphur dioxide and air 
over a temperature range of a few hundred degrees at one atmosphere by 
means of a Fabry-Perot etalon, but so far as can be ascertained, no one with 
the exception of Phillips? who worked with carbon dioxide, has published any 
reliable data on pure gases at pressures higher than one atmosphere. 

In the light of the Lorentz-Lorenz®* equation such information should be 
particularly significant. Furthermore as a result of the work of Debye’ and 
others in which dipole moments have been coordinated with dielectric con- 
stants, data on refractive index become especially valuable because of Max- 
well’s law connecting the dielectric constant with the refractive index. This 
necessitates a knowledge of the dispersion of the medium because dielectric 
constants are measured at very low frequencies compared with optical fre- 
quencies. 

Some years ago Carl Barus® made a few preliminary observations on the 
index of refraction of air with his displacement interferometer, as a result of 
which, it appeared possible to measure both refractive index and dispersion 
simultaneously. This work was only qualitative and no attempt was made 
to secure high precision, yet the experiments showed that such precision 
might be expected if the instrument were properly developed. It was the 
particular good fortune of the writer to be able to assemble an instrument of 
this type with the personal assistance of Professor Barus, and to adapt it for 
this sort of work. 


DISPLACEMENT INTERFEROMETRY 


The optical arrangement of the displacement interferometer is shown in 
Fig. 1. When the system is properly aligned so that the mirror J/’ is made 
accurately parallel to the virtual image of the mirror .\/ as seen through the 
half silvered mirror //S, the interfering light beams, if viewed through the tele- 
scope 7 on the front of which is mounted a ruled grating, display a beautiful 
interference pattern consisting of black concentric oval fringes superposed on 
a continuous spectrum. The optical path difference necessary to produce this 
result is a known function of the thickness, the index of refraction, and the 
dispersion of the glass half silvered mirror. 

These fringes are very interesting. They have been studied at length by 
Barus,’ and Birchby'® has more recently described and explained fringes of a 
similar nature. Like Michelson fringes they expand or contract with every 


4 E. W. Cheney, Phys. Rev. 29, 292 (1927). 

5 P. Phillips, Proc. Roy. Soc. London 97A, 225 (1920). 

° H. A. Lorentz, “Theory of Electrons.” 

7 Meggers and Peters, reference 1, credit Magri, reference 3, with preferring (u?—1)/ 
(u?+1)/(p)=C but this is a mistake resulting from a typographical error in Magri’s publica- 
tion as anyone can see by checking the calculations, 

8 P. Debye, “Polar Molecules,” 1929. 

® Carl Barus, Carnegie Foundation of Washington Publications, No. 149, I, II, III; No. 
229; No. 249, I, II; 1911 to 1917. 

10 W.N. Birchby, Proc. Nat. Acad. 10,, 452 (1924) also 13, 216 (1927). 
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motion of the micrometer screw, but continuous motion causes the whole 
system to move bodily along the spectrum from violet to red as the path from 
the micrometer mirror to the half silver is made shorter. Using for a half 
silvered mirror a plate of glass 0.25 inches thick, a motion of the micrometer 
mirror of 0.17 mm corresponds to a shift of the fringes from the red to the 
blue-green region of the spectrum. This means that movements of the fringe 
system along the spectrum corresponding to displacements of the mirror can 
be calibrated with a good micrometer. Hence the instrument is particularly 
well suited to measure displacements greater than can conveniently be meas- 
ured with a Michelson interferometer, because it is unnecessary to count 
fringes, vet the accuracy is about the same. 
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Fig. 1. Optical arrangement of a linear displacement interferometer. 


The fringes are visible simultaneously in all spectrum orders and remain 
so even for very wide slit openings, the visibility however, being generally 
best for reduced illumination. When sunlight is used as a source, the dark 
Fraunhofer lines constitute fixed lines of reference by means of which meas- 
urements of the displacement of the fringe system can be made. Ordinarily 
with a carbon arc the yellow sodium line is bright enough to stand out against 
the continuous spectrum background and serve as a reference line. 

By a geometrical consideration of the reflected and refracted rays Barus 
derived the fundamental relations of displacement interferometry. But let us 
now consider in quite an independent manner how the instrument can be 
used to measure the index of refraction of a gas. Suppose the interferometer 
to be so adjusted that the ovai fringes are centered on a given spectrum line, 
for example the sodium D line, when a column of gas of length e¢ is in one arm 
of the interferometer. Assume y; to be the index of refraction of this gas, wu. 
to be the index of refraction of the medium in which the micrometer mirror 
is displaced (air), and if there are any other substances in the path let their 
refractive indices be denoted by pe - - * Mn. 

With the gas in the path, the optical path difference from the half silver 
to the mirrors J/ and N is 


md = mie + F(us+ ++ wn) — P (1) 
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where J is the wave-length of the spectrum line on which the fringe system 
is centered, we is the optical path length of gas, F(uw.---y,) is the path 
length of all remaining substances in the path, and P is the optical length of 
the opposite path. F and P are absolutely independent of the nature of the 
gas in question and are not affected by evacuating the tube. 

When the gas is removed from the path (by evacuation of the containing 
tube) and the micrometer mirror is displaced to recenter the fringe system on 
the given spectrum line 


Wr = et wydn t+ Flus--+u,) — P (2) 


where wv for vacuum is of course one, and Ax is the air displacement of the 
micrometer mirror. 

The radial distance between the central spot and the first fringe ring is 
observed to be greater than that between any two successive rings. Hence in 
accordance with Cornu’s principle, the fraction of a fringe shift for the change 
of one wave-length must be a minimum at the center of the system, i.e. 

on’ an 

- =— = 0) 

Or Or 
Thus the fringes near the center of the configuration are relatively coarse and 
far apart, whereas the outermost ones are very fine and close together. 

Dividing Eqs. (1) and (2) by X, differentiating each with respect to X, set- 
ting the resulting equations equal to zero, subtracting one from the other, and 


solving, we get 
Om 
G iti — 1) — \— 
Or 


a OE ea (3) 


The denominator for all practical purposes is one, whence 


Ay Ou 
(u: — 1) = — 42.—- (4) 
e Or 


DISPERSION 


Eq. (4) shows that before the index of refraction of a gas can be deter- 
mined by displacement interferometry some knowledge of its dispersion is 
necessary. Further considerations, however, show that the dispersion can be 
obtained directly with the instrument from the values of the displacement Ax 
taken for two different wave-lengths in the spectrum. 

Consider the Cauchy dispersion relation 





h=A+ 
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The data of Meggers and Peters' show that for air, two arbitrary con- 
stants A and B are sufficient to represent the facts. The writer has taken 
these data for three different temperatures 0°, 15°, and 30°C and plotted 
(u—1) aginst 1/M* on a very large scale (a graph 50X70 cm) as shown in a 


<a 





Fig. 2. Dispersion data for air plotted from observations of Meggers and Peters. 


reduced photograph in Fig. 2. In each case the points lie rather well on 
straight lines, indicating that two arbitrary constants are sufficient to repre- 
sent the data. 

If now (u—1) =A +B/) whence 


Ou 2B - 
seca: TUNE lite. waa (5) 
Or 3 
Then Eq. (4) reduces to 
Ax 2B 
(ms - 1) = — + —. (6) 
€ A“ 


Now writing Eq. (5) twice for two different wave-lengths \; and A, and 
subtracting the one from the other 


) (= “) (7) 
(u: — po) = BL — —- — }. 7 
1 b AY de 


For these two wave-lengths Eq. (6) gives 


; Ax, 2B ; 

(a — 1) = — — — (8) 
e A,” 
; Axo 2B 

(uy — 1) = = - —. (9) 
e dA," 


Where B is assumed to be constant throughout the visible spectrum as the 
data of Meggers and Peters show to be true for air. Solving Eqs. (7), (8), (9) 
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(Aa, —_ Axe) 
ib ae es (10) 


36 - — 


Now consider again Eq. (6) rewritten so as to read 


x 2b 
-=(4h—1)+—-- (6a) 
e * 

Experiments on nitrogen which will be described later, show conclusively 
that Ax is a linear function of pressure when values of Ax are plotted against 
pressure. Therefore, since e is a constant, [(u—1)+2B/\*] must be a linear 
function of pressure, meaning that either each term separately varies directly 
with pressure or that there is some involved relation connecting the two terms 
in such a way that their sum varies directly with pressure if separately they 
do not. Since ?(u—1) is a first approximation to (u2—1)/(u2+2) we should 
expect (u—1) to be directly proportional to pressure by the Lorentz-Lorenz 
relation, and so B must, therefore, be a linear function of pressure. This 
means that B/p is a constant which can be determined directly by the differ- 
ence in the slopes of the lines Ax,/p=constant and Axs/p=constant. 

Dividing Eq. (10) by p it appears that 


(AN A “*) 
B ( Pp p 


he 8 (11) 


p . ( 1 -) 
A Ceeedineation 
\? » 


By here means B per cm pressure. The B used in the previous work can now 
be interpreted to mean B at pressure p. 

It follows then that the displacement interferometer can be used to de- 
termine the dispersion factor B of a gas, first finding By from the difference 
in the slopes of the (Ax, p) lines for two different wave-lengths and then multi- 
plying by the pressure to get the value of B corresponding to that pressure. 
After B has been determined (u—1) can be determined at the same pressure 
by Eq. (6). The value of (u—1) can be reduced to its equivalent at 0° and 
normal pressure if it is not so measured, by assuming the Lorentz-Lorenz 
equation to be at least approximately true, whence 

Me a = = -! (12) 
p Po 


(13) 


II 
" 
| 
-—— 
: 


(uo — 1) 


PRELIMINARY EXPERIMENT 


An interferometer with arms something over a meter long was assembled 
according to Fig. (1), and measurements of the refractive index of dry air 
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were made at atmospheric pressure using for the dispersion correction the 
value given by Fig. 2. The average of some twenty or more independent ob- 
servations gave a resulting value for (u)—1) of 0.0002921 with a mean devia- 
tion of +0.0000007. This is exactly the value given by Meggers and Peters.' 
Hence the conclusion is drawn that the above theoretical considerations have 
heen experimentally justified. 
DESCRIPTION OF THE APPARATUS 

The arrangement of the apparatus is sketched in Fig. 1. A photograph 
of the actual set-up in the basement of the laboratory is reproduced in Fig. 3. 
The interferometer mirrors were made from plates of optical glass 3 inches 
square obtained from Bausch and Lomb. The micrometer mirror J was fas- 
tened to the movable table of a small Michelson interferometer which had a 
well-made micrometer screw about 3 16’’ diameter of 0.5 mm pitch. The 





Fig. 3. Photograph of apparatus. 


Michelson inteferometer was also set up so that it could be used as an aux- 
iliary instrument as indicated in the figure, for the purpose of calibrating the 
screw over the range actually used, by counting Michelson fringes from an 
independent monochromatic source of light. 

The Michelson interferometer and the two other mirrors of the displace- 
ment interferometer were mounted on independent cement piers which rested 
on the ground and were free from contact with the building. They extended 
up through the floor in the form of shafts one foot square and about three 
feet high. 

The illuminating system consisted of an arc light, a lens and a collimator 
all mounted on a table. To produce a zine spectrum for work on dispersion, 
the are was loaded with pieces of zinc wire inserted into drilled carbons. The 
arrangement was such that a concave mirror could be substituted for the lens 
to make it possible to focus sunlight from a hole in the floor of the room above, 
upon the slit of the collimating telescope to produce a solar spectrum. The 
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oval fringes were observed with a 25 cm focal length telescope having a gela- 
tin grating, with approximately 1150 lines to the inch, mounted normally in 
front of the objective. The arc light, the collimating telescope, the table upon 
which they rested, and the Michelson interferometer, were mounted upon 
squares of sponge rubber to minimize the effects of vibration. Even with the 
heavy cement piers completely separated from the building all trouble due to 
Vibration was not completely eliminated, although with them the fringes oc- 
casionally steadied down for short intervals. Ordinarily in the day time the 
quiet intervals were neither as long nor as regularly spaced as in the evening, 
and for this reason most of the readings were taken in the evening. 

The experimental refraction tube containing the gas to be experimented 
with, was supported in the working arm of the displacement interferometer. 
It was made of Shelby cold drawn steel tubing { inch wall and approximately 
one inch bore. The ends were closed by disks of high quality optical glass, 
obtained from Bausch and Lomb, 1.15 inches in diameter and 0.25 inches 
thick, with accurately plane parallel surfaces. A second tube just like the 
first but containing air at room temperature and pressure was used in the 


— = 
H ret 4e7d 
SLIITID LE MM ITD a Li EAL 
Fig. 4. Method of fastening glass ends to the experimental tube. 


other arm of the interferometer to optically compensate for the experimental 
tube. The length of the refraction tube inside the end plates measured 100.30 
cm. 

The packing shown in detail by Fig. 4+, was used to fasten the glass ends 
to the tube, and proved very satisfactory. The end of the tube 7° was care- 
fully tinned and, while still hot, was screwed into a hot threaded sleeve A 
to form a very tight fit. Cis a thin-walled steel cup carefully machined down 
to fit the glass disk G which is cemented into it with a thin layer K of de- 
Khotinsky cement. The steel washer W is also cemented to the glass, and the 
rim of C projects out beyond G just enough to take up the strain when W is 
pressed against it. Pisa very thin sheet of rubber (dental dam) also cemented 
to C. This whole unit P, C, K, G, W, was set in place and the ring R was 
screwed in about flush with A. Then by means of eight equally spaced screws 
S the unit was pressed tightly up against the end of the tube 7. 

The experimental tube was mounted in a water bath, by means of which 
the temperature was kept uniform and constant to 0.01° c. 

A small McLeod vacuum gauge was used for the purpose of registering 
pressures less than 0.1 mm mercury obtained with a Cenco Hyvac pump. 
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Pressures above atmospheric were measured with a mercury column. The 
lower end of this column was a length of 3. 16’’ steel tubing welded through a 
hole in the plug of a commercial iron mercury flask, and just long enough 
nearly to reach the bottom of the container. The upper end was a glass tube 
fastened to a brass centimeter scale, some 350 cm long in two sections, 
mounted on the side of a wooden ladder just long enough to reach from the 
floor to the ceiling beams, and fastened rigidly at each end. 

The pressure line connecting the various pieces of apparatus: the gas 
source, the mercury column, a dial gauge (which was included for rough 
work), and the experimental refraction tube was made up of {”’ steel tubing. 
Connections were made by means of soldered sleeves or carefully built steel 
unions. Each piece of apparatus could be isolated from the rest by special 
high pressure steel stopcocks."' 


EXPERIMENTAL PROCEDURE AND OBSERVATIONS WITH NITROGEN 


Water pumped nitrogen was obtained in commercial containers under 
about 2000 Ibs. pressure, guaranteed pure except for 0.5°, water and 0.5, 
oxygen. The oxygen was removed by passing the nitrogen through a steel 
tube full of sheet copper clippings maintained at a temperature of 500 600°C 
into an auxiliary steel tank where it was stored under a pressure of 600 Ibs. The 
water vapor was removed by phosphorus pentoxide contained in three tin 
cans suspended by copper wires inside this auxiliary tank. 

The nitrogen thus stored was allowed to stand a few days to become dry. 
The system was thoroughly rinsed out with it and observations were carried 
out in the following manner: A vacuum of 10° * mm of mercury was obtained, 
and after a reading of the barometer had been taken, the temperature of the 
bath being kept constant at 30°C, the fringes were very carefully centered 
on the blue (4811A) zine spectrum line. The reading of the micrometer was 
recorded, and then the screw was so turned as to cause the fringes to proceed 
through the spectrum to the vellow (5893A) sodium line. This process was 
continued to bring the fringes up to the orange (6362A) zine line. Then a 
small quantity of nitrogen was admitted to the system, and after the tempera- 
ture and pressure had come to equilibrium, the temperature of the room being 
usually several degrees~below 30°, the reading of the mercury column was 
noted, and readings of the micrometer were made with the fringes centered 
upon each of the three spectrum lines. This process was continued until 
readings had been made at intervals of approximately one atmosphere up to 
the limit of the gauge, i.e. a pressure of about 475 cm of mercury above 
vacuum. At this pressure due allowance was made for the back lash of the 
screw in the opposite direction, and readings were made for settings of fringes 
at each spectrum line, this time in the order of orange. vellow, and blue, for 
decreasing pressures at intervals of approximately one atmosphere, ending 
with a reading for vacuum. 

After several runs of this sort a very careful run was made at a tempera- 
ture of 0°C obtained by packing the constant temperature bath with ice. 


"FE. G, Keyes and Jane Dewey, Jour. Opt. Soc. Am., 14, 491 (1927). 
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Phe observations obtained in these runs are recorded in Tables | and Tl. 
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TABLE I. Data on nitrogen at 30°C. 
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TABLE I]. Data on nitrogen at °C. 
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COMPUTATIONS AND RESULTS 


The calculated values in the above tables were computed from the most 
probable values of Ax/p obtained from an analysis of the data by least 
squares. The observed values of Ax were also plotted against pressure on a 
paper 50 X90 cm of which Fig. 5 is a reduced photograph showing graphically 


“o 


Fig. 5. Experimental results for nitrogen. 


the linearity of the points. Values for By were computed by Eq. (11) for the 
ranges 4811A-5893A and 4811A-6362A at each temperature. 


Dispersion (By X 10") 


Range (A) ec 30°C 
4811-5893 1.90 2.08 
4811-6362 1.87 1.95 Mean 1.95 


Mean value of B (at 760 mm pressure) = 1.48 X10~-™. Although there is a 
very good agreement between the two values obtained at the zero tempera- 
ture, the precision of the measurements is not sufficient to detect by experi- 
ment any temperature effect. For this reason the arithmetic mean of all the 
B values was used in the calculation of refractive index by Eq. (6). The values 
of (u—1) thus obtained at three different wave-legnths were reduced to values 
corresponding to infinite wave-length, by the Cauchy relation 

B 
ia~- te - I) -—- (14) 
The values at 30° were reduced to zero degrees by Eq. (13) producing the 
results shown in the last row of Table III. The values in parentheses were 
obtained by reducing the 30° data to zero degrees for purposes of comparison 
with the data obtained directly at 0°. 


SIGNIFICANCE OF RESULTS AND CONCLUSIONS 


Observation of the graph Fig. 5, and a least square analysis of the data 
show that the displacements Ax follow very closely a linear relation with pres- 
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TABLE III. Calculated results. 








Temperature 0°C Temperature 30°C 
4811A 5893A 6362A 4811A 5893A 6362A 
Ax p 0 .0004122 0.0004040 0.0004018 | 0.0003743 0.0003653  0.0003635 
Prob. Error + .0000001 + .0000001 + .0000001 | + .00000005 + .00000003 + .00000003 
Ax /e -0003124 -0003061 .0003045 .0002836 .0002708 .0002755 
2B/ -.0000128 .0000085 -0000073 .0000128 .0000085 .0000073 
(u—1) .0002996 .0002976 .0002972 .0002708 -0002683 .0002682 
(.0002999) (.0002973) (.0002973) 
(u—1)« .0002932 -0002933 .0002935 .0002644 .0002640 .0002645 
(u—1)<.0 .0002932 -0002933 -0002935 -0002935 -0002930 0002936 
Mean .0002933 





+ .0000002 





sure. The origin lies just as well on the lines as any of the observed points. 
This linearity means, as was pointed out before, that the Lorentz-Lorenz rela- 
tion holds, and that dispersion is a linear function of pressure for nitrogen 
up to seven atmospheres. Theory suggests an inverse linear dependence of 
dispersion on temperature although the data do not indicate it. The values 
of By differ from the mean on the average by only 6°¢ whereas an inverse 
linear temperature function would call for values at 0° to differ from values 
at 30° by 107. In this connection it might be noted that Meggers and Peters 
data Fig. 2 show no consistent temperature effect. 

On the other hand the precision of the measurements is not great enough 
to guarantee By to better than 87. The micrometer settings can be made to 
one third of one percent. At one atmosphere Ax is accurate to this degree, 
but as the pressure increases, so does the accuracy. Analysis shows the errors 
in the slopes to be less than one part in four thousand, yet due to the fact 
that By depends on the difference between two slopes observed under differ- 
ent circumstances it is probably not accurate to better than approximately 
8¢;. An error of 8°7, however, in By corresponds to an error of less than 0.36 
in (u—1). On the other hand a survey of the literature indicates that Bo for 
a gas has not previously been recorded and that B has not been very accu- 
rately determined. 

This information on dispersion enables one to reduce observations of re- 
fractive index at any wave-length to a value for infinite wave-length or zero 
frequency, and for this reason the results are probably of most value. It is, 
of course, the refractive index at zero frequency which should be compared 
by Maxwell’s law with the dielectric constant for a gas. It will be noted that 
in Table III the values of (u—1),,.» all check to a remarkable degree, thus 
further emphasizing the degree of precision of the observations. 

The above analysis of the data shows that the accidental errors in the re- 
sults amount to less than one third of one percent. The only constant error 
of any consequence which affects them is the possible error in the pitch of the 
micrometer screw, which hitherto has been assumed to be 0.5 mm. To elimi- 
nate this uncertainty a very careful calibration of the screw was made over 
the range actually used by counting Michelson fringes as before described. 
The mercury 5461A line, filtered by screens provided by Cooper Hewitt, was 
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used for a source, and fringes were counted over a 2 mm range. The results 
of ten separate counts over a whole turn (approximately 1827 fringes per turn) 
indicated that the irregularities of the screw were negligible and that the 
actual value of the pitch was 0.499 mm in this range, making pw,» equal to 
1.0002927 instead of 1.0002933, and limiting the absolute error to one third 
of one percent. Asa result pu? is equal to 1.000585 to be compared with Zahn’s” 
value of € equal to 1.000581 which he estimates to be accurate to one percent 
in (e—1). The difference between these values amounts to approximately 
two thirds of one percent in (e€—1). 

Incidentally it might be observed that the values for the index of refrac- 
tion at a pressure of one atmosphere and a temperature of zero degrees at the 
three wave-lengths are as follows, corrected for the screw error: 


Refractive Index of Nitrogen at 0°C and 760 mm Pressure 
r 4811A 5893A 6362A 
Me 1.0002991 1.0002969 1 .0002967 


These values are in accord with previously published results obtained by 
different methods. 


ACKNOWLEDGMENTS 


This problem was suggested by Professor F. G. Keyes, Director of the 
Research Laboratory of Physical Chemistry at the Massachusetts Institute 
of Technology, and at the time the work was started, also Consulting Pro- 
fessor of Physics at Brown University. To Professor Carl Barus and Pro- 
fessor C. W. Miller, the writer is also indebted for assistance in various 
matters, particularly with reference to the displacement interferometer. 
The writer is, however, especially grateful to Professor A. deF. Palmer for his 
sustained interest and help in supervising the whole problem. 


2 C, T. Zahn, (International Critical Tables, vol. 6, p. 78.); also Phys. Rev. 24, 400 (1924). 








FEBRUARY 1, 1931 PHYSICAL REVIEW VOLUME 37 


ON THE ULTRAVIOLET ABSORPTION SPECTRUM OF 
ACETYLENE* 


By GrorGE B. KistIAKOWSKY 
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ABSTRACT 
It is shown that the ultraviolet absorption bands reported previously do not be- 
long to acetylene but to some impurity. In pure acetylene, bands are observed only be- 

low 2400A. These bands are arranged into three progressions and their interpretation 

is proposed using information from the work on the infrared absorption spectrum of 

acetylene. 

HE more recent investigations of the infrared absorption spectrum of a- 

cetylene have revealed a very simple structure of these bands! leading to 
the conclusion that the acetylene molecule has all four atoms in a straight line 
and is therefore characterised by rotational quantum levels typical of a dia- 
tomic molecule. It seemed promising to undertake a study of the ultraviolet 
absorption spectrum of acetylene in the hope that here also a more complete 
analysis could be achieved than is possible with the majority of polyatomic 
molecules. Not much work has been done on acetylene spectrum. Stark and 
Lipp* investigated the absorption spectrum in gaseous acetylene and reported 
bands between 2200-1900A without giving their wave-lengths. On the other 
hand, Henri and Landau® studying gaseous acetvlene and using a 40 cm ab- 
sorption tube found numerous bands at longer wave-lengths which they di- 
vide into three groups: bands with fine structure from 3157 to 2872A; partly 
sharp, partly diffuse bands in the region 2960 to 2495A and baiids with fine 
structure between 2327 and 2236A. 

In the present work absorption tubes up to 2 meters length were used. As 
a light source a hydrogen discharge tube was employed and the photographs 
taken first with a medium size and later with a large Hilger El spectrograph. 
Acetylene was taken from a small “Prestolyte” tank and was purified* by 
passing it slowly through a series of five Milligan wash bottles containing 
water, chromic acid in sulfuric acid, mercuric chloride in hydrochloric acid, 
copper nitrate in nitric acid and sodium hydroxide solutions, respectively. 
It was finally dried over phosphorus pentoxide. 

Acetylene thus purified was found to be completely transparent down to 
2400A and even using plates of the most contrast no traces of the some 70 
bands reported by Henri and Landau could be detected. Exposures were also 
made with acetylene prepared from commercial calcium carbide and only 
crudely purified. Here a continuous absorption was found extending from 


* Contribution #18 from the Loomis Laboratory. 

1 Levin and Meyer, Journ. Opt. Soc. Am. 16, 137 (1928); Hedfeld and Mecke, Zeits. f. 
Physik 64, 151 (1930); Mecke, Ibid. 64, 173 (1930). 

? Stark and Lipp, Zeits. f. phys. Chem. 86, 36 (1913). 

3’ Landau, C. R. 156, 697 (1913). 

4 Beilstein. Handbuch d. Organischen Chemie, Vol. 1. 
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ULTRAVIOLET ABSORPTION OF ACETYLENE 


shorter wave-lengths up to about 3200A. It became evident that the spec- 
trum described by Henri and Landau is due to some impurity and a further 
search for it was abandoned. 

In the short ultraviolet a number of bands was observed. The wave num- 
bers of the band heads will be found in Table I. Most of these bands are 
extremely faint and, in addition, a continuous absorption sets in at about 


TABLE I. Absorption bands of acetylene. 








Ist progression 2nd progression 3d progression 
Estimated 


Wave number intensity Wave number Intensity Wave number Intensity 
42056 1 42961 1 44084 5 
44100, 
42156 1 43074), 5 44205 3 
43088 44220 
42433 1 43184 3 44304 3 
43199 44318 
42504 1 43290 1 44475 3 
426048 2 43441° 2 
426060 43455 
426098) 3 43502) 2 
42712 43518 
43674) 10 
43687 | 
43730) 10 
43744 


Unassigned 44441 (3) and 44514 (3) 


2300-2350A and becomes so strong below 2250A that not more than 5 per- 
cent of the light of this wave-length is transmitted by a 2 meter layer of acety- 
lene at 760 mm pressure. In view of this, no great precision can be claimed 
for the wave-length determinations, the error for the fainter bands being 
probably of the order of 5 cm~!. The bands are shaded to the red and have 
double heads with a spacing of about 14 cm~!. The resolution obtained was 
not sufficient to show rotational structure near the heads, but in the stronger 
bands (43674-43687 and 43730-43744 em~') evidence of structure was ob- 
served some distance away. Even in this spectral region the bands reported 
by Henri and Landau must be attributed to impurities for, with the excep- 
tion of 42961 cm, none of them coincides with the now observed ones within 
less than 20 cm~'. It is of course not excluded that the bands here reported 
are also to be ascribed to impurities. On the other hand every precaution 
exercised in purifying acetylene seemed not to effect them in the least. 

The bands recorded in Table I can be arranged into three progressions at 
a distance of 1032 and 1012 cm~! apart. These progressions can be described 
by the same equation, as will be best seen from Fig. 1 in which the wave- 
numbers of the shorter wave-length band heads (after subtracting from the 
second and third progressions 1032 and 2044 cm~, respectively) have been 
plotted against an arbitrary number v7. A quadratic expression in v is, how- 
ever, certainly insufficient to describe the progressions since they show a pro- 
nounced convergence for the lowest values of v, but for v=5 to 10 the bands 
become nearly equidistant. It is also unusual that both the first and the se- 
cond progression break off abruptly after the strongest band (v =10). A tenta- 
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tive interpretation of the observed regularities is possible using information 
obtained from the vibration-rotation bands of acetylene. Mecke' has pointed 
out that a normal acetylene molecule possesses three modes of vibration par- 
allel to the molecular axis, of which one is optically active in the infrared 
(vy, = 3312.83 cm), and two transversal vibrations of lower frequencies 
(0; = 729.27 and @ = 1328.5). It is quite likely that the three progressions here 
reported are due to transitions into successive quantum states of either of the 
longitudinal vibrations of the excited molecule—probably v.—while the bands 
within each progression are caused by transitions to successive quantum 
states of one of the transversal modes of vibration. Against this interpreta- 
tion speaks mainly the abnormal intensity distribution within each progres- 
sion, but perhaps this, as also the rather unusual law of binding force neces- 
sary to account for the observed spacing of the bands can be attributed to the 
transversal character of the vibrations involved. 








= 


Fig. 1. A plot of three progressions of the acetylene bands. 





Although a calculation of the fundamental frequencies of the excited mole- 
cule is impossible at present, it can be pointed out at least that they are prob- 
ably lower than those of the normal molecule, the bands being shaded to the 
red, thus indicating a larger moment of inertia and presumably weaker bind- 
ing force in the excited molecule. Assuming that the rotational structure of 
the acetylene bands is similar to that of diatomic molecules and that the dou- 
ble band heads are due to R and Q branches, one can estimate from their sepa- 
ration and the known moment of inertia of the normal molecule,> the moment 
of inertia of the excited one. One obtains thus 25.6+0.3X10-* gm cm? as 
compared with 23.509 X10-*° gm cm? of the normal molecule. 

The origin of the continuous absorption, mentioned in the beginning of 
this note, remains uncertain. It is likely, though, that this absorption is due 
to residual traces of impurities, because Stark and Lipp? were able to observe 
absorption bands down to 1900A. 

It isa pleasant duty to thank Mr. Alfred L. Loomis for his kind permission 
to work in the Loomis Laboratory and for the facilities put there at my dis- 
posal. 


5 Assuming that the head of the Q branch coincides with the zero line of the band, that the 
equation of the R branch is Av = (B’+B")K+(B’ —B")K? and that for its head K = —B’+B"/ 
2(B’ —B"). 
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SPACE CHARGE vs. IMAGE FORCE IN THERMIONIC EMISSION 


By R. S. BARTLETT AND A. T. WATERMAN 
SLOANE Puysics LABORATORY, YALE UNIVERSITY 


(Received December 7, 1930) 


ABSTRACT 

This paper is a joint report on work by the two authors, preliminary to further and 
more detailed publication on different phases of the problems here discussed. Atten- 
tion is called to certain deficiencies in the Schottky image force as an explanation of 
the thermionic work function, both on theoretical and experimental grounds. Evi- 
dence is advanced to show that space charge is more important than image force over 
most of the region in which the image force is thought to be valid. It is further shown 
that space charge with Fermi-Dirac statistics is able to account for observed phenom- 
ena of thermionic emission, including the effect of external fields. Numerical ex- 
amples indicate some of the consequences of this point of view. 


HE most commonly accepted explanation of the thermionic work func- 

tion is that which attributes the major part of it to the “image” force of 
the electron (Debye, Schottky, Langmuir). Schottky! gave an analysis to show 
that space charge and structure effects, in the layer outside that in which the 
surface atomic structure of the emitting body plays an important part, were 
negligible compared with the image force. To do this it was necessary to es- 
timate the electron density in this layer; this he calculated on the assumption 
that the image force law was correct, showing that under these conditions 
space charge and structure effects were negligible. 

On the other hand, if the assumption is made that the potential and there- 
by the charge density in the above-mentioned layer is determined by space 
charge instead of by image force considerations then it may be shown that 
the space charge is more important than the image force. Thus we must look 
further for evidence by which to choose between the two explanations. As a 
numerical illustration of the difference involved, a calculation of the field in- 
tensities from these two contrasting standpoints at various distances from 
a plane tungsten cathode at 2300°K gives the following results in the absence 
of external field. E; being the field intensity in v/cm computed from the im- 
age force, E, that computed for the space charge field, we have, at a distance 
of 10-7 cm from the cathode, E;=3.6X10°, E,=4.0X10°; at a distance 5X 
10-7 cm, £;=1.4X10°, E,=8.010°; at 10-° cm, E;=3.6X104, E, =4.0 
10°. Thus, for distances from the surface greater than about 10-7 cm, the 
field intensity computed from space charge considerations is greater than 
when computed by the image force, and this difference is accentuated as the 
distance from the surface increases. The concentration of electrons shows a 
similar behavior. The space charge treatment of course assumes a continuous 
distribution of electricity whereas the image force assumes the escape of only 


1 Schottky, Phys. Zeits. 15, 872 (1914). 
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one electron. As Schottky correctly points out, since in practice we find a sit- 
uation intermediate between these two view-points it is necessary to examine 
which of these two extremes is most nearly approached in practice. 

A reconsideration of this matter seems desirable, not only in view of the 
above fact, but also because the earlier work on the subject postulated a Max- 
well distribution law in the electron atmosphere. Since the advent of the Som- 
merfeld electron theory it becomes desirable to find whether the degeneracy 
of the electron gas at and presumably also near the emitting surface affects the 
analysis. Among the chief reasons for the general acceptance of the image 
force explanation have been (1) the observed constancy of the work function 
for a given pure metal under moderate applied accelerating fields, i.e. the 
thermionic saturation current, (2) for pure metals the fairly accurate verifi- 
cation of the variation of thermionic current with stronger applied fields 
(Schottky effect) as computed by Schottky, (3) the estimation of a reasonable 
value for the work function, with the opportunity of characteristic values for 
different materials. Its chief drawback at the present time appears to be in its 
application to emission from coated filaments, i.e. in failing to account for 
lack of saturation and for the Schottky etfect observed in these cases without 
complicated, rather ad hoc, assumptions. 

It is clear therefore that before any superiority may be pointed out a space 
charge analysis should also lead to confirmation with experiment in the above 
directions. In deriving the thermionic current, extending through the Schott- 
ky effect, by use of space charge considerations alone, if Poisson’s equation 
is combined with a Maxwellian distribution of electron velocities the results 
show definite lack of agreement for pure metals. But if a Fermi distribution 
is substituted for the Maxwellian, the calculated dependence of current on 
temperature and applied field is found to be in reasonable agreement with 
experiment. This may be taken to indicate the constancy of the work func- 
tion calculated by this method, but it seems more reasonable to deal with 
directly observable quantities, and compare calculated with observed cur- 
rents under varying conditions. Furthermore, the variation of the current 
with the nature of the emitter can be explained more readily by space charge 
than by image force. According to the latter the field becomes constant with- 
in a certain critical distance from the surface, this distance being so chosen 
that the total work function thus calculated agrees with the experimental 
value. As justification for this step it is pointed out that close to the surface 
structure effects would be important. Now, assuming the space charge meth- 
od valid right down to the surface, different currents from different metals 
are explained by the different electron concentrations within the metals, and 
the work function comes out to be of the correct order of magnitude according 
to the Sommerfeld theory (i.e. = W,) without further assumption. This leads, 
incidentally, to a simple explanation of the observed interdependence of A 
and b, the constants in the usual thermionic emission equation. Moreover it 
appears that the use of Fermi-Dirac statistics takes account, in part at least, 
of structure effects. For a simple calculation will show that an electron gas 
is degenerate when the electrons are so closely crowded that the mutual po- 
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tential energy of neighboring pairs is large compared with their kinetic energy 
of thermal agitation, and obeys classical laws when the potential energy is 
small. For example, in an electron gas at 1600°K for which the Fermi A is 
10-*, i.e. the gas nearly in the classical state, the average kinetic energy 
is roughly twenty times the potential. 

Thus it appears that the space charge calculation based on Fermi-Dirac 
statistics is a more proper method of dealing with pure metals, and it may be 
further pointed out that the same calculation based on classical statistics, 
assuming a low electron concentration in the emitter, seems capable of ex- 
plaining results obtained with coated filaments. This latter point is being in- 
vestigated further. 

Following out their preliminary investigations on this subject® the authors 
have carried the work further, one of them taking exclusively the problem of 
deducing the thermionic current under applied accelerating fields from zero 
to fields operative in the Schottky effect, the other treating exclusively the 
case of the statistical equilibrium resulting under retarding applied fields (in- 
cluding zero). Both deal only with the case of an infinite plane emitting sur- 
face. 

The former (R. S. Bartlett) has handled the problem after the manner of 
Langmuir and Fry by postulating the emission of electrons with velocities 
distributed according to the Fermi statistics, these electrons being subject 
only to the space charge field, and obeying the law of continuity of current. 
This calculation is rendered troublesome by the necessity of graphical compu- 
tation of the Fermi integrals for the region of transition from a degenerate 
to a classical state, and thus are not at present in a form suitable for detailed 
report in this communication. The general results have been mentioned 
above, and appear encouraging, especially in pointing the way to a single 
explanation of thermionic currents from zero applied field up through fields 
of Schottky intensity. 

The latter author (A. T. Waterman) starting with Poisson's equation and 
the Fermi analogue A =A,e~"*'*? of the familiar Boltzmann relation, has 
evaluated the potential and the electron concentration at points outside the 
surface, where the atmosphere is in the classical state (assuming a degenerate 
state within the metal). In the presence of a retarding field FE the solution 
takes different forms on either side of a critical distance, which, expressed in 
terms of potential difference between this distance and the surface, is V’ = 
2 log E/B where 8? = 167(27m)*2h-3(kT)*”? e #7 where W is Sommerfeld’s W’;, 
the thermodynamic potential of the electrons within the body. Thus for V 
>V’ (thus including zero field), V (at distance x) = —2k7T/€ log (1+ €8/2kT 
x) and the number of electrons per cc, n=k7T/27é€x*®. For V<V', V=—kT 
/e(€Ex/kT —log 4E*/8*) and n= E?/2rkT e~*®*'*7, It will be noted that on 
the space charge conception the concentration of electrons in equilibrium 
with a hot body depends only upon the distance and the temperature,—in 
particular not upon the nature of the hot body, for distances at which the 
electron gas is in the classical state (roughly > 10-7? cm). Examination as to 


Abstracts 52, 53 Phys. Rev. 35, 668 (1930). 
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the above-mentioned critical distance introduced by an applied retarding 
field, which is a measure of the distance inside which the applied field has no 
appreciable effect on the electron atmosphere, results in the observation that 
applied retarding fields of the order existing in cold extraction currents would 
control the electron atmosphere extending into the region where the electron 
gas becomes degenerate very near the surface. A prediction of the work is 
that the ratio of the anode potentials which would give equal thermionic 
currents under zero applied field from two different emitters under otherwise 
identical conditions (temperature and geometry of tube) should be V;/V2= 
log B:/log B:. If Sommerfeld’s electron theory is assumed, then this ratio 
should show appreciable difference from unity, e.g. Ve./ Vw =0.65; on the 
classical theory the ratio should always be very nearly unity, e.g. Vea/ Vw = 
0.98 in the instance quoted. 
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ORIGIN OF THE AXIAL ROTATION OF THE SUN 


By Ross GuNN 
NAVAL RESEARCH LABORATORY, WASHINGTON, D. C, 
(Received December 24, 1930) 


ABSTRACT 


In earlier papers the solar atmosphere has been shown to rotate rapidly asa result 
of the electric and magnetic forces which act on the atmospheric ions. At the equator 
the ionized layers move much faster than the solar surface proper but this relative 
movement decreases with latitude in the manner required by observation. The fast 
moving atmospheric layers transfer momentum to the sun proper and its rotation is 
a necessary consequence of the observed atmospheric motions. The resultant torque 
is calculated and the angular acceleration shown to be adequate to account for the 
observed rotation. 


HE almost universal prevalence of rapid rotation in celestial objects has 

been pointed out again and again by astronomers, but the origin of such 
rotation has been somewhat of a mystery. Celestial mechanics shows how a 
star may rotate faster and faster as it contracts due to condensation, but 
gives no clue as to the origin of the initial rotation, which must be comparable 
to the final value. In an important paper Jeans! discussed the relation of the 
radiation of matter to the rotation of a star and he concluded that the radia- 
tion of the star carries off angular momentum in such a way that it “generally 
lessens the angular momentum per unit mass of the star.” He thus pointed 
out that the common assumption of a constant angular momentum for a ra- 
diating star is invalid. We must, therefore, for the most part, attribute the 
observed high rotational velocities to changes in the mean density together 
with an initial rotation of the star or else to some phenomena not adequately 
described by celestial mechanics. 

In a series of papers which deal with the electromagnetic effects operating 
in the solar atmosphere, it has been shown? that the anomalies of the sun’s 
apparent rotation arise from an electromagnetic drift of the ionized atmos- 
phere and are not directly dependent on the motion of the surface of the sun 
proper. It was shown that the solar atmosphere rotates in the same direction 
as the sun proper but much faster and that this superposed drift results from 
the interaction of the atmospheric ions with the observed solar magnetic field 
and a radial electric field. 

The sun may be thought of as a highly viscous mass of gas bound together 
by gravitational, radiative and electromagnetic forces in such a manner as to 
form a semi-rigid system. Radiative viscosity acts in such a manner as to 
equalize differences of rotation in different layers but Jeans! calculations indi- 
cate that these forces alone are incapable of equalizing the rotations in a time 


! Jeans, Monthly Notices R. A. S. 86, 328, 444 (1926). 
2 Gunn, Phys. Rev. 35, 635 (1930); 36, 1251 (1930). 
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so short as the age of the sun. The sun's observed magnetic field is asymmetric 
with respect to the axis of rotation and hence shells of highly conducting solar 
matter rotating about the geographical axis will induce electromotive forces 
and systems of electrical eddy currents which will strongly oppose motion 
relative to the magnetic field. A quantitative calculation of the electromag- 
netic damping forces due to axial slippage is difficult but we have seen that 
in regions of high electrical conductivity the forces due to radial contraction 
are tremendous’ and it seems probable that electrical forces contribute in an 
important manner to the apparent viscosity of the sun. To be sure that these 
forces have time to act we must examine whether the electrical time constant 
of the induced current circuit is sufficiently short compared to the solar age 
for current systems to build up. It is convenient to note that the largest time 
constant 7 possible is the constant for the whole sun which is given by 7= 
LY where L is the effective inductance and © is the total conductivity. On 
making this calculation it results that 7 is of the order of 10° years. This is 
very small compared to the age of the sun and we assume that inductive ef- 
fects are unimportant. 

On account of the high effective viscosity of the interior it seems satisfac- 
tory to assume as a first approximation that the sun is a semi-rigid body ro- 
tating on its axis at a period corresponding to the observed period of rotation 
of its magnetic pole. With the solar radius known it is evident that the periph- 
eral velocity of any point on the solar surface proper can be calculated. 
The difference between this velocity and the velocity of the same point as 
determined by an earth-bound observer watching the solar atmosphere is 
evidently the atmospheric drift velocity. This superposed drift velocity U 
of the ionized atmosphere is readily determined from observation or can be 
calculated from? 


u = ———— ~ (1) 


where £ and B are the electric and magnetic field intensities, \ the mean free 
path of the ions, and R the radius of the helix generated by an ion as it spirals 
around the impressed magnetic field. This radius is calculated from 


mV 


= 


Be Be 





(2mkT)*!? 
x (2) 


where m is the mass of the ion, V the component of velocity perpendicular 
to B and e the ionic charge. 

It is clear from Eq. (1) that the drift motion of the atmosphere drops ab- 
ruptly as the ion pressure increases to such a value that the free path be- 
comes of the order of R. Thus a more or less sharply defined transition layer 
exists between the atmosphere and the solar surface proper. In this critical 
transition region between the sun and its atmosphere there is a constant 


3 Gunn, Phys. Rev. 35, 107 (1930). 
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interchange of momentum, and the faster moving layers transfer their 
momentum to the slower. At the surface of the sun the temperature is low 
and the radiative and electromagnetic viscosities are small compared to the 
molecular viscosity, so that we need only consider the forces due to the latter. 
The tangential force applied to the periphery of the sun by the more rapidly 
rotating atmosphere is 
A 
F = s— (3) 
S 
where 7 is the mean coefficient of viscosity of the transition layer, « the mean 
difference in velocity of the layers, A the effective area in contact and S their 
separation. To a sufficient approximation we may assume that in the transi- 
tion layer the coefficient of viscosity is independent of the magnetic field and 
write 
(3ZmkT)*!* 


n (4) 


2(2)'!*x0" 

where #7 is the mass of the hydrogen atom, Z the mean atomic weight of the 
ions composing the layer, k the Boltzmann constant, 7 the absolute tempera- 
ture, and o the kinetic theory diameter of the ions. To a sufficient approxima- 
tion we may assume that the effective contact area A is a ring around the 
equator of width R where R is the radius of the sun and that the mean veloc- 
ity diflerence is not appreciably different from the equatorial value. The 
torque Q applied to the solar mass is 


0 =Ia= — —R? (5) 


where 7 is the effective moment of inertia and a is the angular acceleration. 
Assuming that the solar density is uniform we get from Eq. (5) 
5(3ZmkT)' * uR 
om. (6) 
2(2)'/*62 SM 
where ./ is the mass of the sun. Taking from earlier work Z =3.3, m=1.66 
X10-*% gm, k=1.37XK10-", 7 =6,000, ¢=10-8 cm, R=7X10" cm, MW =2 
X 1083 gm, w=5X10' cm/sec and S=5X10° cm we find a=2.3X10~-* rad/ 
sec.” 

We have seen in an earlier paper that the solar magnetic field is main- 
tained by a system of electrical currents which flow in rings about the sun's 
axis and that this system of currents is not primarily related to the rotation, 
but to the sun's radial and axial symmetry. This circumstance suggests that 
the magnetic field and hence motions of the solar atmosphere have not 
changed in order of magnitude for a considerable portion of the sun’s life, 
which various estimates put at 7 X10" years or 2X 10°' seconds. By multiply- 


* Gunn, Phys. Rev. 34, 335, 1621 (1929). 
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ing the angular acceleration @ by the life time we find to a crude approxima- 
tion that the final angular velocity of the sun should be 4.6 X 10>° radians per 
second, if its initial rotation was negligibly small. The observed angular ve- 
locity, as computed from the motion of the magnetic pole is 2.3 X 107° radians 
per second. The agreement is rather better than one might expect from a cal- 
culation based on steady states. 

Were allowance made for the increased radiation during the youthful 
period and a smaller moment of inertia used, due to the concentration of 
mass toward the center, then the elapsed time necessary to account for the 
observed rotation would be considerably reduced. This suggests that in the 
earlier stages of development the sun was unmagnetized and was a typical 
pulsating star.® It is clear that the solar rotation can be accounted for by the 
mechanical drag of the solar atmosphere on the sun proper if the velocity of 
the atmosphere exceeds that at the surface by 0.5 km. sec. This is the precise 
value required to explain the anomalous rotation of the sun; this having 
been adequately accounted for in an earlier paper by the author. The mech- 
anism of rotation which we have considered doubtless has universal appli- 
cation in all stars having a magnetic field. It seems probable that this mech- 
anism, which constantly supplies angular momentum to stellar systems, may 
assist in explaining certain puzzling facts in connection with spiral nebulae 
and binary systems. Numerical application to other star systems seems 
meaningless at the present time due to lack of data, and we can only point 
out that most other stars will behave like the sun. Ordinary mechanics has 
been found inadequate to explain the rapid rotation of heavenly bodies. The 
present investigation shows that the rotations can be accounted for by the 
momentum transferred to the body proper by its own highly ionized atmos- 
phere. We have shown previously that the motion of the solar atmosphere 
is controlled by the magnetic and electric fields existing in the atmospheric 
layers and that these fields are probably controlled by the rate of radiation 
of matter. We may then conclude that the rotation of all heavenly bodies is 
intimately related to the dissolution of matter and in a static, non-radiating 
universe rotating systems would not exist. 
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THE NATURE OF THE CONDUCTIVITY 
OF INSULATING OILS* 


By Karv F. HEeRzFELD 
DEPARTMENT OF Puysics, THE Jouns HopKINs UNIVERSITY 
(Received December 23, 1930) 


ABSTRACT 


It has been shown in a previous paper how the methods developed in the study of 
conduction of electricity through gases could be applied to determine the nature of the 
ions present in liquids and solids. In the present paper these methods have been 
applied to the measurements of the space charge set up in insulating oils by a steady 
electric field, which were recently performed by Whitehead and Marvin. The pro- 
cedure is as follows: the shape of the curve shows that the current must be consider- 
ably below three-fourths of the saturation current. It is possible to fit the curve closely 
by assuming the Langevin constant to be 8=36+10. From this the number of ions 
present in the normal state and their mobility can be calculated, which turns out to be 
3.5: 107° + 50°% cm,'sec per volt/cm. With the assumption that the mobility is in- 
versely proportional to the viscosity, the mobility of these ions in water would be 1.2 
10°*cm._ sec per volt ‘cm or in the usual units 12. If one assumes that Lorenz's results 
for organic ions can be applied these ions should have approximately 110 atoms. 
From the boiling point of the oil one would conclude that the molecule of the oil con- 
tains about 45 atoms. 


NUMBER of papers' have recently appeared from the Engineering 

School of The Johns Hopkins University, investigating the conductivity 
of transformer oils. The last of these make possible an application’ of Mie’s* 
formulas to discover the properties of the ions. 

We start with the space-charge distribution‘ in oil “B”, when a potential 
of 1500 volts is applied between large condenser plates 2 cm apart. The dis- 
tribution of potential is nearly symmetric around the center, proving that the 
mobility of the ions of opposite sign must be nearly equal. We are going to 
neglect whatever dissymmetry might be present. 

Fig. 1 shows that the line is very nearly straight in the neighborhood of 
the center, i.e. the field very closely constant for quite a distance. It follows 
from this that the current must be considerably below ¢ of the saturation cur- 
rent and the approximate formulas have to be used which apply, according 
to Mie, to weak currents. 


* Presented at The Third Annual Meeting and Conference, Committee on Electrical Insu- 
lation, Division of Engineering and Industrial Research, National Research Council, Bureau 
of Standards, November 7, 8, 1930. 

1 J. B. Whitehead and R. H. Marvin, A. I. E. E. Trans. 48, 299 (1929); J. B. Whitehead, 
Jl. Franklin Inst. (IV) 208, 453 (1929); J. B. Whitehead and R. H. Marvin, A. I. E. E. Trans. 
49, 641 (1930). 

2 K, F. Herzfeld, Phys. Rev. 34, 791 (1929), 

3G. Mie, Ann. d. Physik 13, 857 (1904). 

4 An enlarged drawing of this figure, which has been used for Fig. 1 of the present paper, 
has been kindly put at my disposal by Dr. J. B. Whitehead. 
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The next step is the calculation of 8(3°' is Langevin’s constant). For this 
purpose, one has to compare the field in the center and the field near the 
electrodes. The field in the center is 

ho = 450 + 10 volt/cm. 


The field near the electrodes is not measured directly; from the drawing it 
was estimated to be h,, =2400 volt ‘em. This would lead to 


hy,” 2400\2 
can PP! = -( = 28.5. 
hy? 450 


Subsequently, c=40 proved to be the value which suited the curve best, but 
calculations were also made with c=30 for comparison, which turn out to 
fit the curve appreciably less. 


c = 40 gives B = 36.3, (28)'* = 8.5. 


It turns out that it is more useful for the actual calculation to put Mie’s for- 
mulas into a form different from the one that had been given previously. 




















A Distance of center of test B 
electrode from A (cm) 


Fig. 1. 
Calling, as was done there, xx the location of the bend in the graph of the 


field, we measure the location in the condenser relative to xx and find for the 
flat (central) part® 


kh? = W?{1 + (e — 1)(e"+e%)] (1) 
(c — 1)\(e — 1 — (28)!'"*) fe — 1 — (28)! l—yx 
28 —c+1 c—1 l— xK 


in ve only /—x has to be replaced by /+x. 


® In formula (50) of the paper cited in footnote (2) 4/° should be stricken from the denom- 
inator. 
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With the high value of 8 given above the exponent is always so large that 
the deviation from a constant  =// for the field or from a straight line V= 
—IIx+ V)/2 for the potential is negligible in the range of this formula, that is 
within —x«K <x <x. That means /ip =//, there is no change in the amount of 
ions in the central part. In the central part Ohm’s law is preserved with the 
initial conductivity. In the range between the electrodes and xx, (1) holds, 
but there y, and ys have a different form from the one in the central part. yes 
as before is derived from y, by substituting ]-+ x for /—x. We have 


c— 1— (28)'"7 l— x 1 ¢—1-— (28)'? l—-x 
vi = + (28)! — (1 — - ~) (3) 
238 =>. _ + 1 1 — VK 2 28 l — XK 


But here too, with the large value of 8, the coefficients are so large that 
(c—1)e™ is appreciable only for small values of 1—x/l—x,, and than e~¥ 
with the argument /+.x/1/—.xx is negligible. With this in mind, it is useful to 
measure, as Mie does, the distance from the electrode instead of the center 
and write 


h? = he?[1 + (ce — 1)e] (1’) 


¢—1— (28)? x 1 ¢—1— (28)' x’ 
n= 2 « (1 +— ————-—}, (9) 
28 =. + 1 VK 2 28 XK 


The subsequent procedure is as follows: 

One draws a curve of / as function of x/xx, from x/x« =0 (electrode) to 
x/xx =1 (bend) and proceeds from thereon with =). The problem is to 
determine xx, or to put it differently, to determine J) xx. 

We accordingly integrate graphically beyond the point x/xx =1, taking 
into consideration the scale of our drawing and changing it continuously as 
we go on | because V=/hdx =xxfhd(x XK) |, until we get to a point where 
V = V,/2=750 volts. This must be the center. To show a few examples, we 
get with different assumptions for the position of the bend: 


YK 8 10.5 11 
l 12 12 1? 
V 680 750 765. 


It seems therefore safe to say that 


XK 7 
oan Os Oe. 
l § 


The points calculated with this value and c =40 are marked with crosses in the 
figure. Formula (53)* gives then 
J VK 


(23)!'*(28 — ¢ + 1) 
ces —-— = 0.2 + 0.05. 
joa lL (¢ — Ile — 1 — (28)'*) 


6 See footnote (2). 
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We put then the current 7=61.6X10-" amp/cm? with 6 probably near 
unity. That makes j,,,=68X10-"'+25°% amp/cm? or, if we consider that 
jsat iS the total amount of charge generated per cm’, that is in a cylinder of 
21=2 cm length, that gives the number of pairs of univalent ions, generated 
in 1 cc per second 


q = = 1.25 X 108-6 + 25%. 


Furthermore, we have 


SrK, 
B = —— = 40 + 10. 
ae 
If we want K in cm/sec/ volt/cm instead of per est. unit/cm, we have 
to divide by 300. It is this number that shall be meant by A from now on. 
We assume €= 2.5 + 10° 


— = 2.8 X 107 + 35%. 

Q 
Now ¢g/a=n¢ (my being the equilibrium number of the ions of each sign which 
is present in the central part of the condenser even under the conditions of 
our experiment). Therefore 


es ’ ——o 

— = Kn,* = 3.5 X 10" + 50%. 

a 
On the other hand, the conductivity in the central part, when Ohm's law is 
valid gives 


——— = aK = 1.1 X 10% 
450.2.e 
or 
no = 3.2 X 10! + 50% 
and 


K = 3.5 X 10-% + 50% in cm/sec /volt/cem. 


7 The current density was measured on March 27 and found to be 1.6X10™ amp. cm?; the 
space charge curve was taken April 4-6, a remeasurement of current on May 2 gave 2.9X10-" 
amp/cm?, the curves in Fig. 12 (in the last paper mentioned in footnote 1) give for the initial 
conductivity 5.1X10-" ho which, for 450 volt/cm, would amount to a current of 2.25 
X10-" amp/cm?*. But during April a considerable manipulation with the oil makes the first 
given current more trustworthy. There is another fact pointing in this direction. The ratio of 
the initial field to the central field in the steady state should be equal to the ratio of the initial 
current to the steady current, if the conductivity (number of ions) remains unchanged in the 
centra! part. But in the later measurements the ratio of the currents is only 1.45 (private com- 
munication), while the ratio of the fields is 750/450 =1.67. This is clear if in the meantime the 
number of ions has increased because then one would be further from saturation and the ratio 
of the two fields at the time of the second current measurements would be less. 

On the other hand, the measurement on March 27 (Fig. 11) points toa ratio of currents 1.7. 
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To get a reasonable comparison for this value, one remembers that ac- 
cording to Walden,* the mobility of large ions in different solvents is inversely 
proportional to their internal friction 7. Professor J. C. Hubbard, of this 
Laboratory, was kind enough to determine 7 for this oil and found it 0.35 
(water 0.01), which would make the mobility in water 1.2 X 10-4 cm/sec/volt/ 
cm or, in the usual units (multiplication with 96540 coulomb) equal to 12.+ 
50°. 

According to Lorenz*® one could now estimate that an anion of mobility 
12 contains about 110 atoms, an ion of mobility 18, 65 atoms. 

M. F. Hamburger, of the Department of Electrical Engineering, has 
kindly measured the boiling point of the oil and found it to be about 290°C. 
Accordingly, one would estimate that the oil has molecules containing about 
15 C-atoms, therefore about 45 atoms. If the ions are products of oxidation 
(last paper footnote 1) the anion should contain about two molecules (be the 
anion of an acid coming from the oxidation of a polymer of double molecular 
weight). The equal mobility of the action would have to be explained as due 
to the fact that the H®* ion has two oil molecules (or a molecule of the above 
mentioned isomer) attached to it. 

If the interpretation of Dr. Whitehead, that the increase in conductivity 
above the Ohm-value upon reversal is due to an ion of a mobility K =10-4 
cem/sec/volt/cm is correct, this ion would turn out, with our present data, 
to be very quick. Multiplying namely with 10° (units!) and 35 (ratio of of 
oil and water) we get 350 which is a speed only attained by H* in water. We 
might have to divide this by 2 or 3 as small anorganic ions move in water 
slower as expected, and reduce the estimated K to 1/2 10-*. Then we could 
still have the mobility of a small anorganic ion. 

If this is accepted, then the picture is as follows: There are present a large 
number of slow, organic ions which are continuously generated and move on 
to the electrodes where they are discharged without appreciable accumula- 
tion They are responsible for the space charge and the long time phenomena. 

Besides, there are quick, small anorganic ions which are not discharged at 
the electrodes, but accumulate there, without practically influencing the 
space charge, and are the cause of short time processes. Offhand, three origins 
could be imagined: (a) They are present as impurities, but are not supple- 
mented continuously. That is impossible, because then upon reversal of the 
field the current would not go beyond the initial current of the original field. 
(hb) They are generated continuously; then their number must be much 
smaller than the number of the slow ions, because the latter determine the 
space charge (8). (c) Their origin could be a secondary process on the elec- 
trodes (OH’?). 


§ P. Walden, Zeits. f. anorg. Chem, 113, 85 (1920), 
* R. Lorenz, Raumerfiillung u. lonenbeweglichkeit, Leipzig, (1922) p. 79. 
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POTENTIAL DUE TO A BURIED SPHERE 


By J. H. Wess 
WILLIAMS COLLEGE 


(Received December 16, 1930) 


ABSTRACT 


A series solution of the problem dealing with the electrical potential due to a 
buried conducting sphere is given. The general method of attack consists in finding a 
solution of Laplace’s equation V?&=0, subject to proper boundary conditions. The 
condition requiring the vanishing of the normal derivative of potential on the surface 
of the half space is satisfied by placing an image sphere in an upper half space. The in- 
troduction of this image sphere also introduces the complication of the interaction of 
the potential due to two spheres. This is handled by means of spherical harmonics and 
the solution of the problem is made possible by means of a special transformation in 
spherical harmonics. 


INTRODUCTION 
M’* Y problems of this type have arisen recently in connection with the 
subject of geophysics. Since the introduction of electrical methods for 
ore exploration, there have been numerous attempts to find some means by 
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Fig. 1. Buried and image sphere. 


which to predict the nature and scope of irregularities in homogeneity of 
structure beneath the surface of the earth, when point electrodes on the sur- 
face are maintained at definite potentials.! In most of these attempts, cer- 
1A. Petrowsky, Phil. Mag. [5] 31, 927 (1928); 28, 334 (1928). 
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tain approximations have been made. In fact, in the more complicated prob- 
lems, approximations must be made. The present problem of the sphere is 
worked out, making no approximations, as a first step toward finding out the 
allowable approximations which may be made in the more complicated prob- 
lems, e.g. a buried ellipsoid or cube. 

Assume a sphere of radius a and conductivity o, buried in an infinite half 
space of conductivity a2 (see Fig. 1). At the center of the sphere, locate the 
common origin of a set of cartesian coordinates x, y, s and a set of polar coor- 
dinates p,6,6. Place an electrode at the common point P (h/2, 0-c) and Po 
(po, 9%, do). Note that @)=7/2. The potential due to the electrode alone is 
I/4zroz in which r is the distance from Po to any point P and J is the current. 
It is desired to find the potential at every point in space due to the buried 
sphere S. Designate by « and U the inside and outside potential, respective- 
ly, due to the sphere alone. Then the functions u and U must satisfy La- 
place’s equation 


v°*@ = 0, (1) 


at all points in space, and are further subject to the conditions: 





ou OU . ou, 
-— —-— =-— (k-1)— > p=, (2) 
Op Op Op 
a= U, p=a, (3) 
U = 0, p= 2 (4) 
and 
Om 
—=0, += h/2, (5) 
Ox 


in which U, denotes the potential due to the electrode and k represents the 
ratio of conductivities ¢,/¢. Condition (2) results from the well-known equa- 
tion: 

ab a 


Pintle ete 
On, ON» 





in which ® is the total potential and 1 and n. are the normals pointing into 
the regions whose conductivity is a; and o2 respectively. In the above equa- 
tion 


u+U,.= ®, 
inside the sphere, while 

U+U,.=© 
outside the sphere. Then since 

0 0 
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while 


the last equation becomes 


Ou al, ou al’, 
— if + +i +> — = 0, 

Op dp Op Op 
from which Eq. (2) follows at once. The complication which arises in the 
solution of the above problem is due to the boundary condition (5). In order 
to satisfy this condition, an image sphere S; is placed in the upper half space, 
which is imagined filled with material of conductivity o.. Then, from sym- 
metry, the normal derivative of potential will vanish on the boundary 


h 


) 


It may be seen that the problem can be viewed as either one in which a single 
sphere S is situated in an infinite half space, subject to the boundary condi- 
tion (5), or a problem in which two spheres S and S;, (S; being the image in 
plane JN of S), are situated in an (whole) space. In the problem which in- 
cludes the image sphere, the boundary condition (5) is automatically satisfied 
while at the same time, the complication arises that the interaction between 
the two spheres must be taken into account. Analytically, this results from 
the fact that, when the image sphere is included, the applied potential con- 
sists not only of U,., the potential due to the electrode, but also U;, the po- 
tential due to the image sphere. This new potential will be called U,. Accord- 
ing to the uniqueness theorem, the solution of the problem with the image 
sphere will also be a solution of the problem in hand, since such a solution 
satisfies Laplace’s equation and boundary conditions (2), (3), (4) and (5). 
The solution of the problem, using the image sphere is the one here given. 
The following solutions of Laplace’s equation may be written immediately 


x p m ’ 
Mu, = > ( ) V,(u, >), (0) 


m—U a 


*) a m+ 
%2= > (- ) Vn'(u, >), (7) 


m=U p 


in which Y,,(u,@) stands for the ordinary spherical harmonic and yu for cos @. 
Because of condition (4) and because the potential of the sphere must be finite 
at p=0, m and w#: are chosen as u and U’, the inside and outside potentials 
respectively. If the coefficients A.,,,, and B,,,, in the above expressions, can be 
so chosen that «and U satisfy all the boundary conditions, the solution of the 
problem will have been found. 

Eqs. (6) and (7) satisfy Laplace’s equation, condition (4), and will satisfy 
condition (3) provided the coefficients in (6) and (7) be set equal, making 
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Vin(u, >) = V»'(u, o). 


The next step is to apply condition (2). Before this is possible, a suitable 
expression for U, must be found. U, consists of two parts 


U,=U.+ Ui, (8) 


a part U, due to the electrode and a part U’; due to the image sphere S;. From 
symmetry, the potential at any point P due to sphere S; is the same as the 
potential at point P; due to sphere S, where? P; is the image of P in the plane 
MN, i., e., 


a ( m+1 
U(P) = U(P) = >> ( “) Vm (midi) « (9) 


p: 


m= 


Thus the expression for the potential U,, at P, can be written, 


; T x p m : * 
Ue=—— F(Z) tal 60,0, 0) + D ( 
Po 


drpooe m=0 


a m+1 
) Y m(mi, oi), (10) 
Pi 


m=vU 


in which L,, is a Laplacian, a special form of the spherical harmonic Y,, Since 
¢=¢;, ¢; may be replaced by ¢. Hereafter D will be written for the constant 
term J/4702. The first member of Eq. (10) is the part of the potential due 
to the electrode, and the second member is the part due to the image sphere 
S;. The coefficients in the second member of (10) are identical with the coeffi- 
cients in the expressions (6) and (7) for the potential of the sphere S. The 
boundary condition (2) will be employed for determining the values of these 
coefficients which appear in (6), (7) and (10). This will be accomplished by 
substituting for «, U and U, their values from (6), (7) and (10) in Eq. (2) and 
equating coefficients of like tesseral harmonics. Before this is possible, how- 
ever, U’, must be expressed as a sum of tesseral harmonics with arguments yu 
and @ whose coefficients contain (besides the constants A, and B,,,, to be de- 
termined) successive powers of p. The first term of (10) is already essentially 
in the desired form. The second term of (10) is at present, however, expressed 
in terms of p;, u; and @. It will be expressed, as desired, in terms of p, wand @ 
by the establishment of a relationship of the type: 


oo a m+1 oo 
> (<) Vin(wisd) = Dofmlo’)¥ m(u’, >), (11) 


m=O Pi m=0 


and subsequently replacing p’ and yw’ by their values in terms of p and wu. When 
the second term of (10) has been thus transformed, the solution of the prob- 
lem is then found at once by comparing on the two sides of Eq. (2) coeff- 
cients of like tesseral harmonics. 
Now a relationship between tesseral harmonics, which makes possible a 
transformation of type (11), has been found. This relationship is 
Pag) (= 1)*"* po" | & & re 


— p ; 
Se ae etal ln os ek Sie. tll 
p;"*' (m — n)! Oh=-* kh* Oun’* 22, Bett! 


2 The subscript 7 is used throughout to refer to image points. 
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in which the quantities 0’ and 6; are related a shown in Fig. 1. Differentiat- 
ing with respect to uw’ and /, the result is 
Pn" (ui) 7 (m + k)! p" 
———= }} : Pi"(w’) = Yeinm. (13) 
p,"t! ken (mM — n)(n + Rk)! mre ; 
For the special case n =0, (13) reduces to the form,” 


P m( ui) < (m+ hk)! p’ 


——= >> —Pr(u') = Deron. (14) 


pet so «60 miki mt 


Substituting (13) and (14) in the second member of (10) one obtains 


D va m r 


p 0 m=QQ p Q ” ( 


— pki oe COs NO on b, om sin Ng) Lau | 
By an interchange of indices m and k, one tinds 


D x m x 
U. = y 2 (*) Lm ( Mo, Po, M, 0) + >» | tee mn 


PO m=0 \Po 
x 
+ = 2 (Ane cos vd + BB, sin 16)etune | (15) 
eo 
? m A 


which is the desired result. It is now possible to apply the boundary Eq. (2). 
The values of vu, U and U, to be substituted in Eq. (2) have been given above. 
This boundary Eq. (2) involves the first derivative of uv, U and U’, with re- 
spect to p. The second member of (15) for U’,, gives the potential in terms of 
p’ and yw’ at point P;, due to the sphere S. From symmetry it may be seen, 
however, that if p’ and yw’ be replaced by p and yp, the second member of (15) 
will then give the potential at point P, due to the sphere S,. If now the values 
of u, U and U, be substituted in Eq. (2), p set equal to a, and coefficients of 
like tesseral harmonics equated, the following sets of linear equations result: 


a m x 
Gal an an Koa(—) one > GatmA nt (16) 
0 


p k=n 
fa m x _ 
a a = Kan(- : = DoanrmBrs (17) 
Po k=n 
(n=1,2-- ©), (m=n,n+1,--- @), 
(n= 0), (m,k=1,2---@), 


in which 





an = 


k k + 1) 1 . 
+ sy ae 
(k— 1) m(k —1) 


* See Whittaker and Watson, Modern Anal. 4th. Ed., p. 400, Ex. 7. 
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? 2D/a\™ (m — n)! 
Ram SS (= = nooP m"(uo), (19) 
po \po (m+n)! 
and 
(m + k)! o\ eter 
ankm = ————_——— (“) ' (20) 
(m+ nk —n)i\h 


It should be noted that there results, for each value of , an infinite set of 
linear equations. The different sets are identical in form, however, and so it 
will be necessary to demonstrate the possibility of solution for only one of the 
sets. For simplicity, the set (16), for which 2 =0, is chosen, and it will be 
shown by an argument due to von Koch‘ that this system of equations has 
a unique bounded solution for the Ao,;. In the argument, a, will be replaced 
by its minimum value unity. If the argument holds for this value of a@,,, it 
will hold for all possible values of a,,. (It may be noted that the dependence 
of a, upon m is only slight.) Eqs. (16) may then be written: 


a m 
dim- lot + dom-1ye2 + cits + (1 + Um A om Sila -x.(~) (21) 
m Po 
where, 
(m= 1,2--- oo), 
In order to change to von Koch's notation, replace the unknowns Aq - + + Aox 
by x1 °°: X,%, ete. and write Eq. (21) as 
a i 
GiX, + Aint, + aigvs +--+ + (1+ aia: ++ = x (+) , (22) 
Po 
in which 
(¢=1,2---@), 
Divide (22) by (1+a;,) and set 
dik 
1 + d;; 
b:; = 0, 


kK; (<) 
eum — oom = Cj. 
1+a;; Po 


(i+ k)!fa\it* 
bir < - es . (<) , 
i!k! h 
; 6" (24) 
Cj <. K( -), 
Po 


*H. von Koch Jahres bericht der Deutschen Mathematiker Vereinigung 22, 285-291 
(1913). 


— = Dix, (23) 


Then 
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in which K is the upper limit of K;, namely D/po. Eqs. (22) become, 


=e YD bine = Ci, (25) 


Let 


substitute and obtain 


L x 
(= » £22 . | . 
x4 = Do biere : + dines —_ ¢;“ 
k=1 k=1 


Let 

x, = x,.¢ L o,, 
and so on, continuing this process of successive approximations indefinitely. 
Using the notation 


(2) 


bs. 


Yo; bux 
v=1 


xe 
(3) (2) 
bix > dis Dork, etc., 
t=] 


then according to von Koch's argument, the series 


¥; = C; 4- C; 1 a c;° oe ea 


(3) 


nua Tha + Thea Tila +:--, (26) 
k=1 k=1 k=1 


is absolutely convergent and represents the unique bounded solution of the 
system (25), provided }> 2, /b;,| <1 for all values of i and provided c;<K 
a fixed constant. The first of these conditions is fulfilled, provided a/h <1/2, 
as may be seen from the following summation, 


— * st) ( d y" (<)" | 
by. = —---——. { — = —-— —{— , 27 
2» 2d itk! h | h—a h ai 


The second condition is fulfilled, as from (24) 


; a\' 
ci <A ( ) 
Po 


where 


















POTENTIAL DUE TO A BURIED SPHERE 


and 


D 
K=s—.- 
Po 


The first five equations of the system (22) for the special case 
ad 1 ad 1 : z 
= ‘ — . A i = K ° 
h 2 Po 2 


will be written down; 





1.1250; + O.187x2 + .125x34+ .078x,+ .047x5--- = 4K | 
187x; + 1.187%. + .156x3+ .117%,+ .082x5--- = 3K | 
125x; ++ .156x2 + 1.156x, + .136x,+ .110x,--- = 1/8K | (28) 
O78, + ALTx2 + £13623 + 1.13604 + 123.05 ° = 1/16K | 
O47x1 + .082x2 + .109x3 + 11234, + 1.12345--- = 1/32K | 


The calculated values of x; and xs, using the first four terms of formula (26), 
are 


.3716K, 
1456K. 


v1 


Xo 


It is very labotious, however, to use formula (26) and it has been found that 
the values of the unknowns x; may be found more easily by solving the first 
of the Eqs. (28) for x, retaining only the first term; using this value in the 
second equation to obtain an approximate value of x2 by retaining the first 
two terms only; then correcting the value of x, by substituting the value 
found for x2 in the first equation and so on. By this method, using three of 
the equations in three of the unknowns, it is found that 


x, = .3725K, 
vo = .1458A, 
v3 = 0475 XK , 


in good agreement with the values as calculated by the formula (26). It 
should be noted that the values of the coefficients fall off rapidly in the pres- 
ent case. If the ratio of a/h<1/2 or a/py<1/2, the convergence will be still 
more marked. 

All that has been said in regard to Eqs. (16) may be said of the infinite 
set of Eqs. (17). The convergence, as » grows, becomes more rapid. For m 
greater than five, the A, and B,,; can be found to an accuracy of about two 
percent by making the approximation in (16) and (17) 


1 ; a m 
Aan =e Ran er ; 
Am Po 
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1 : ei” 
| ae = - ee . 
Om Po/ 


For » less than five, and a/h =1/2, the values of the coefficients 1, and B,,, 


can be found with an accuracy of the order of one percent by solving the first 
five equations, in the first five unknowns, of the infinite sets of Eqs. (16) and 
(17). 

The expression for the total potential, due to the buried sphere and the 
electrode, at all points without the buried sphere is 


i) a m+1 = a m+1 
Pp = 3 ( -) Y »(u, o) ot ry (=) ef u’, @) 
p 


m=0 p m=0 


(29) 


D = /p\™ 
aloes +o ( -) Lem(Mo, Qo, KM, o), 


PO m=0 Po 


in which the coefficients in each of the spherical harmonics Y,,(u, ¢) and 
Y,(u’, 6) have the values found by solving the sets of linear Eqs. (16) and 
(17) for given values of a/h and a/po. 

The first and third members of Eq (29). are expressed in terms of 9, 6 and 
@. The second member is expressed in terms of p’, 0’ and @. If actually there 
existed the two spheres S and S;, then the first and second members of Eq. 
(29) would express the potential due to the spheres S and S; respectively. If, 
as is the case here, there is only the buried sphere S in the infinite half space, 
which is bounded by the plane /.V, then the first two members of (29) are 
necessary for expressing the potential due to the buried sphere alone. The 
second member of (29) simply serves to accomodate the boundary condition 
expressed by Eq. (5). 


NUMERICAL COMPUTATION 


As an example, the potential (®— U’,) will be computed for a point on the 
surface of the half space directly above the center of the sphere. The poten- 
tial will be computed for various depths of the sphere and for the ratio of 
conductivity, k, of the two media equal to ten. The calculation will be made 
first by means of approximation formulae described below, and finally by 
means of the exact formula (29). An electrode of radius a), at which a current 
I is being introduced, is placed at the point whose cartesian coordinates are 
(h/2, 0, (3)'a). The potential is calculated for the point whose cartesian co- 
ordinates are (h/2, 0, 0). It will be assumed that the current is removed at 
some electrode so far distant from the buried sphere, that its influence upon 
the sphere may be neglected. It should be mentioned, however, that if the 
second electrode were not so situated, the potential of the buried sphere, 
under the influence of this second electrode, could be computed by means of 
the formula (29) and this potential added to the potential calculated for the 
first electrode. 

First approximation. As a first order approximation to the potential, the 
field due to the electrode is considered uniform over the buried sphere and 















POTENTIAL DUE TO A BURIED SPHERE 301 





























equal to its value at the center of the sphere, and the interaction between the 
buried and image sphere is entirely neglected. The potential at all points 
without the buried sphere then is given by the expression 


a\?* 
ep — Ue = ( ) [1 cos 6 a (4 11 COS @ + By sin o) sin 6| 
p ; 

(30) 


1\2 
,* (<) [.1, cos 6’ + (Ay, cos @ + By sin @) sin 6], 


p 
in which the coefficients have the values 


Di aAi- k) a 
A, = — ——— —coshr, 


po (k + 2) Po 
and 


lin) D(1—k) a j COS Go| P 

-=— ————_ — > sin Oo. 

Bul po (k+ 2) polsin go! 

Second approximation. As a second order approximation to the potential, 
the field over the buried sphere due to the electrode will again be considered 
uniform but the first order cross effect due to the interaction of the two 
spheres will be taken into account. By first order cross effect is meant that 
the field, due to sphere S;, is considered uniform over the buried sphere S and 
equal to the value it has at the center of the sphere S. Such a formula can 
be obtained by applying boundary condition (2) and then dropping all terms 
containing (a py)? and (a/h)* and higher powers of these ratios. If this be 
carried out, the resulting expression for the potential is identical with Eq. 
(30), but the coefficients have the values: 

D (1 — k) a 
A, = — ———_________. — cosh, 
po [(k+2) + (a/h)*%(k —1)] po 
and 
-lin) D (1 — k) a {COS do) 


- gon sin 6. 
By, po [(k+ 2) + (a/h)(k —1)] polsin oo 


Third approximation. As a third order approximation, the field over the 
buried sphere due to the electrode will be considered uniform so that terms 
in (a/po)? and higher powers may be dropped. However, the field due to 
sphere S;, will no longer be considered uniform over sphere S. That is, terms 
in (a h)* and higher powers will be retained. The formula found under these 
conditions is: 


a a m+ x fa m+1 
?P ned ee <= p» ( ) V n(u, ?) + > “) Vn (u’, ?), 


m-=Q p m=O p 


in which the coefficients must be determined from the infinite sets of equa- 


tions of the type (22). Since terms in (a/po)* are to be neglected, the right 
sides of Eqs. (22) after the first are to be set equal to zero. 
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The values of the potential (®— U’,) as calculated by these three approxi- 
mate formulae and finally from the exact formula (29) are shown in Fig. 2. 
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Fig. 2. Potential curves. (a) first approximation; (b) secona 
approximation; (c) third approximation; (d) actual. 


To obtain the total potential ® simply add to the values on the curve, the 
constant term due to the electrode, 
U, = I/4no2(3)'*a = M/(3)"?. 
This last equation serves to define 7, which occurs in the data. In order that 
the data may be readily put into practical units, the following calculation is 
to be noted. The potential due to the electrode is 
U, = D/r, 
where ¢ is the distance from the electrode. If, in particular, r=a,, a, being 


the radius of the electrode, 
U. = T/Arooa, 
or 
I = (Ue)a dmooas. 


Therefore, J/ in the above data becomes 
M= a;/a(Ue)a, 


and since (Ue)a, is the potential of the electrode in volts, ./ will also be in 
volts. 

From the curves, it is clear that any of the approximate formulae are suf- 
ficiently accurate for computing the potential, due to the buried sphere, in 
cases where the depth of the center of the sphere is greater than three times 
the radius of the sphere. In the region for which h/2 is less than a, only the 
exact formula can be used with any degree of accuracy. 

In conclusion, I should like to express my sincere thanks to Professor 
Warren Weaver of the University of Wisconsin for his helpfulness and in- 
terest in the above problem. 
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MAGNETIC SUSCEPTIBILITIES AND IONIC 
MOMENTS IN THE PALLADIUM AND 
PLATINUM GROUPS 


By ALBERT N. GUTHRIE AND L. T, BouRLAND 
Puysics LABORATORY, UNIVERSITY OF ILLINOIS 


(Received December 15, 1930) 
ABSTRACT 


Susceptibility measurements have been made on Ru, Rh, Pd, Os, Ir, and Pt and 
some of their compounds over a range of temperatures by means of an improved Curie 
balance. 

Pd and Pt are found to follow approximately a Curie-Weiss law. Ru, Rh, Os, and 
Ir are found to have paramagnetic susceptibilities which increase with temperature. 

RuCl; and IrOQ: are found to follow a Curie-Weiss law which gives for the ion Ru® 
a moment of 9.8 Weiss magnetons per mole and for Ir‘ 13.8. The spectroscopically pre- 
dicted value for each ion is 29.4. The paramagnetic susceptibility of Rh,Os is found in- 
dependent of temperature while that of RuO, increases rapidly with temperature. 
RhCl; is found to have an unmeasureably small susceptibility over the temperature 
range investigated. IrCl; and PtCl, are diamagnetic with susceptibilities independ- 
ent of temperature. 


HE determination of ionic magnetic moments in the Pd and Pt groups, 

by susceptibility measurements on compounds of these elements, is of 
particular interest in connection with the Hund method! of predicting these 
moments from spectroscopic theory. Hund? has employed his method for 
the moments of the trivalent rare earth ions and obtained remarkable agree- 
ment with the measured values of several investigators. In the case of the 
ions of the iron group, the elements Sc to Ni of the periodic table, the method 
has been unsuccessful. Laporte and Sommerfeld* calculated the moments to 
be expected for the iron group ions in two cases, infinite separation and zero 
separation of the J-levels of the ground multiplet predicted by the Hund 
method. Several of the measured values are found to lie outside of these lim- 
its. Bose has pointed out that the measured moments of the iron group ions 
are best explained by assuming that they are due to the spin of the electrons 
alone, the orbital moments being ineffective. Stoner® has called attention to 
the striking difference in the structures of the rare earth and iron group ions, 
and the effect that this difference should be expected to have on the ions with 
the conditions under which their moments are measured. In the iron group 
ions the incompleted group of electrons which gives the ion its mechanical 
and magentic moments is that of highest total quantum number. In the rare 





1F. Hund, Zeits. f. Physik 33, 345 (1925). 

2 F. Hund, Zeits. f. Physik 33, 855 (1925). 

’Q,. Laporte and A. Sommerfeld, Zeits. f. Physik 40, 333 (1920). 
*D. M. Bose, Zeits. f. Physik 43, 864 (1927). 

® E. C. Stoner, Phil. Mag. 8, 250 (1929). 
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earth ions this group of electrons is not that of highest total quantum num- 
ber. The interaction between the paramagnetic ion and surrounding ions and 
molecules will presumably effect primarily the electrons in the group of high- 
est total quantum number; in the iron group ions the electrons responsible for 
the magnetic moment, in the rare earth ions electrons that are not responsible 
for the magnetic moment. In general this interaction is a strong orbital inter- 
action while the spin moments are left relatively free. In a qualitative way 
this is an explanation of the apparent ineffectiveness of the orbital moments 
in the iron group ions and their almost complete effectiveness in the rare 
earth ions. 

There exist few data on the susceptibilities of the elements of the Pd and 
Pt groups from which ionic moments can be calculated. Bose and Bhar* have 
made measurements at room temperature on some chlorides and complex 
salts of some of these elements and have calculated ionic magnetic moments 
assuming that the susceptibilities of the compounds follow the Curie law. 
Cabrera and Duperier? made measurements on RuCl, RhCls3, PdCle, OsCl., 
IrCl3, and PtCl over a range of temperatures from about 20°C to 100°C. 
They found that RuCl; approximately followed a Weiss law which gave a 
moment for Ru’ of 12.9 Weiss magnetons. The susceptibilities of the other 
chlorides were found to remain constant or increase with temperature. 

For this investigation the compounds RuCls, RuQ., RhCl;, RheOs, IrCls, 
IrQO., and PtCl. were obtained from the American Platinum Works, Newark, 
N.J. These compounds were prepared with special regard for purity by Dr. 
S. Streicher of their chemical laboratory. 

Through the kindness of the Bureau of Standards samples of the six me- 
tals Ru, Rh, Pd, Os, Ir, and Pt of very great purity were made available for 
a new determination of their susceptibilities. Several investigators have 
measured the susceptibilities of these metals but the results have been in 
very poor agreement. 


DESCRIPTION OF APPARATUS 


The susceptibility measurements were made by means of a Curie balance. 
In this apparatus the usual mechanically operated torsion head was replaced 
by an electrically operated one. The ease and speed of making measurements 
with this arrangement far surpasses that attainable with the usual mechani- 
cal types. In Fig. 1 is shown a sketch of the essential parts of the electrically 
operated torsion head. A coil of small copper wire, A, wound on a flat, hard 
rubber spool is suspended by a fine wire, B, with its plane parallel to a mag- 
netic field produced by two fixed coils, C, carrying a current. To the coil A 
is attached the suspension, D, which supports the balance arm and to which 
is applied the torque to overcome the force on the sample. A plane mirror, £, 
is attached to the coil A and the twist in the suspension D is read from it 
by means of a telescope and scale placed about one meter away. In order to 
fix the axis about which the coil A rotates suspension guides, F, are em- 


6° D, M. Bose and H. G. Bhar, Zeits. f. Physik 48, 716 (1928). 
7B. Cabrera and A. Duperier, C. R. 185, 414 (1927). 
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ployed. These guides consist of thin brass strips attached to the box which 
supports the balance and pierced with small holes through which the suspen- 
sions pass. The coil A is damped by means of vanes G suspended in stano- 
lax. A current, which is controlled so as to give the proper torque to the sus- 
pension, is carried through the coil A by means of long thin spirals of wire 
which do not hinder the coils free rotation. The other parts of the balance 
are of the usual design.”¢ 

The apparatus was calibrated by making measurements on distilled water, 
using as its susceptibility with respect to air at 20°C and 76 cm pressure the 
value —0.749X10~°. 
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Fig. 1. Sketch of the essential parts of the electrically operated torsion head. 


The inhomogeneous magnetic field was furnished by a Dubois electro- 
magnet with its pole pieces set at an angle of about twelve degrees with a line 
joining their centers. The field strengths obtainable were in the range 1400 
to 6400 gausses. 

The samples were heated by means of a nichrome resistance furnace. A 
junction of a Cu-constantan thermocouple was mounted so as to occupy the 
same relative position in the furnace as the sample on which measurements 
were being made. The other junction of the thermocouple was kept in a bath 
of ice and water. Temperatures were read by means of a D’Arsonval gal- 
vanometer in the thermocouple circuit. A galvanometer deflection of one mm 
corresponded to a temperature change of about 2.5°C. The thermocouple was 
calibrated by means of the boiling points of well-known substances. 


EXPERIMENTAL RESULTS 


In Tables I, II, III, are given the results of measurements on the vari- 
ous substances indicated. 

The susceptibility of Rh was found to decrease with increasing field 
strengths. At room temperature the change was from 1.60X10-* at 1420 
gausses to 1.20 10-® at 6390 gausses. This variation was assumed to be due 
to ferromagnetic impurities and the method first used by Honda® was used 
to calculate the susceptibility of Rh from these measurements. 


7@ The modification of the Curie balance as used by Foéx was considered but calculation 
showed that the torsion method could be made more sensitive for our measurements. 


§ K. Honda, Ann. d. Physik (4) 32, 1027 (1910). 
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x X 10° 
T°K Ru Rh* Pd RuQ, Rh.Os 
298 0.427 1.08 $.15 1.21 0.408 
333 431 1.09 4.79 1.24 409 
380 435 Bi 4.39 1.28 410 
433 443 te 4.03 B.de 403 
480 452 1.14 3.73 1.306 403 
523 457 1.15 3.52 1.40 402 
573 .400 1.16 3.27 1.43 401 
623 475 [87 3.05 1.46 398 
673 487 1.18 2.85 1.49 397 
723 .496 1.19 2.06 
* Calculated from the measurements by the method of Honda. 
TABLE II. 
; : : : x X108 
T°kK Os Ir Pt IrO, 
298 0.052 0.133 0.982 0.997 
348 .138 947 984 
398 .059 .141 .925 974 
473 .146 .876 971 
548 .005 .151 .831 .942 
623 - .159 795 920 
698 .070 .167 745 891 
TABLE ILI. 
x X 108 
T°K RuCl, RhCl; IrCl; PtClh 
298 7.21 0 .00 —Q.114 —0.151 
323 6.64 .00 
348 6.106 .00 —0.114 —0.149 
373 5.77 .00 
398 5.44 .00 —.149 
423. 5.18 00 
5 4.91 .00 


For comparison the available data on the specific susceptibilities of the 
metals Pd, Rh, Ru and Pt at 48°C and the investigators who obtained them 


are given in Table IV. 


TABLE IV. 





Investigator 





Curie 

Honda 

Owen 

Onnes and Oosterhuis 
Foéx 

Kopp 

Authors 


Pd Rh Ru Pt 
35.2 : aa a 

5.8 1.14 0.56 1.10 
5.2 1.08 0.43 0.90 
5.3 _ aioe i 
5.27 1.028 
5.29 ie , 
5.24 1.08 0.426 
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In Fig. 2 are plotted the values of 1/x for Pd, Pt, RuCl; and IrO, as a 
function of temperature. 
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Fig. 2. Values of 1x for Pd, Pt, RuCl; and IrO, plotted as a function of temperature. 


DISCUSSION OF RESULTS 


Pd and Pt are found to follow approximately a Curie-Weiss law, the de- 
viation being accounted for by an underlying diamagnetism independent of 
temperature. From purely formal theory, it would be expected that the dia- 
magnetic component of the susceptibility would be greater for Pt than for 


. Pd but the experimental results indicate the opposite. 
Ru, Rh, Os, and Ir are found to have paramagnetic susceptibilities which 
increase with temperature. This is, of course, inconsistent with Langevin’s 
theory. 


RuCl; and IrO, are found to follow a Curie-Weiss law. For the former the 
| equation is: 


2340 10-¢ , 
x = —— X 10°‘ c.g.s./gram 
T + 37 





and for the latter 
x = ——— X 10‘ c.g.s./gram. 


These results give for the magnetic moment of the Ru’ ion 9.8 Weiss magne- 
tons per mole, and for the Ir‘ ion 13.8. The predicted normal state for both 
ions is °S which gives 29.4 Weiss magnetons per mole for the moments. In the 
t iron group, good agreement between measured and predicted moments has 
, been found for all ions with predicted S normal states. In these groups, there 
is no agreement whatever and the behavior of some of the compounds is such 
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that no moments can be determined from measurements on them. If the 
susceptibility does not follow a Curie-Weiss law, there is no theory which will 
give the elementary magnetic moments. 

A deviation from the Curie-Weiss law is predicted for compounds con- 
taining a paramagnetic ion with other than an S normal state due to the in- 
crease with temperature of the number of excited ions with moments dif- 
fering from that of the ions in their normal state. With an S normal state 
there are no excited ions at ordinary temperatures. In these researches only 
the compounds containing the paramagnetic ion with a predicted S normal 
state are found to follow a Curie-Weiss law but the deviations of the other 
compounds are opposite to (excepting RuO:) and much greater than can be 
accounted for by the theory. 

The behavior of those compounds which do not follow a Curie-Weiss law 
is not due to a permanent chemical or physical change due to heating for at 
room temperature the susceptibilities are found to be the same after as be- 
fore heating. 

The disagreement between the measured and the spectroscopically pre- 
dicted ionic moment may be due to error in prediction of the normal state. 
Hund’s rule, combined with present ideas of the electronic arrangement in 
atoms and ions, predicts for Ru* and Ir‘ a ®S normal state. In the same way 
the prediction for the normal state of the Pd atom is a *F. The atomic spec- 
trum of Pd has been analysed by Beals, Beckert and Catalan, and McLennan 
and Smith.’ They all have concluded that the normal state for the Pd atom 
is a'\S which means that the Pd atom in its normal state has zero magnetic 
moment. That the moment of the Pd atom is zero in its normal state has 
been verified by Copley and Guthrie!’ by the molecular ray method. 

The question of the validity of Hund’s rule in the wave mechanics has 
been examined by Slater.'' He has found that the wave mechanics gives the 
same normal state as Hund’s rule although their predictions differ higher up 
in the spectrum. 

The interpretation of the magnetic moment determined from suscepti- 
bility measurements on compounds as the true moment of the paramagnetic 
ion seems not to be well founded. The ion in chemical combination, sur- 
rounded by other ions and molecules, must be a different thing from the ion 
with which the spectroscopist deals in the Zeeman effect. Therefore, in spite 
of some past successes it appears hopeless to expect always to find agreement 
between these two. 

The writers wish to express their appreciation to Professors E. H. Williams 
and Jakob Kunz for their encouragement and advice during the progress of 
this work. 


® See F. Hund, Linienspektren und Periodische System der Elemente (1927). 
) Not yet published. 
"7. C. Slater, Phys. Rev. 34, 1293 (1929). 
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CHANGES IN ELECTRICAL RESISTANCE DUE 
TO MAGNETISM AND HARDNESS 


By S. R. WILLIAMS AND RIcHARD A. SANDERSON 
DEPARTMENT OF Puysics, AMHERST COLLEGE 


(Received December 19, 1930) 


ABSTRACT 


A series of twenty-two rods or strips of nickel, to which had been imparted differ- 
ent degrees of hardness by cold rolling, had their resistance measured when subjected 
to a longitudinal magnetic field. There were eleven degrees of hardness, two samples 
for each degree. The longitudinal magnetic field increased the resistance in all cases. 

The different degrees of hardness did not vary the specific resistance of the nickel, 
but it did cause a large variation in the change of resistance due to a magnetic field. 
The measurements were carried out at a constant temperature of 35°C. 


[I< HIS summary to Part I of an extensive study of the electrical conductiv- 
ity in strong magnetic fields, Kapitza' says, “It has been shown that the 
physical change produced in a conductor by hardening and annealing has a 
strong influence on the phenomenon of change of resistance. The influence of 
the impurities is also very marked and was studied.” 

Kapitza has brought together as fine a collection of the elements for study 
in strong magnetic fields as has ever been assembled, but his greatest difficulty 
was in procuring these samples in a pure state. It has been pointed out? that 
very worthwhile results may be obtained by comparing magnetic phenomena 
in the same specimen, even though the purity of the sample may not be 
vouched for. 

Some years ago there was begun in this laboratory a study’ of a series of 
nickel rods or strips, twenty-two samples in all, with eleven degrees of hard- 
ness in the series. The rods were given different degrees of hardness by cold 
rolling to various degrees of reduction. The rods were all of the same original 
heat and had the same chemical analysis, which is herewith given: copper, 
0.16°%; nickel, 98.8%; iron, 0.56%; sulphur, 0.008%; silicon, 0.06°%; carbon, 
0.09%; manganese, 0.23%. 

While these rods could not be considered as pure nickel, yet they could be 
said to be identical specimens so far as their chemical constitution was con- 
cerned. It was, therefore, possible to study the effect of different degrees of 
hardness on the electrical conductivity of a definite material when magnet- 
ized longitudinally. 

The rods came in flat strips and were furnished through the courtesy of 
the International Nickel Company. There were, originally, eleven strips, 91 


! Kapitza, Proc. Roy. Soc. 119, 358 (1928); 123, 292 (1929); 123, 342 (1929). 
2 Williams, Phys. Rev. 34, 258 (1912); 4, 498 (1914). 
3 Williams, Trans. A. S. S. T. p. 885, 1926; Science 65, p. 306 (1927); 66, p. 358 (1927). 


309 








310 S. R. WILLIAMS AND R. A, SANDERSON 


cm long, 4.8 cm wide and 0.61 cm thick, all annealed to the same degree. Ten 
of these were then given the percentage cold reduction shown in Table I. 


TABLE I. 


Strip number . 3 + 5 6 7 8 9 10 11 
Percentage reduction 0 9.7 18.9 28.9 29.5 50.0 59.5 69.0 79.0 89.1 93.3 


Each one of these eleven strips was then cut in two, longitudinally, and 
given a common length of 58 cm. This made a total of twenty-two specimens 
in all. The cross-sectional dimensions of the completed strips are given in 
Table II. They are the averages of ten to twelve readings taken along the 
length of the rods with a micrometer caliper. The two sets of strips were num- 
bered; 


and 
by, Bo, Bs, ae, a By 


wherein the rods A and B with the same subscript formed one of the eleven 
original strips. 


Tasce IT. 
Cross-section of strips 
Strip number Series A Series B 
1 0.953 X0.601 0.952 X0.603 
2 .956X .552 .955X .549 
3 .955X .496 .952X .495 
+ 942 XK .435 .944X .434 
5 .945X .372 .943X .370 
6 .950X .306 .904X .306 
7 .900 XK .248 .901X .249 
8 .9560X% .193 .956X% .193 
9 .950X% .131 .957X .131 
10 .957X .069 .956X .069 


11 955 .043 .956 X 


.044 


The scleroscope hardness, (universal tup or hammer) was furnished by 
the International Nickel Company and is given in Table III. 


Taste III. 











Strip number 1 2 3 4 5 6 7 8 9 10 11 
Scleroscope hardness 17 27 31 35 40 43 45 47 50 51 49 





The resistance was measured between two fiducial marks 51.8 cm apart 
on each rod. The measurements were made by means of a Kelvin double 
bridge. The complete outfit was one manufactured by Leeds and Northrup 
Company and found to be most satisfactory. The bridge was made sensitive 
enough to indicate easily a change of resistance less than 0.000001 ohm. 
Higher sensitivity created more difficulty in making observations. 
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In Fig. 1 is shown the arrangement of apparatus. The longitudinal mag- 
netic field was obtained by means of a solenoid, S,one meter long. The internal 
diameter of free space was 5.5 cm. There were 6025 turns which gave a con- 
stant of 75.3 gausses per ampere, both by calculation and by experimental 
test. In order to keep the strips at a constant temperature while being tested, 
they were placed in an oil bath. The container for the oil consisted of a cop- 
per trough, 77, passing through the solenoid and insulated from it by a coat- 
ing of shellaced cloth. By means of brass gas-pipe fittings this trough was 
connected with a rectangle G, of ¢ inch brass pipe. The pump, P, kept the 
oil in continuous circulation around the rectangle. This pump consisted of 
a propeller blade inside of a tube and driven by the small motor, M. Trans- 
former oil was used for circulating in the trough. Inasmuch as heat was 
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Fig. 1. Arrangement of magnetizing solenoids and constant temperature bath. 


slowly given up to the system by the solenoid, it seemed best to keep the 
temperature of the bath somewhat above that of the room. A temperature 
of 35° C, was chosen as a suitable point from which the bath would always 
be cooling off. The moment it got below 35°C, the thermostat would turn 
on more gas to the burner, located under one end of trough, and keep the 
temperature up at 35°C. The bath did not vary over 0.2 degrees during the 
measurements. This kept the temperature sufficiently constant for making 
comparable resistance measurements. The temperature of the bath was 
determined by a sensitive, calibrated mercury-in-glass thermometer. Some- 
times it was found that manual control of the gas was helpful. The ther- 
mostat for controlling the gas supply at the burner, ZL, which heated the bath, 
was a simple mercury-in-glass device. All the values for resistances, given 
in this paper, are for the temperature of 35° C, unless otherwise stated. 

As found by other investigators,‘ the resistance of nickel is increased when 
magnetized longitudinally. This was true for all of the twenty-two specimens 
examined. The effect of hardening may best be studied in Fig. 2, where for 


‘W. E. Williams, Phil. Mag. 4, 430 (1902); 9, 77 (1905). 
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the field strength of 75.3 gausses the changes in resistance per ohm is plot- 
ted against hardness numbers. The effect of increasing the hardness is to de- 
crease the change in resistance produced by the magnetic field. Fig. 2 is also 




















Fig. 2. Change in resistance per ohm for the different nickel 
strips at a field strength of 75.3 gausses. 


of interest in comparison with a curve published in a paper® presenting the re- 
lation between the change in length due to a magnetic field and the hardness 
of this same set of strips. This curve is shown in Fig. 3 with a field strength 
of 58.8 gausses. 














‘ig. 3. The changes in length for the different nickel 
strips at a field strength of 58.8 gausses. 


The most rapid change in length and in resistance comes in the first few 
reductions. In some way or another the changes in resistance due to a mag- 
netic field are very closely related to the changes in length due to a magnetic 
field. There seems to be a common mechanism operating in both effects. 

In Fig. 4 is shown the way in which the change in resistance per ohm varies 
with increasing magnetic field. Curve 1 is the average of the values for the 
two rods, A;, and B,, while curve 2 is the average for the two rods, Ay and By. 
This figure represents the greatest difference obtained in the series of strips, 
due to differences in hardening. It confirms in a very marked way the re- 
sults obtained by Kapitza, i.e., that hardness plays a very important rdéle 
in this effect. 


® Williams, Trans. A. S. S. T. p. 885 (1926). 
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The various values of dr/R in the curves represent the average of at least 
eight observations. Demagnetization of the rods occurred between each 
reading. This was accomplished by sending a decreasing a.c. through the sole- 
noid. 

In getting the changes in resistance due to magnetism and hardness, it was 
also easy to find the specific resistance of these strips. The values obtained 
are given in Table IV, and which are the averages of the A and B series. Three 

















Fig. 4. Change in resistance per ohm as magnetizing force increases. 


sets of values are given which were taken bydifferent observers using different 
methods and about two years apart. The first set is by Mr. Lusk of the Raw- 
son Electrical Instrument Company. Mr. Lusk checked his results by two 
entirely different methods. The second set is by Mr. G. K. Schoepfle, a for- 
mer graduate student in this department. The third set is furnished by the 
junior author of this paper. 

This set of readings seemed rather astonishing at first because as one 
crushes the nickel strips by rolling, and therefore the crystals comprising 














TABLE IV, 
Strip number Specific resistance 

Lusk Schoeptie Sanderson 

1 0 .00001012 0 .0000 1083 0 .00001099 
2 .00000993 .00001081 .00001102 
3 .00000993 .00001035 .0000 1093 
4 .00000995 .00001085 .00001089 
5 .0000 1000 .0000 10060 .0000 1086 
6 .00000995 .00001071 .00001082 
7 .0000 1000 .00001055 .00001079 
8 .00001008 .00001070 .00001071 
9 .0000 1032 .00001090 .00001076 
10 .00001033 .00001051 .00001079 
11 .00001038 .00001102 .00001073 
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the strips, it would suggest the possibility that the specific resistance would 
also be changed in a corresponding manner. The specific resistance appears 
to be constant. 

The values given by Lusk were obtained about five years ago and at room 
temperature. The values by Schoepfle were obtained about two years later 
and also at room temperature. In all probability there had been some aging 
process going on in the meantime. To release local strains set up in the rods 
at the time they were cut and machined into shape, they were annealed for 
one hour in an oil bath at 304°C. This process took place about seven years 
before the readings were taken for this paper. 

It was also possible with this equipment to determine the temperature co- 
efficient of the strips. This was not worked out completely, but it was care- 
fully measured on the rods, A; and Ax, those rods which had the widest varia- 
tion of dr/R for a given magnetic field Measurements were made at 35°C and 
at 85°C. This gave A; a temperature coefficient of 0.0042 and A, the value of 
0.0043, indicating that there was little or no change in the temperature co- 
efficient due to difference in rolling. 

Although hardening by cold rolling did not change the specific resistance 
of this series of nickel strip, yet it did affect the change in resistance due 
to a magnetic field. This seems to be the most significant fact emerging from 
the results of this investigation,—no change in specific resistance due to 
varying hardness, but a decided change in the variation of resistance with 
magnetic field. 

Attention should be called to the excellent paper by McKenhan® which 
has appeared since these results were obtained. 


® McKeehan, Phys. Rev. 36, 948 (1930). 
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THE ROLES OF DISCRETE AND CONTINUOUS THEORIES 
IN PHYSICS* 
By ARTHUR RUARK 
UNIVERSITY OF PITTSBURGH 
(Received December 4, 1930) 


ABSTRACT 

The differential character of the principal equations of physics implies that phys- 
ical systems are governed by laws which operate with a precision beyond the limits of 
verification by experiment. This appears undesirable, from an axiomatic standpoint. 

Postulates of real, finite and rational description.—\WVe may safely assume that all 
results obtained experimentally can be described by using real and finite functions of 
real, finite and rational variables. It appears desirable, therefore, to use the calculus of 
finite differences in formulating the results of our experiments. In dynamics, the use of 
this calculus yields terms corresponding to the fact that the configuration is known 
only when we observe it and must be arbitrary at intermediate times. For example, if 
a “uniformly moving particle” is observed once a second, we find experimentally that 

A? x/Af =0, 
whence 
v = Ft+G. 

Here F and G are arbitrary Fourier series having the period At. The general difference 
equations of dynamics are formulated. The Lorentz transformation is reinterpreted, 
relativistic difference equations of motion and of light propagation are set up, and geo- 
desics in a discrete manifold are discussed. 

Extensions of these ideas to electrodynamics are outlined. 


I. INTRODUCTORY 


T IS often said that while the fundamental equations of physics must be 

capable of describing every feature of our experimental results, at the same 
time they must not introduce extraneous or unobservable features. On the 
basis of this criterion, many branches of physics are not in a satisfactory po- 
sition. Broadly speaking, the difficulties to which we refer originate in the 
use of mathematical concepts which are unsuited to the description of our 
experimental knowledge, and which introduce complications having no coun- 
terparts in nature. 

Let us begin with an example. Consider, on the one hand, the information 
we obtain by performing an experiment once, or at most, only a few times. 
Usually, this information consists of a table of related numerical data. Con- 
sider, on the other hand, the statistical knowledge we obtain by correlating 
the results of a very large number of similar experiments. We draw curves 
and calculate averages, and finally, we may attempt to represent our results 
in an idealized fashion, by writing empirical equations. Up to this point the 

* A preliminary account was given before the American Physical Society in February 1929. 
See Phys. Rev. 33, 637 (1929). 
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procedure is beyond criticism, if we keep in mind the experimental basis of our 
equations, but when we carry the idealization process one step further, and 
derive differential equations from the empirical ones, we introduce conceptual 
difficulties. A differential equation connotes much more than the data from 
which it arises; herein lies its strength and,—unfortunately,— its weakness. 
For instance, Newton’s equations imply that a dynamical system is governed 
by laws operating with a precision beyond the limits of detection by experi- 
ment. If, at a time ¢, we know the position and velocity of a body, these equa- 
tions predict its condition at the time ‘+d, but in order to verify the predic- 
tion, one would be obliged to distinguish between two positions, x and x+dx, 
by making two measurements separated by an infinitesimal time interval. 
But as a matter of fact, we have no desire to find the condition after the time 
dt; we wish to predict the state, within limits, after a finite time interval, and 
this we do by integrating the equations of motion. Each coordinate may then 
be expressed as a continuous function of the time, and we may calculate the 
configuration at any time, that is, at an “infinite number of instants.” The 
information yielded by the equations is much in excess of what we can test 
experimentally. No one will deny the great practical value of Newton’s equa- 
tions, and clearly, the point we wish to emphasize is this: from an aesthetic 
and logical standpoint, it is desirable to have available a mathematical al- 
gorithm which makes it possible for us to summarize precisely the informa- 
tion obtained from a series of measurements, without excess and without defi- 
cit. It is our aim to show how this may be accomplished with the aid of the 
calculus of finite differences. 


II. CRITIQUE OF THE USE OF INFINITESIMAL ANALYsIs IN PuysiIcs 
Tue PRINCIPLE OF REAL AND FINITE DESCRIPTION 


Most of the purposes we have in mind may be accomplished by adopting 
a set of postulates which, for brevity, will be collectively designated as the 
principle of real and finite description. For our present purpose! we may state 
the postulates as follows: 


All results obtained experimentally can be described in terms of real and finite 
functions of real and finite variables. These functions need only be defined for 
discrete values of their independent variables. Both the functions and their inde- 
pendent variables need only assume rational values. 

As an illustration, we consider a series of experiments on the motion of a 
single atom, or on a small number of atoms. The principle of real and finite 
description indicates clearly that our numerical results should not be dis- 
cussed by means of continuous analysis. But if the series of experiments were 
to be extended indefinitely, until observations had been made on a great many 


1 For the sake of simplicity, we shall speak as though a measurement consisted in assigning 
a definite number xo to a variable at a definite time, fp. The real situation is that we know x lies 
in the range x» + xo when / lies inthe range tp + Ato. Thus we may say that x is an “indefinite func- 
tion” of the “indefinite variable” ¢. It is easy to devise an algebra for dealing directly with the 
possible connections between such variables; but these matters will be postponed to a later 
occasion. 
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individuals, and over the entire ranges of the independent variables, then we 
would have the material necessary for building up a continuous, statistical 
theory,—in fact, a branch of wave mechanics. The situation in regard to the 
motion of a macroscopic body is quite similar. On the basis of a small amount 
of experience with its motion, we should not write a continuous law. (The 
true status of the laws of motion is connected with Ehrenfest’s theorem, and 
is discussed in section V.) In the light of these ideas, the use of continuous 
analysis in physics can be clearly understood. Continuous laws should be 
introduced only when we deal with statistical results, obtained from vast 
numbers of experiments. 

While the situation in regard to the use of infinitudes in physics is gener- 
ally unsatisfactory, an exception must be made in the case of wave mechanics, 
where it is required that the wave functions shall be finite, continuous and 
single-valued at all real points of the coordinate space used by Schriédinger. 
The validity of these conditions has sometimes been questioned on the ground 
that, where the wave function y is zero, 1/y is infinite and there does not 
seem to be any reason why we should keep W from becoming infinite, rather 
than 1,y. The writer feels that this criticism is invalid, for the square of the 
absolute value of Y is measurable, indirectly at least, and therefore must not 
become infinite for real, finite and rational values of its variables. How wy be- 
haves for irrational values of its variables is none of our concern. 

Again, it may be objected that in wave mechanics integrals extending over 
the entire ranges of the coordinates x, vy and z (from —* to+~) are con- 
stantly used with valuable results. Just as in the case above, this appears un- 
desirable since we should think of the coordinates as measurable quantities, 
and should therefore limit them to finite values. Since in many problems the 
¥-function approaches zero very rapidly for large values of the coordinates, 
it will usually occur that the values of matrix components will be little af- 
fected by restricting the range of integration to exclude infinite values, and 
the practical consequences of the theory will be unaltered. Even in scattering 
problems, where sinusoidal wave trains are used to represent incident beams 
of particles or protons, it would be better physics to speak of these trains as 
having their origin at a great, but finite, distance from the scattering center. 
Further, it appears desirable to use wave equations which are modified to take 
account of general relativity. It appears plausible from astronomical evidence 
that our four-dimensional space-time is limited in extent, and if this is the 
case, integrals of the type /y*fy dv would naturally be extended only over a 
finite range. Furthermore, from our present stand-point, it is obvious that 
we should not treat the electron or the proton as point-singularities at which 
the electromagnetic and gravitational forces become infinite. Indeed, all 
such artificial mathematical difficulties are excluded by strict adherence to 
our postulates. 

The situation in regard to the use of imaginaries in physics also requires 
comment. It is generally recognized that the introduction of exp iw? in vibra- 
tion problems is merely a mathematical artifice, but many a physicist would 
feel surprise at the proposal to replace the quantum condition, pg—qp= 
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h/2zi, by an equation which is free of the imaginary unit. However, from 
an axiomatic and aesthetic standpoint, imaginaries should be excluded from 
any discussion of fundamental physical questions, and this can always be 
done by rationalizing our equations, even when we are dealing with quantities 
which do not commute. (In the present instance we may use (py—gp)?= 
—h?/4x° or (pg—qgb)(bq—gp)* =h?/4x*. Then, in solving any specific prob- 
lem, we get the Heisenberg matrices and their imaginary conjugates. In 
practice, this doubling of the solution will introduce no difficulties, and may 
even be advantageous.) 

After these explanatory remarks we are ready to apply the principle of real 
and finite description to dynamics and electromagnetic theory. We shall find 
that when we apply it in dynamics we are naturally led to equations of finite 
differences, replacing the differential equations of motion. The solution of these 
difference equations brings in arbitrary functions which reduce to known con- 
stant values only at times when the configuration of the system is actually ob- 
served. This corresponds to the fact that we are not able to state the con- 
figuration of the system, except at the instants of observation. We proceed at 
once to amplify these statements. 


Ill. THe DIFFERENCE EQUATIONS OF DyNAMICS 


Let us analyze the information obtained by an experimental study of a 
body moving along the x-axis. For the moment we pass over difficulties aris- 
ing from errors of observation and assume that the observer can determine 
times and distances with any desired precision. The result of a set of position 
measurements will be a table showing the values of x at each instant of ob- 
servation, such as the following: 

hi, to, ts, ee » bn; 
(1) 


V1, Vo, ¥3,°°* » Une 


With a view to simplicity, suppose the intervals between successive ob- 
servations have a constant duration, At. This restriction is unnecessary, and 
later we shall remove it. Perhaps the most natural method of analyzing the 
data is to study the difference coefficients of the series of x-values. Writing 
Xe— x, =Ax,, A’x,; =Ax.—Ax, etc., we form the table (2). 





vy 
Ax, 
At 
A*x 1 
Xe 
Ai? 
Axe 
jeome (2) 
Al 
A°x 2 
X3 dete 
Al? 
A X3 
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X4 


Let us suppose we find that 


Ax; . 
— =constant, i=1,---,n2—1, (3) 
At 
and 
A®x; ) ' 
<= m i=1,---,n—2. (4) 


from our present standpoint, such a motion will be defined as uniform, and 
the constant value of Ax,/At is defined as the velocity. Either one of the rela- 
tions (3) and (4) enables us to construct the original table of x-values,? but 
in order to attain a degree of generality equivalent to that provided by the 
equation d°x/df =0 we shall deal with (4). Taken as they stand, the Eqs. (3) 
and (4) state a regularity of the data and say nothing about the values of 
Ax /At and of A’x/Af at instants which lie outside the period covered by the 
experiment. Further, they say nothing about the values of A’x/Af for times 
intermediate to the instants of observation. They contain nothing more than 
the isolated facts revealed by the experiment, and we must recognize the ele- 
ment of hypothesis introduced when we use them as a basis for prediction. 
The most natural method of generalizing (4) is to write 


A*x 
AP 





= 0, (5) 


and to consider this equation, provisionally, as applying to any series of ob- 
servations like the above, obtained at time intervals At and extending over a 
finite time. Now, the general solution of (5) is not identical with the general 
solution of d’x/di#?=0. In fact, it is 


x =Ft+G, (6) 


where F and G are trigonometric series of period Af,’ let us say 
7 : t ; t t = 
F = Ag + A, sin 2n— + Aosin 4r— + ---+ B, cos 2m— + <>, (7) 

Al At Ai 


the expression for G being similar. The coefficients in these series are arbi- 
trary constants, except for one restriction. We require that F shall reduce to 
the observed velocity v, and G to the observed initial value of x, let us say xo, 
whenever ¢ is a multiple of At. The coefficients in (7) are not Fourier coeffi- 
cients which cause the series to represent the function x)+v/, in the interval 

2 In using equation (3) to reconstruct the x-values, we must specify x at one of the instants 


of observation, let us say ¢,; and in using (4), we must specify one value of x and one of Ax/At, 
3 Boole, Finite Differences, reprint by Stechert, p. 80. 
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from {=0 to t=At. On the contrary, when an arbitrary choice of the form 
assumed by x between observations has been made, they are the Fourier co- 
efficients of that form, and the series will generally be infinite. However, the 
series may contain only a finite, but very great, number of terms, if we specify 
the values of x at only a finite number of instants. The latter type of choice 
is in fact the only type which is in agreement with the principle of real and 
finite description. With these understandings, Eq. (6) completely reproduces 
the experimental results. At intermediate times, when the system is not un- 
der observation, the arbitrary trigonometric terms come into play, which 
means that the position of the system is not known.* The state of affairs 
is illustrated by Fig. 1, in which the solid line represents the classical motion, 
while the dotted curve shows a path which might possibly be represented by 
(6) and (7). The dotted curve might even have a finite number of disconti- 
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Fig. 1. 


nuities in each interval At. When the specified set of observations is com- 
pleted, we have no physical evidence which excludes the possibility of such 
a curve as the dotted one in Fig. 1. Of course, we could have taken addi- 
tional observations, at the times Af/2, 3A#/2, etc., let us say, obtaining an 
equation similar to (6). Under these circumstances, however, the state of 
our knowledge is different, for now the periodic function x must have the 
same values as the function vf+x 9, whenever fis a multiple of one-half the 
original Af; at all other times it is arbitrary. The number of observations at 
times between ¢; and /¢, cannot be increased indefinitely; therefore the Eqs. 
(5), (6) and (7), since they yield the observed data only and nothing more, 
represent a closer approximation to our experimental knowledge of the mo- 
tion than the corresponding differential equation and its solution. However, 


‘ Ina paper entitled “The Limits of Accuracy of Physical Measurements,” Proc. Nat. Acad. 
Sci. 14, 322 (1928), the writer has shown that there are definite natural limitations on the accu- 
racy with which we can measure a coordinate or specify a time, quite independent of all con- 
siderations of experimental skill. If we are dealing with a free particle of mass m, we have 
AXmin =h/me and Atmin=h/mc*. It is interesting to note that if the time interval between ob- 
servations is reduced to the order of Atmin, the first trigonometric term in (7) has a frequency of 
the order mc?/h; that is, the frequency of the matter waves for a particle of mass m. 
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(6) still suffers from several defects: (1) it definitely prescribes the value of x 
at instants lying outside the range ¢, to ¢t,, although the observations were 
made only within that range; (2) it is scarcely arbitrary enough at interme- 
diate times, for the trigonometric terms represent a periodic function; (3) 
it is assumed that the observations are equally spaced. In removing these 
defects, we may as well speak of a general system, exposed to forces of any 
kind whatever. Let all the coordinates, q;, - - -qg, be observed at the instants 
ti, to, - - -t,. We desire to construct a set of functions q(t), go(t), - - - which will 
reduce, at the observation times, to functions which adequately represent the 
observed values of the coordinates and which will be arbitrary at all other 
times. The solution is immediate.® 

We need only write down the polynomials obtained by applying La- 
grange’s method of interpolation, let us say 


qi = filt, C1; C2 aii ), ge = felt, C1, C2, atte >» (8) 


where the c’s are the constant coefficients of the various powers of t. We re- 
place these constants by any functions of ¢t which are entirely arbitrary except 
that they must reduce to the experimental values c,, C2, etc., at the instants of ob- 
servation. \Ve shall call any such function A (¢), or simply A. The Eqs. (8) 
are now replaced by 


| oo fill, Aj, As, ar ), ime falt, A, As, ). (9) 


In these equations we have a faithful picture of experimental facts. It is 
necessary, however, to reconsider the concept of mass, since we naturally 
desire to use it as soon as we consider accelerated motions. In order to obtain 
a definition of inertial mass suitable for use in the finite difference equations 
of dynamics, we must describe a procedure by which a number may be as- 
signed to the mass of a body, without any mention of instantaneous velocity, 
which is a differential coefficient. When we neglect the dependence of mass 
on velocity, the original definition of Newton is quite satisfactory. That is, 
if two freely moving bodies collide head on, their masses are inversely as the 
changes of velocity which they experience. Since the velocity measured is 
Ax/At rather than dx/dt, we see that the value of a mass, resulting from such 
an experiment, can be formulated entirely in terms of finite difference coeffi- 
cients and a constant; namely, the number assigned to the mass of a standard 
body. 

The construction of the Eqs. (8) is especially simple when the time inter- 
val between observations is constant. Just as in the illustration of the uni- 
formly moving particle, we are led to finite difference equations. From the 
standpoint of the principle of real and finite description, these equations, or 
their integrals (in the sense of the calculus of finite differences), expressed in 
the form (9), are the true laws of motion. 

If, for the moment, we suppose there are no experimental errors, and that 
the observations are consistent with Newton’s laws, we can easily see what 
the finite difference equations will be. This supposition is made only for pur- 


5 Boole, reference 3, p. 39. 
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poses of discussion, and implies no recession from the principles laid down in 
section II. The general difference equation for any problem of particle 
dynamics, At being constant, is obtained by noting that mA*x/Af is the time 
mean of the time mean of the force, taken over the interval At. \We may verify 
directly that 


A?x 1 t+At t+At d2x 1 t+At t+At 
m— = - f f m dt dt = —|{ f F ,dt dt, 
Al? At? t t dt? At? t t 


A’x a 
= F,, (10) 
At? 








or, 





ml 


where the bar denotes a time mean over the interval mentioned. This equa- 
tion, and two similar ones, constitute a formulation of the laws of motion, to 
which we should doubtless be led if, in all cases, we were to carry through a 
process like that which led us to Eqs. (3), (4), and (5). In general the right 
member of (10) will depend on Af as well as on ¢, because the upper limits of 
the integrals contain At. The analogue of (10) in curvilinear coordinates is 
very complicated. While the writer has found it possible to make consider- 
able progress in working out a theory for second order difference equations 
analogous to the transformation theory of ordinary dynamics, the results are 
so involved as to be of little interest. Therefore, we content ourselves with 
stating Eq. (10). 


IV. RELATIVITY DyNAMICS 


In discussing special relativity from the standpoint of the principle of real 
and finite description, we meet at once with a kinematic problem; namely, 
the relation between the coordinates in different systems of reference. \When- 
ever we pass to a new system of dynamical variables, the transformation 
equations themselves should be in accord with the principle. In other words, 
they may contain difference coefficients, but they should not contain deriva- 
tives. It is satisfying to note that the Lorentz transformation may be inter- 
preted in such a way that it obeys this criterion. If the moving system of 
reference, x’, y’, 2’, t’, has velocity V( = Bc) parallel to the X-axis of the ob- 
server’s system, x, y, 2, ¢, this transformation takes the form 


Vx 
xs’ = k(x— Vi),y = 97,2) = 32,0 = iC - =), k=1/(1 — B%)'"2, (11) 
c 


On the usual theory, V represents the instantaneous velocity of the moving 
axes, dX /dt, but for the present purpose a different definition is needed. We 
may suppose the two sets of axes are attached to isolated bodies between 
which there is no detectable interaction. If then we find experimentally that 
AX/At is a constant V, we postulate that the Eqs. (11) hold true, and with 
this understanding they contain no derivatives. Quite similarly, c represents 
the quotient of a measured distance travelled by a light beam between two 
stations, and the measured time required for its journey. In accordance with 
the spirit of this paper, c must therefore be expressed as a difference quotient. 
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If we assume, momentarily and for purposes of discussion only, that the 
relativistic differential equations of motion are valid, the difference equations 
of motion for a rapidly moving particle may be obtained very easily. The 
relativistic equations of motion are 


d 
—(mv,) = F,, etc. (12) 
dt 


where m =n/(1—8?)'”? and Bc=v. If the square of the element of arc in a 
four-dimensional space is written as 


ds? = c*dt? — dx? — dy*? — dz’, 
(12) may be rewritten in the form 


d*x —_ dt 
mo— = F,—- (13) 
ds* ds 
Suppose now that the observations are uniformly distributed along the 
world line of the particle, so that we have As = constant, in place of the former 
assumption, Af=constant. We proceed just as we did in obtaining (10) with 
the result 


A*x _ dt 
moc—- = F,—; 
As? ds 


where the right member denotes 


1 s+As stAs dt 
— i) F,— dsds. 
As? J, " ds 


There is a way of expressing the law of light propagation without recourse 
to infinitesimals or to integrals. Since, in the theory of special relativity, the 
law is ds? = df? —dx* —dy* —dz* =0, we replace this by 


As? = c2Af? — Ax? — Ay? — Az? = 0. (15) 


In general relativity also, the law of light propagation may be taken as 
As =0. 

It is much more difficult to formulate an analogue for the geodesic law 
{ds =0, which applies to the free motion of a particle. If we follow the most 
obvious course and assume that the law should take the form 


As = an extremum, (16) 


we at once encounter a question of interpretation. 

On the one hand, we may assume that we are dealing with a continuous 
space-time manifold, even though our measurements are discontinuous in 
character. On this basis As means a finite interval between adjoining observa- 
tions, marked off along a geodesic line, where the word geodesic is taken in 
its ordinary sense. The significance of (16) is then simply that of a theorem 
deducible from Einstein’s law of motion. 
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On the other hand, we may suppose that the space-time manifold is sim- 
ply a discrete aggregate of observed point events, analogous to a two or three 
dimensional network, like a crystal lattice. The problem of defining a geo- 
desic in such a lattice is not an easy one, as Silberstein has pointed out in his 
“Theory of Relativity.” He suggests that the “shortest distance” between 
two points a and d in such a lattice (Fig. 2) should be defined as the smallest 
number of steps which will carry one from a to b. As to the concept of a step, 
it is necessary to think of the lattice points as joined by a network of lines. 
A step is then defined as the transition from a given point to its neighbor on 
one of these lines. In Figs 2a and 2b we have indicated two different methods 
of drawing such line-networks and it is obvious that the shortest distance de- 
pends to a great extent on the nature of these lines. Further, even with such 
a simple point lattice and network of lines as that shown in Fig. 2a, the paths 









































| 
d £ ig > 
c a e f . Iq 
la P 
2a 2b 


Fig. 2. 


which correspond to shortest distances are not unique. Thus, the two paths 
a,c, d---b, anda, c’, d’--- 6 both involve the same number of steps, and 
both may be considered as geodesics between a and b for the particular net- 
work shown. Although the shortest distance has a definite value, the path 
which must be traversed in order to achieve the journey in this shortest dis- 
tance is not uniquely determined, save in special cases. After all, we should 
expect that some feature like this would emerge when we deal with discrete 
quantities. 

In view of the existence of a limit to the accuracy with which a coordinate 
can be measured,® we should expect that the concept of shortest distance has 
an indefinite meaning in the atomic domain. The contribution of the present 
discussion is to show that a similar difference exists even for macroscopic 
measurements, as soon as we abandon the idea of a continuous background 
for the description of events. 


V. Discussion OF EHRENFEST’S THEOREM 
In sections II and III we emphasized that Newton's laws of motion should 
be interpreted statistically—that our knowledge of a single body can never 
be extensive enough to justify us in supposing it follows the laws of motion 
exactly. There is a theorem of wave mechanics discovered by Ehrenfest,’ 


6 Ruark, reference 4. 
7 Ehrenfest, Zeits. f. Physik. 45, 455 (1927); Ruark and Urey, Atoms Molecules and 
Quanta, p. 652, McGraw-Hill, 1930. 
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which throws further light on this situation. Let y denote the wave function 
for a particle, so that *dz is the probability for the particle to have its co- 
ordinates in the element dv of the coordinate space. Then the x-component of 
the center of gravity of the probability distribution is given by X = fy*xpdz, 
provided the wave function is so normalized that fy*~dv =1, and the theorem 
may be stated as follows: 


” a ay” 
aX = ( - — Vwdo. (17) 
OX 


The expressions for Y and Z are similar, and V is the potential energy. 
Now let us consider a wave packet, having a very small extension compared 
with any distance in which the potential energy changes appreciably. In 
evaluating (17), we may use the value of —dV/0x at the centroid of the packet, 
in every element of the integral, and the result is 


” ov oV 
mX = ( — ) [ vrvas = ( _ ~) , (18) 
OvJ ye Ox /x 


On the basis of an outworn interpretation of wave dynamics, the individual 
particle is thought of as a wave packet, and yas the amplitude of the waves. 
If this were correct, Ehrenfest’s theorem would replace the ordinary law of 
motion, and it would appear foolish to use finite difference equations in de- 
scribing the motion of a particle. But, on the probability interpretation, the 
meaning of Ehrenfest’s theorem is this: We consider a great number of 
similar experiments, in each of which a single particle moves along the x-axis, 
and if we measure the x-coordinate of each particle at a definite instant, we 
obtain the distribution of YW* appropriate to that instant. At any future time 
the distribution can be calculated from the wave equations and we have a 
picture of the statistical behavior of the particles at that time. In other words, 
the wave equations enable us only to follow the time changes of the proba- 
bility distribution, and there is room for uncertainty in the motion of the in- 
dividual particles. Only the centroid obeys the approximate law of motion 
(18). Now, from a knowledge of the probability distribution obtained by a 
great number of experiments on individual isolated particles, we can draw 
no definite conclusions as to the motion of the centroid of a great number of 
particles, bound together in a rigid body and set in motion all together. This 
remark should make it clear that there is nothing in Ehrenfest’s theorem 
which contradicts the conclusions of sections IT and ITI. 


VI. DIFFERENCE EQUATIONS IN ELECTRODYNAMICS 


The use of difference equations in electromagnetic theory arises naturally 
from our discussion of mechanics. Details will have to be reserved for a later 
occasion, but the leading ideas will be indicated here. To begin with, the 
electromagnetic field is defined entirely in terms of dynamical measurements. 
For our present purpose, the forces between charged bodies, magnets, and 
conductors carrying currents, are the fundamental things. We may employ 
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a small test body carrying a charge e, to map out a field of electric force. 
Again, if the test body is endowed with a “velocity” As/At, chosen to suit 
our convenience, a field of magnetic force may be studied. Thus we obtain 
the distribution of six quantities, E,, Ey, E., I/,, ,, and H,, which might be 
called components of electric and magnetic intensity if the names were not 
preémpted by analogous entities incapable of being measured, and defined 
in terms of the behavior of infinitesimal test bodies. The quantities E,, - - -JZ, 
are defined by the very processes used in measuring them. There is no implica- 
tion that they obey Maxwell's equations. Indeed, they do not, because, from 
the usual standpoint, they represent average values of the field intensities, 
measured by a finite test body and not by a mathematical abstraction. These 
quantities may well be referred to as the physical field intensities. Their prop- 
erties are simply the properties of numbers which fit a modified Lorentz 
equation, 


1 
eF=e (E+- (vH)). (19) 
c 


In this equation v means As/At rather than ds/dt; e is measured by experi- 
ments whose interpretation requires no use of infinitesimals and eF, is a force 
measured, let us say, by readings on the scale of an actual spring balance, or 
by equating it to the value of mA*x/Af for the test body. The field equations 
obeyed by the physical field intensities may now, in principle, be found by 
trial. They involve the properties of the test body as well as those of the bod- 
ies producing the field—a point which is important in itself. The same con- 
sideration applies in mechanics but was not emphasized in section III, for the 
sake of simplicity. 


VII. ConcLusIon 


In conclusion, this paper serves to introduce a method of talking consist- 
tently about the results of experiments, and attention may be directed again 
to the refinements hinted at in the first foot note. The writer’s aim has 
been to show how physical facts can be formulated by the use of difference 
calculus, that is, essentially with the aid of rational numbers. Considering 
the almost universal use of idealized concepts such as infinitesimals and irra- 
tionals, it seemed worth while to examine the possibilities of proceeding other- 
wise, and the results may be summarized in a broad sense,as follows: a method 
has been outlined by which we can describe experimental results mathe- 
matically, without any implications concerning the state of a system at times 
when it is not being observed, or concerning the causal connections custom- 
arily implied by the use of continuous analysis. 

The writer wishes to thank Doctors Hutchisson, Muskat, Rashevsky and 
Worthing for helpful discussions. 
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Nuclear Spin 


On the basis of measurements by Rasetti' of 
the Raman effect in No, Heitler and Herzberg? 
have suggested that the electron in the nuc- 
leus has lost not only its spin, but also its 
power to determine what the statistics of the 
nucleus shall be. (One could say, perhaps, 
that the electrons in the nucleus behave as if 
they had no spin and obey, therefore, the 
Bose-Einstein statistics, just as photons do.) 
The purpose of the present note is to point out 
that the assumption that the protons alone con- 
tribute to the resultant spin is sufficient to ex- 
plain all the known facts. Furthermore, it may 
be that there is a building-up principle analo- 
gous to that of the outer electrons, although 
the evidence® is insufficient to show exactly 
what this would be. 


Tasie I. Builling-up process. 


————— : 





Protons 


Spins 

Element P Z IT added added 

H 1 1 1/2 

He 4 2 0 3 (+-)- 

Li 7 $ 3/2 3 +++ 

Cc 12 6 0 5 —--—--(+-) 

N 14 7 1 2 ++ 

O 16 8 0 2 -- 

F 9 9 1/2 3 (+-)+ 











The table gives in the last two columns the 
number of protons and spins, respectively, 
which must be added to the element in the 
row above to give the element in the row it- 
self. 

Thus, if we suppose that each additional 
proton contributes a spin of + % to the result- 
ant (which is the simplest assumption), the 
elements may be ordered as shown. Starting 
with He‘, J =0, three parallel spins will produce 
Li’, J=3/2, and five more, two forming a pair, 
can result in C®, J=0. Two additional paral- 
lel spins give N“, J=1; two more, opposite in 
direction to the first two, O"*, 7=0; and three 
more, two forming a pair, can yield J=}4 for 


F!%, The value of nuclear spin for Na® is in 
doubt, but if we assume J =5/2, then, by in- 
terpolation, Ne*® should have 7=1 and Ne”, 
I=2. One would also expect, since J=0 for 
C” and S®*, a value of J= 1/2 for B" and P®. 

The difficulty mentioned by Goudsmit for 
the cadmium spectrum also disappears, if one 
adopts the above hypotheses. The value of 
I=1/2 may be attributed to Cd™ or Cd!", 
but there exist no criteria for choosing be- 
tween the two isotopes. (Perhaps both give a 
hyperfine structure.) 

It seems, therefore, that the experimental 
evidence to date indicates the validity of the 
rule that even numbers of protons imply in- 
tegral values for J, and odd numbers half-in- 
tegral values. This is true for Mn, Br, Pr and 
Bi. 

One would expect aluminium, phosphorus, 
chlorine, potassium, scandium, vanadium, 
cobalt, copper, gallium, arsenic, and other 
elements with isotopes having only odd num- 
bers of protons to have half-integral values 
for the nuclear spin. 

It remains an open question as to whether 
or not the orbital angular momentum of the 
protons plays a réle, and also whether or not 
closed shells of protons exist, as do closed 
shells of extra-nuclear electrons. If there are 
such closed shells of protons, then this might 
provide the explanation of inverted hyperfine 
structure, the inversion being due to the ab- 
sence of a proton from such a shell. 

James H. BARTLETT, JR. 

University of Illinois 

January 15, 1931. 


1F, Rasetti, Proc. Nat. Acad. Sci. 15, 515 
(1929). 

2 Heitler and Herzberg, Naturwiss. 17, 673 
(1929). 

3 R. S. Mulliken, Trans. Faraday Soc. 25, 
634 (1929), 

* Pauling and Goudsmit, 
(1930). 
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The Cathode-Ray Tube in X-Ray Spectroscopy and Quantitative Analysis 


Under this title the Journal of the American 
Chemical Society has published a paper by us 
in the January number. A Coolidge tube was 
used, having a soldered aluminum window 6 
mm in diameter and operated continuously at 
0.2 m.a. and 85 kv D.C. 
made with metal targets placed above the 


Experiments were 


window and cooled with a jet of hydrogen. 
For the recording of the relative intensities of 
the Ae doublet, primary attention was given 
to a method involving the application of an 
ionization chamber and an amplifying set. 
The latter made use of the new FP54 tube 
(Metcalf and Thompson, Phys. Rev. 36, 1489 
(1930)). The smallest measurable grid current 
was 8.7 (10) amps. With a copper target 
and 0.4 mm slits, the peak of the Ka doublet 
gave a net detlection of 90 cm with an accur- 
acy of 0.5 cm ona scale 2.6 m distant, corre- 
sponding to an ionization current of 3 (10)~% 
amp. 

The loss in electron velocity in passing 
through the window agreed with the Thomson 
formula. For molybdenum and iron the inten- 
sity was found to vary as the square of the 
difference between the electron velocity and 
This relationship held 
maximum velocity employed, 


the exciting voltage. 
up to the 
amounting to four-fold the exciting voltage 
for the case of molybdenum and eleven-fold 
for iron. Experiments as yet unreported are 


being made on titanium and calcium and the 
exponent has been found to be 1.3 for the vol- 
tage range extending from ten to sixteen-fold 
the exciting voltage. The exponent accord- 
ingly decreases less rapidly with increase in 
voltage than has been observed in’ similar 
experiments with x-ray tubes. 

Curves showing the relation between chemi- 
cal composition and the relative intensities of 
the unresolved Ka doublets were obtained for 
the lat- 
ter representing the case in which the intensity 


copper-nickel and iron-nickel alloys 


of radiation from one constituent is increased 
as the result of excitation by the other. 

Other experiments were carried out, using 
photographic recording by means of a See- 
mann spectrometer. For a molybdenum tar- 
get, a 1° oscillation of the crystal and a slit 
wilth of 0.05 mm, a density of the Ka; line 
amounting to 0.4 was obtained by an expo- 
minutes. 


sure of sixteen Analyses by this 


method were made on silver-cadmium and 
tin-antimony alloys and the results were con- 


firmed by chemical examination. 


GorTON R. FONDA 
GEORGE B. CoLLINs 


Research Laboratory 
General Electric Company 
Schenectady, New York 
January 16, 1931. 


The Hyperfine Structure of Ionized Lithium 


The hyperfine structure of the low *P state 
of Li* is of the same order of magnitude as 
the triplet separation itself. One might ex- 
pect large deviations from the interval rule 
for hyperfine splitting, corresponding to large 
contributions of one triplet level to the per- 
turbed eigenfunctions of another. Since the 
Li* hyperfine structure has not yet been ex- 
plained conclusively,'! we have estimated the 
effect of these deviations on the calculated 
positions of the components. The resulting 
corrections are, however, small, and are not 
sufficient to give a better agreement than ob- 
tained by Giittinger. 

Our method is simply to consider the transi- 
tion between the extreme cases wherein the 
hyperfine splitting is either small or large com- 


! H. Schiiler, Zeits. f. Physik 66, 431 (1930). 
P, Giittinger, Zeits. f. Physik 64, 749 (1930). 


pared to the triplet separation. The coetii- 
cients of general equations for the levels were 
determined by use of the vector model in the 
extreme cases, and by using the usual identifi- 
cation of the lines for the isotope Lis. The 
nucleus was considered as coupled only to the 
1s electron in the *P. About 6 percent was sub- 
tracted from the coupling constant as deter- 
mined from the 4S, due to the interaction with 
This correction alone 
has some effect on the results of Giittinger, as 
is indicated in the third column of the table. 

In the case of 7=1}, one of the coefficients 
of a cubic equation cannot be determined, and 
is represented by x in Table I. Its order of 
magnitude should be one-tenth cm. Its 
evaluation would be intricate, but no value 


the 2s electron there.* 


2 G. Breit and F. Doermann, Phys. Rev. 36, 
1732 (1930). 
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TABLE I. 


Line® Schiiler Giittinger 
(exp.) 
1 (0.5) 18231.31 ck) ~ 30.38 
2 (1) 30.91 (2) 30.91 
3 (2) 30.26 (3) 30.25 
4c (2) 27.82 
4b (6) 27.72 (5.5) 27.74 
4a (2) 27.62 
5b (12) 27.52 (11) 27.56 
5a (3) 27.42 
6c (35) 27.22 
6b (2) 27.15 (2.5) 27.10 
8c (2) 25.96 
8b (6) 25.91 (2.5) 25.89 
13b (.7) 25.56 
13a (1.5) 25.46 (2) 25.36 
9 (4) 25.25 (5.5) 25.23 
10 (1) 24.93 (2.5) 24.90 


would bring much better agreement than was 


found in Giittinger’s approximation. In the 


3’ The components 7, 11 and 14 are attrib- 
uted to the Li® isotope and have been omitted 
in this table. They are at 26.43, 24.33 and 
29.48. 





i=1} i=} 
Extreme Intermed. Intermed. 
coupling coupling coupling 

31.31 31.31 

30.91 30.91 (1) 30.91 

30.26 30.26 (2) 30.26 

27.79 27.83 —x 

27.88 27.97-—3.1x (5) 27.56 

27.69 27.70—x 

27.58 27.58—x (9) 27.39 

27.48 27 .57—3.1x 

27.29 27 .30—x 

27.14 27 .18—x (1) 26.91 

26.03 26.12 

25.94 25.91—x (1) 25.80 

25.63 25.72 

25.44 25 .44—x (2) 25.51 

25.29 25 .26—x (5) 25.15 

24.98 25.07 (1) 24.86 








| 


case 1=}, our results give agreement just 
slightly worse than those quoted by Schiiler. 


S. GoupsMIT 
D. R. INGLIs 
Department of Physics, 
University of Michigan, 
January 19, 1930. 


An Attempt to Detect Axiality of X-Ray Emission 


J. Stark! has recently reported that fluores- 
cent K-series x-rays are not emitted with 
equal intensity in all directions from certain 
crystals, and that these x-rays are at least par- 
tially plane-polarized. 
dibromoanthracene, a monoclinic crystal of 
needle-like habit with the b-axis parallel to the 
length of the needle, he found about 20 per- 
cent more bromine A-radiation in directions 


In the case of y-y- 


90° from the b-axis than in a direction 10° 
therefrom. The scattering of this radiation by 
aluminium indicated that it was polarized 
with the electric vector parallel to the b-axis. 
These effects are so important from the the- 
oretical side that independent confirmation 
seemed desirable. 
successful attempt to detect differences in the 
intensity of zinc A-radiation emitted in dif- 


This note reports an un- 


ferent crystallographic directions from a zinc 
crystal. 

Zine was chosen because the axial ratio for 
its hexagonal crystals, c/a=1.80, is so much 


1Stark, Ann. d. 
(1930). 


Physik [5] 6, 637-602 


greater than that (1.63) to be expected in a 
close-packed arrangement of spheres that it is 
reasonable to suppose the atomic axes are not 
free to rotate about axes in the (001) plane. 
This is one of the criteria suggested by Stark 
for selecting materials likely to show the new 
effect. His objection to the use of metals, that 
the atoms therein are presumably free to ro- 
tate, does not therefore seem valid with re- 
spect to zinc. 

The primary x-rays used were from a tung- 
sten target tube operated at 55000 volts 
(peak) or from a molybdenum target tube op- 
erated at 44000 volts (peak). The zinc mono- 
crystalline plate, 0.4 cm thick was rotated 
about its normal into various azimuths while 
its secondary radiation in a fixed direction 90° 
from the primary beam was measured by bal- 
ancing its ionization with that of part of the 
primary beam. The absorption coefficient of 
the secondary radiation was found to be that 
of zinc K-radiation, indicating that the scat- 
tered radiation at 90° was negligible in com- 
parison with the fluorescent. The normal to 
the plate was found by Laue and rotation 








330 


photographs to lie at 20° from the hexagonal 
axis of the crystal. so that if @ is the angle be- 
tween the normal and the axis of the principal 
ionization chamber the angles between this 
axis and the hexagonal axis of the crystal vary 
between 0+20° and @—20° during the rota- 
tion of the specimen. For the three settings 
tested, @ was 52°, 62° and 77°. In no case was 
the variation in observed intensity as much 





LETTERS TO THE EDITOR 


as one percent. We conclude that at least as 
regards intensity the effect reported by Stark 
is unimportant in zinc. 


L. A. PARDUE 
L. W. McKEEHAN 
Sloane Physics Laboratory, 
Yale University, 
January 23, 1931. 


Rotational Fine Structure in Raman Spectra 


In recording the rotation spectrum of a 
molecule by scattered light one of the major 
difficulties is the very great intensity of the 
unmodified line as compared with those lines 
arising from changes in rotational energy. 
Long exposures tend to broaden the strong 
line and obliterate the weaker ones which lie 
very close to it. An ingenious arrangement 
employed by Rasetti (Zeits. f. Physik 66, 646 
(1930)) is partially effective in overcoming 
this handicap. He introduced mercury vapor 
into the spectrograph to absorb from the scat- 
tered light those unmodified frequencies which 
arise from the normal state of the atom. The 
observations are thus practically confined to 
\2537 as parent line. Even this expedient is 
not wholly successful, since \2535 is super- 
posed upon the rotation pattern and appears 
very intense in comparison with it. 

Recently Manneback (Zeits. f. Physik 62, 
224 (1930)) has computed the relative intensi- 
ties and polarizations for Raleigh and Raman 
scattering of lines with and without rotational 


displacements. He shows that if the incident 
light is plane polarized and the scattered light 
is observed in a plane normal to the electric 
vector, the rotation lines are almost com- 
pletely depolarized. Since the unmodified line 
is not depolarized, this suggests the introduc- 
tion of a Nicol prism between the camera and 
the scattering chamber. Any convenient 
source and any frequency for the incident line 
could then be used, and not only the parent 
line but also neighboring lines scattered with- 
out change in frequency would be eliminated. 
In all probability the background illumination 
would also be materially lessened. With this 
improvement in contrast considerably smaller 
total exposures should prove adequate, and 
the necessary time might not be much ex- 
tended, even though the illumination would 
be considerably fainter. 
E. F. BARKER 
Department of Physics, 
University of Michigan, 
January 28, 1931. 


Biological Effects of Gamma-Rays 


In connection with the development ir va- 
rious laboratories of x-ray tubes operating at 
voltages of the order of a million volts, the 
question of adequate protection of those en- 
gaged in the work against the very penetrating 
gamma-rays emitted by such tubes has be 
come of considerable importance. It has 
seemed valuable to perform some experiments 
with highly filtered gamma-rays from radium 
in order to ascertain the dangers from expo- 
sure to the penetrating radiation, which may 
be considerable in amount, passing through 
whatever shielding may be used around the 
tube. Estimates doubtless might be made 
from previous work, but for comparison it ap- 
peared desirable to have data on exposures to 
gamma-rays from radium from which all soft 
components were similarly filtered out. 


A group of 63 white rats (with known ante- 
cedents) has been exposed to the radiation 
from 6 grams of radium filtered through one 
mm of platinum, one mm of brass, 16 mm of 
lead, and 5 mm of celluloid at a distance of 41 
mm to the nearest side of the rat. For four ad- 
ditional rats the celluloid was omitted. The 
rats were put in celluloid boxes 47 mm high 
placed above and below the radium box. The 
radium was spread out over a surface of 3 by 
10 cm area to give more uniform exposure. 
With this large source-area and distance, the 
exposure was reasonably similar throughout 
the whole body of the rat. In addition, 16 rats 
were exposed to 2.5 grams of radium at the 
same distance but with the lead filter re- 
moved; for eight of this group the celluloid 
was also omitted. This group of 16 thus re- 
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ceived approximately the same exposure to 
the hard components of the radiation but with 
the addition of that amount of soft gamma- 
radiation which penetrated the platinum and 
brass. The celluloid was used in order to re- 
duce the beta-ray intensity to about the level 
corresponding to the rate of beta-ray produc- 
tion in the tissues. Exposure-times were va- 
ried from one-half minute to 17 hours in a 
roughly geometric progression. The 2.5 grams 
without the lead filter should give an “ery- 
thema dose” in about four hours at the aver- 
age distance of 6 cm according to ordinary 
clinical data. 

All rats died which were exposed six hours 
and longer to the 6 grams of radium with fil- 
ters. With a 16-hour exposure, death ensued 
in three or four days, due at least in part to in- 
jury to the gastro-intestinal tract. With six- 
hour exposures the rats showed similar but 
less severe gastro-intestinal symptoms from 
which all recovered in about three days but 
died one to three weeks later from a progres- 
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protective measure for those working with 
high-voltage tubes. Breeding experiments are 
in progress, and the evidence so far indicates 
that the sterility dose approximates the lethal 
dose under these conditions. Abnormalities 
have to date appeared only in the litter of one 
female fertilized four weeks after a three-hour 
exposure. 

To provide a suitable index of the total ex- 
posure of investigators working with high- 
voltage tubes, a series of Eastman extra fast 
(CC) dental x-ray films was enclosed in a lead 
cassette having walls 6.9 mm thick, and this 
cassette was exposed in the mid-position of 
the rat (6.2 cm to the film) for periods of 10-5 
to 10—' lethal exposure-time as determined for 
the rats. An exposure of 1X10-* lethal dose 
gives a point just at the lower end of the 
straight portion of the density-log exposure- 
curve, that is, 7X10 lethal dose gives a den- 
sity which is easily confused with fog due to 
development, whereas 3X10~ gives a very 
obvious blackening of the film. Densities as 























Fraction of Percent Percent F oe 

lethal transmission transmission Average nny 9! T 

exposure-time incident film emergent film B10 
None 100 100 1.00 0.0 
0.00005 100 100 1.00 0.0 
0.0001 83 97 0.90 0.045 
0.0003 65 83 0.74 0.13 
0.0007 46 o4 0.55 0.26 
0.001 33 45 0.39 0.41 
0.003 9.5 12 0.11 | 0.96 
0.007 1.8 1.8 | 0.018 1.74 
0.01 0.71 0.82 | 0.0075 2.42 
0.05 0.028 0.033 0.0003 | 3.5 
0.10 0 0.013 0 3.7 


.028 


.0002 





sive anemia affecting both red and white 
blood-cells. One of the six rats exposed to 
three hours died from a similar anemia after 
11 days. The lethal exposure lies thus be- 
tween three and six hours, or about four 
hours. Two of four rats exposed 2.5 hours to 
2.5 grams of radium without the lead filter 
died from anemia after 19 and 37 days. The 
lethal exposure for these conditions is thus 
about three hours. 

Blood-counts were made on 66 rats before 
and 0, 1, 3, 7, 14, 21, and 35 days after ex- 
posure. The minimum exposure which showed 
an appreciable effect on the blood-count was 
approximately 20 minutes for the rats exposed 
to the 6 grams with filters. Since a single mas- 
sive dose is the most favorable for detecting a 
change in blood-count and since ten times 
this exposure is above the lethal dose, blood- 
count observations alone are not a satisfactory 


measured on a photoelectric photometer for 
the two films in each placque are given below. 
The placque was exposed as labeled with the 
“black paper toward tube” (radium). 

For these film-exposures 2.5 grams of ra- 
dium were used with the 16-mm lead filter (in 
addition to the cassette), and 12 hours was 
taken as the lethal exposure-time. 

W. G. Wuairman, M.D. 
M. A. Tuve 
School of Hygiene and Public Health 
(W.G.W.), 
Johns Hopkins University, 
3altimore, Maryland. 
Department of Terrestrial Magnetism 
(M.A.T.), 
Carnegie Institution of Washington, 
Washington, D. C., 
January 22, 1931. 
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BOOK REVIEWS 





Vorlesungen uber Wellenmechanik. A. |.\Nbf. 
gesellschaft, Leipzig 1930. Price RM 9.50. 


Pp. 132, figs. 15, Akademische Verlags- 


This small volume represents the lectures delivered at Ohio State University during the 
academic year 1929-30. In general character it is more of a Handbuch than a self-sufficient 
introductory treatise. The number of topics touched upon is surprisingly large, but the treat- 
ment of the individual items is correspondingly shortened. 

After a brief introduction on the general foundations of atomic theory there is a section 
on the dualism of the wave and corpuscle viewpoints, and one on the uncertainty principle. The 
treatment is along the usual lines, but the compactness of the material is so great as to detract 
very considerably from its value. The third section on quantum statistics gives a treatment of 
the Bose-Einstein method but only a mention of that of Fermi. 

The major portion of the book is of course devoted to Schrédinger’s equation and its solu- 
tion. Only the simplest cases are discussed, but no attempt is made to give actual methods of 
solution, results only being stated. One may call attention to the treatment of the Pauli princi- 
ple and the symmetry properties of wave functions, which is rather more carefully stated than 
is usual. The last section on the relativistic wave equation contains a paragraph on Dirac’s 
electron theory and even one on quantum electrodynamics. Perhaps one of the best uses of 
this book will be as a guide in the development of more detailed lecture notes. 


E. L. Hew 


Handbuch der Astrophysik. Band II, Erste Hiilfte. Edited by G. Eneruarp, A. Koxt- 
SCHUTTER AND H. LUDENDORFF. Pp. 430, figs. 51. Julius Springer, Berlin, 1929. Price bound 
RM 69. 

The importance of photometry in astrophysical work is well illustrated by the fact that 
this volume of 430 pages forms but a half of the section of the Zandbuch which is to be devoted 
to this subject. The first and longest chapter deals with photometric theory. The photometry 
of the stars, and of self-luminous bodies in general, is rather summarily treated, for example, 
nothing is said about the measurement of the brightness of the sun relative to the stars, or of 
the comparison between the astronomical standard of light (expressed in stellar magnitudes) 
with the terrestrial standard (meter-candles). On the other hand, the photometry of diffusely 
reflecting bodies—the moon, the planets, Saturn's rings, the zodiacal light, ete., and the extine- 
tion of light in the earth's atmosphere are discussed in great detail, with the inclusion of im- 
portant investigations by the author (Professor E. Schoenberg) some of which are new, while 
others have been difficult of access. This complete monograph, accompanied by extensive tables 
should be of much value to all workers in this field. There are few references, however, to work 
of later date than 1922, and it is to be regretted that no mention is made of Danjon's work on 
the earth-shine upon the moon, or of Hubble's on the illumination of nebulae by neighboring 
stars. 

The chapter on spectro-photometry (by Dr. A. Brill) deals with a subject so new that it is 
largely an account of the methods employed by individual investigators. These are clearly 
presented and the summary is brought up to the end of 1927— including some of the modern 
work on the contours of individual spectral lines. 

A short chapter on colorimetry (by K. Bottlinger) contains an admirable survey of the 
meaning of color-indices and related magnitudes, and a critical discussion of the methods and 
results of observation up to 1926, The final chapter on photoelectric photometry, (by H. Rosen- 
berg) contains a thorough discussion of photoelectric and selenium cells. The latter have been 
practically abandoned in favor of the former, which are of ever-increasing importance. He 
justly remarks that their introduction into observatories has doubtless been delayed by the 
unfamiliarity of astronomers with the details of delicate electrical equipment, and his full de- 
scriptions of instrumental design and wiring should be of real aid here. 

HENRY Norris RUSSELL 
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